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Abstract 

An overview of the physics of spinor and dipolar Bose-Einstein condensates (BECs) is pre- 
sented. Mean-field ground states, Bogoliubov spectra, and many-body ground and excited states 
of spinor BECs are discussed. Properties of spin-polarized dipolar BECs and those of spinor- 
dipolar BECs are reviewed. Some of the unique features of the vortices in spinor BECs, such as 
fractional vortices and non-Abelian vortices, are delineated. The symmetry of the order parame- 
ter is classified using group theory, and various topological excitations are investigated based on 
homotopy theory. Some of the more recent developments such as the Kibble-Zurek mechanism 
of defect formation and the Berezinskii-Kosterlitz-Thouless transition in a spinor BEC are also 
discussed. 
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symbol definition 

^•TfAx^ *■') annihilation operator for atomic pairs in |:F, At) 

A pair amplitude of the spin-singlet pair, {Aoo)o 

ajr s-wave scattering length of the total spin-'F channel 

d unit vector indicating the spin quantization axis 

T, M total spin of two colliding atoms 

/, m atomic spin and its projection 

f (r) = (f I, fy, fj) vector of / X / spin matrices 

Fi'') = 2mn ^In^mn^Pn Spin Operator 

Pir) = Yjnm Vm^mn^n Spin expectation value 

/('■) = Hmn Cn^mn^n Spin expectation value per atom 

i, j, k indices for x, y, and z in the coordinate space 

M atomic mass 

total magnetization of the condensate 

m, n indices for the magnetic sublevel 

n(r) = 2m 'f^m'f^m number density 

p linear Zeeman energy 

q quadratic Zeeman energy 

s = F/\F\ unit vector spin along the spin expectation value 

t, u unit vector satisfying s = i xu 

^m{r) field operator 

il/mir) order parameter 

im = 4'ml yfn normalized spinor order parameter 

H, V, A indices for x, y, and z in the spin space 



Table 1 : Table of symbols 



abbreviation 


definition 


AF 


antiferromagnetic 


BA 


broken axisymmetry 


BEC 


Bose-Einstein condensate 


BN 


biaxial nematic 


C 


cyclic 


CSV 


chiral spin vortex 


DDI 


dipole-dipole interaction 


DW 


density wave 


F 


ferromagnetic 


PL 


fiower 


GP 


Gross-Pitaevskii 


MI 


Mott insulator 


P 


polar 


PCV 


polar-core vortex 


SF 


superfluid 


ss 


supersoUd 


UN 


uniaxial nematic 



Table 2: Table of abbreviations 
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1. Introduction 



Most Bose-Einstein condensates (BECs) of dilute atomic gases realized thus far have internal 
degrees of freedom arising from spins. When a BEC is trapped in a magnetic potential, the spin 
of each atom is oriented along the direction of the local magnetic field. The spin degrees of 
freedom are therefore frozen and the BEC is described by a scalar order parameter. 

On the other hand, when the BEC is confined in an optical trap, the direction of the spin can 
change dynamically due to the interparticle interaction. Consequently, the order parameter of a 
spin-/ BEC has 2/ + 1 components that can vary over space and time, producing a very rich 
variety of spin textures. A BEC with spin degrees of freedom is referred to as a spinor BEC. This 
article provides an overview of the physics of spinor and dipolar BECs. 

Bose-Einstein condensation is a genuinely quantum-mechanical phase transition in that it 
occurs without the help of interaction. However, in the case of spinor BECs, several phases 
are possible below the transition temperature T^, and which phase is realized does depend on 
the nature of the interaction. With the spin and gauge degrees of freedom, the full symmetry 
of the system above is SO(3) x U(l). Below the transition temperature, the symmetry is 
spontaneously broken in several different ways; the number of possible broken symmetry states 
is four for the spin-1 case, at least five for the spin-2 case, and the situation has yet to be fully 
understood for the spin-3 case. This wealth of possible phases makes a spinor BEC a fascinating 
area of research for quantum gases. 

The spin is accompanied by a magnetic moment that gives rise to magnetism. Depending 
on the phases of individual systems, spinor BECs exhibit a rich variety of magnetic phenomena. 
Moreover, because the systems are superfluid, we can expect an interplay between superfluidity 
and magnetism. For example, a ferromagnetic BEC has spin-gauge symmetry; if you rotate the 
spin locally, the system tries to undo it by generating a supercurrent. The basic properties of the 
mean-field spinor condensates are reviewed in Sec.|3] 

The magnetic moment of the atom leads to the magnetic dipole-dipole interaction. Although 
the strength of this interaction is the smallest among all the energy scales involved in the system, 
it plays a pivotal role in producing the spin texture — the spatial variation of the spin direction. 
The magnetic dipole-dipole interaction plays key roles in determining the static and dynamic 
properties of spin-polarized dipolar BECs. It also serves as a key parameter in determining the 
quantum phases in an ultralow magnetic field (below 10 jiG). The physics of dipolar BECs is 
discussed in Sec.|4] 

The Bose-Einstein condensed system responds to an external perturbation in a very unique 
manner Even if the perturbation is weak, the response may be nonperturbative; for example, 
the phonon velocity depends on the scattering length in a nonanalytic manner. The low-lying 
excitations of spinor and dipolar BECs are described by the Bogoliubov theory, as discussed in 
Sec.|5] 

One of the hallmarks of superfluidity manifests itself in its response to an external rotation. 
In a scalar BEC, the system hosts vortices that are characterized by the quantum of circulation, 
K - h/M, where h is the Planck constant and M, the mass of the atom. The origin of this 
quantization is the single-valuedness of the order parameter As mentioned above, however, in 
the spinor BEC, the gauge degree of freedom is coupled with the spin degrees of freedom. This 
spin-gauge coupling gives rise to some of the unique features of the spinor BEC. For example, 
the fundamental unit of circulation can be a rational fraction of k, and when two vortices collide, 
they may not reconnect unlike the case of the scalar BEC. Vortices of spinor BECs are discussed 
in Sec.|6] 
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The direction of the spin can vary rather flexibly over space and time. Nonetheless, the global 
configuration of the spin texture must satisfy certain topological constraints. This is because the 
order parameter in each phase belongs to a particular order parameter manifold whose symmetry 
leads to a conserved quantity called a topological charge. Such a topological constraint deter- 
mines the nature of topological excitations such as line defects, point defects, skyrmions, and 
knots. The order parameter can be classified using group theory, and possible topological excita- 
tions are determined based on homotopy theory. Elements of group theory and homotopy theory 
with applications to spinor BECs are reviewed in Sec. |7] 

Efibrts are underway to magnetically shield the system so that the ambient magnetic field 
is reduced to below 10 /iG. In such an ultralow magnetic field, the spin degeneracy, which is 
usually lifted by Zeeman eff'ects, sets in significantly, and we have to confront the question of 
how the degeneracy can be lifted by minute spin-exchange interactions. It is in this region where 
the many-body spin correlations dramatically alter the ground states of spinor BECs and their 
response to an external magnetic field. The present understanding of this problem is described in 
Sec.g] 

Spinor BECs offer topics that are of broad interest beyond ultracold atomic gases. Prime 
examples are the dynamics of a quenched BEC and the Berezinskii-Kosteritz-Thouless transition 
in two dimensions. These subjects are discussed in Sec.|9] 

It seems almost certain that we will face more interesting questions than have been answered 
so far. In fact, several new possibilities are emerging both experimentally and theoretically: mag- 
netic crystalization and non-Abelian vortices, to cite only a few examples. Moreover, under an 
ambient magnetic field below 10 yuG, not only do manybody spin correlations become significant 
but the magnetic dipole-dipole interaction also plays an equally important role. We therefore 
have good reasons to expect that new symmetry breakings will occur in this regime. In view 
of these ongoing developments, it now seems appropriate to consolidate the knowledge that has 
been accumulated over the past few years. In this paper, we provide an overview of the basics 
and recent developments on the physics of spinor and dipolar BECs. Some of the remaining and 
emerging problems are listed in Sec.fTOl 

One of the major topics that is not treated in this review is fictitious spin systems such as 
a binary mixture of hyperfine spin states. Such systems do not possess rotational symmetry in 
space nor is the fictitious spin associated with the magnetic moment. However, they have some 
intrinsic interest such as interlaced vortex lattices and vortex molecules caused by the Josephson 
coupling. A comprehensive review of this subject is given in Ref. lit]. 

This paper is organized as follows. Section |2] describes the fundamental Hamiltonian of 
the spinor BEC. Section [3] develops the mean-field theory of spinor condensates and discusses 
the ground-state properties of the spin-1, 2, and 3 BECs. Section |4] provides an overview of 
the dipole BEC and the spinor-dipolar BEC. Section |5] develops the Bogoliubov theory of the 
spinor and spinor-dipolar BEC. Section|6]discusses various types of vortices that can be created 
in spinor BECs. In particular, fractional vortices and non-Abelian vortices are discussed. Sec- 
tion Q examines the topological aspects of spinor BECs. The symmetry of the order parameter 
is classified based on group theory, and possible topological excitations are investigated using 
homotopy theory. Section|8]reviews the many-body aspects of spinor BECs. Section|9]discusses 
the Kibble-Zurek mechanism and the Berezinskii-Kosterlitz-Thouless transition in spinor BECs. 
SectionfTOlsummai'izes the main results of this paper and discusses possible future developments. 
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2. General Theory 

2.1. Single-particle Hamiltonian 

The fundamental Hamiltonian of a spinor BEC can be constructed quite generally based on 
the symmetry argument. We consider a system of identical bosons with mass M and spin / that 
are described by the field operators ^, where m - f,f- 1, ■ ■ ■ , -/ denotes the magnetic 
quantum number. The field operators are assumed to satisfy the canonical commutation relations 



where 5„,„ is the Kronecker delta. 

The noninteracting part of the Hamiltonian comprises the kinetic term, one-body potential 
U{r), and linear and quadratic Zeeman terms: 



where f, is the z-component of the spin matrix whose matrix elements are given by (fz)mn — niSmn, 
and hence, (ff )„,„ - rrp-S,„,i', and p - -gpeB is the product of the Lande hyperfine g-factor g, 
Bohr magneton ps = ehjlm^ (nie is the electron mass, and e > is the elementary charge), and 
external magnetic field B that is assumed to be applied in the z-direction. The coefficient q of the 
quadratic Zeeman term is calculated using the second-order perturbation theory to be 

q = — , (3) 

where E^f - E,^ - Ei is the hyperfine splitting and it is given by the difference between the 
initial (Ei) and intermediate (E^) energies. For example, for the case of **^Rb, g - -1/2 and 
Ehf ^6.8 GHz for the f - I hyperfine manifold and g - 1/2 and E\if ^ -6.8 GHz for the f - 2 
hyperfine manifold. The value of q can be tuned in the negative direction by using a linearly 
polarized microwave field due to the AC Stark shift iQilstl- 

2.2. Interaction Hamiltonian 

2.2.1. Symmetry considerations and irreducible operators 

Next, we consider the interaction Hamiltonian. We discuss the case of the magnetic dipole- 
dipole interaction later as it does not conserve the total spin of the system, and hence, it requires 
a separate treatment. Because the system is very dilute, we consider only the binary interaction. 
Upon the exchange of two identical particles of spin /, the many-body wave function changes 
by the phase factor (-1)^^. By the same operation, the spin and orbital parts of the wave function 
change by (-1)^^^^ and (-1)"^, respectively, where ;F is the total spin of the two particles and X, 
is the relative orbital angular momentum between them. To be consistent, (-1)^^ must be equal 
to (-1)'^'^^^ X (-1)"^; hence, (-1)^^"^ = 1. Thus, T + £. must be even, regardless of the statistics 
of the particles. In the following, we consider only the s-wave scattering (X = 0), and therefore, 
the total spin of two interacting particles must be even. The interaction Hamiltonian of the binary 
collision can therefore be written as 



[0'„(r),(Al(r')] =5,„„(5(r-r'), 
[0',„(r),(A„(r')] = [^I(r),^I(r')] =0, 



(1) 




(2) 





(4) 



r=0,2,-,2/ 
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where V^"^^ is the interaction Hamiltonian between two bosons whose total spin is 9^. 

The interaction V*^' in Eq. (HI can be constructed from the irreducible operator Afr^(r, r') 
that annihilates a pair of bosons at positions r and r' . The irreducible operator is related to a pair 
of field operators via the Clebsch-Gordan coefficients {T, M\f, niuf, 1112} as follows: 

/ 

ArM(r,r')^ ^ <5", M|/, mi; /,m2>-Am,W<A,«. (»•')• (5) 

m] ,m2=-/ 

Because the boson field operators ij/„,M) and iff„^(r') commute, Ajryv((r) vanishes identically for 
odd f. In fact, when - 1 and f - I, the Clebsch-Gordan coefficient is given by 

(-1)'-'"' 

(1, A1|l,mi; l,m2) = — 6,n^+m2,M 

V2 

X [<5A(,l(^mi,l + ^mi,o) + ^A(,0 '«1 - 5A1,-l(<Jm,,0 + ^m,,-l)]- (6) 

Substituting this into Eq. (|5]l, we find that Aiyvi(r) = 0. Similarly, we can show that Ayr^(r) = 
if ;F is odd. Incidentally, for the case of fermions, we obtain Ajrj^(r) = for even ;F because the 
fermion field operators anticommute. 

Because the interaction is scalar, it must take the following form: 

y(^' = ^ fdrf dr'v^'^\r,r') ^ A^^^ir,r')ArM(r,r'), (7) 

^ ^ M=-r 

where v^'^\r,r') describes the dependence of the interaction on the positions of the particles. 
When the system is dilute and the range of interaction is negligible as compared to the inter- 
particle spacing, the interaction kernel v^^* may be approximated by the delta function with an 
effective coupling constant gjr: 

v'-'^^(r,r')^gr6(r-r'), (8) 
where gf- is related to the i-wave scattering length of the total spin-;F channel, oyr, as 

Sr - (9) 

2.2.2. Operator relations 

It follows from the completeness relation 

Y^\r,M}{r,M\ = i (10) 



that 



where 



r 



Z Z A^rM^r,r')ArM(r,r') h(r)h(r') (11) 
r=o,2,-,2f M=-r 



f 

h(r) = 2 4'lir)kr(r) (12) 

m=-f 



is the total density operator and :: denotes normal ordering that places annihilation operators to 
the right of creation operators. 

Another useful relation can be derived from the composition law of the angular momentum: 



F^-F2^\ [(Fi + F2f -F\- F]\ = if L - /(/ + 1), 



(13) 



where Ftot = Fi + F2 is the total spin angular momentum vector. Operating this on the com- 
pleteness relation (flOl ). we have 



1 



:nr + i)-f{f + i) 



\r,M){r,M\. 



(14) 



Applying this formula to the case of f - 1, we obtain 

2 



Fi-F2^ |2,A1)<2,M| -2|0,0>(0,0| = i -3|0,0><0,0|, 



(15) 



M=-2 



where we have used the completeness relation (fTOl i to derive the second equality. It follows from 
this that the following operator identity holds: 



3Ajo(r,r')Aoo(r,r')+ : Fir) ■ Fir') :=: «(r)«(r') :, 
where F - iF^, Fy, F.) is the spin vector operator defined by 



(16) 



F^ir)^ ^ ii^)mn^mir)^nir) in^x,y,z). 

m,n=-f 

Here, (f^),„„ iju - x,y, z) are the (m, n)-components of spin matrices f^. 
For the case of / = 2, we have 

f I ■ f 2 = -6|0, 0)<0, 0| - 3 ^ |2, M><2, M| + 4 ^ |4, M}{4, M\. 



M 



M 



(17) 



(18) 



The last term can be eliminated by using the completeness relation (fTOl i. giving the following 
operator identity: 

2 

^ : Fir) ■ Fir') : + J] ^m^''' '•'M2M(r, /) + yAjo(r, r')Aooir, r') = ^ : n(r)n(r') : . (19) 

A1=-2 

Similarly, for the case of / - 3, we obtain 
21 



1 

IT 



— : Fir) ■ Fir') : +— Ajo(r, r')Aoo(r, r') 



+ TT Z AlM('-^'-'^^^Mir,r') + J] A;^(r, r')A4M(r, r') = ^ : «(r)«(r') : . (20) 

M=-A 
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2.2.3. Interaction Hamiltonian of spin- 1 BEC 

We derive the interaction Hamiltonians for the / = 1 case. For f - I, the total spin T of two 
colliding bosons must be or 2, and the corresponding interaction Hamiltonian (|7]i gives 



y(0) 
y(2) 



drAUr)AQQ{r) 



(21) 



^ dr [: 



A2^(r)A2Ai(r) 



gi 
2 



ir) : -A;o(r)Aoo(r) 



(22) 



where Ayr/^ir) = Ayrj^(r,r), and Eq. (fTTT i and Aiyvt('') = are used to derive the last equality. 
Combining Eqs. (1211 1 and (l22l i. we obtain 



I : n2(^) : +^^L_IiA;,o(r)Aoo(r) 



(23) 



Equation ( l23T l is the interaction Hamiltonian of the spin-1 BEC. We use the operator identity (fTST i 
to rewrite Eq. ( |23] ) as 



y 



where 



i J flfr[co:nV):+ci : f V) :] . 



^0 + 2^2 82 - go 

Co = , ci = — - — . 



(24) 



(25) 



In Eq. (l23T l. Aqo is the spin-singlet pair operator The explicit form of this operator can be 
found using the Clebsch-Gordan coefficient for f - M - 0: 



{0,0\f,muf,m2}^6, 



{-l)f- 



mi +1112,0 



V27TT 



(26) 



Substituting this into Eq. (|5]l, we obtain 

1 ^ 
V2/ + 1 „,=-/■ 



(27) 



2.2.4. Interaction Hamiltonian of spin-2 BEC 

For f -2,T must be 0, 2, or 4, and the interaction Hamiltonian d?) for ;F = 4 gives 



AM 



(r)A^M(r) 



y(^) . 'i\drY,A 

^ M=-4 

r \ ^ 

= Yjdr: h\r) : -A;,(r)Aoo(r) - ^ A^^(r)A2A,(r) 



M=-2 



(28) 
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where Eq. (fTTl l and Aim - ^3M - are used to obtain the last equality. Summing Eqs. (|2TI) . 
(I22I 1. and ( |28] |. we obtain 



y 



y(0) + yP) + y(4) 



-I 



dr 



— : n^(r) : +- 
2 ^ ' 2 



-A5o(r)Aoo(r) + ^^-^ ^ A^^(r)A2Ai(r) 



A1=-2 



We can eliminate the last term by using the operator identity ( fT9l ). obtaining 
V = ^ Jflfr [co : n2(r) : +ci : F\r) : +C2A^^{r)A^{r)\ ., 



where 



4^2 + 3^4 gA - g2 IgO - 10^2 + 3^4 

Co = > c\ - — - — , C2 - 
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(29) 



(30) 



(31) 



We use the same notations cq and ci for both spin-1 and spin-2 cases because there is no fear of 
confusion. 



2.2.5. Interaction Hamiltonian ofspin-3 BEC 

In a similar manner, the interaction Hamiltonian of the spin-3 spinor BEC is given by 

V = ]- \dr cq: h\r) : +ci : F^(r) : +C2Ajo(r)Aoo(r) + C3 ^ A^^(r)A2At('-) 



A1=-2 



, (32) 



where 



9g4 + 2g6 g6-g4 1 IgO - 21^4 + 10^6 1 1^2 - 18^4 + Tg6 ,,,, 
Co - 7-; , Ci = 7-; , C2 = 7-; , C3 = 7-; . (33) 
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3. Mean-Field Theory of Spinor Condensates 

3.1. Number-conserving theory 

The mean-field theory is usually obtained by replacing the field operator with its expectation 
value (i^m). This recipe, though widely used and technically convenient, has one conceptual 
difficulty; that is, it breaks the global U(l) gauge invariance, which implies that the number of 
atoms is not conserved. However, in reality, the number of atoms is strictly conserved, as are the 
baryon (proton and neutron) and lepton (electron) numbers. In fact, it is possible to construct the 
mean-field theory without breaking the U( 1 ) gauge symmetry 

m. 

To construct a number-conserving mean-field theory, we first expand the field operator in 
terms of a complete orthonormal set of basis functions {(pmiir)}' 

>l^m{r) = ^ a„i(p„i(r) (w = /, / - 1 , ■ ■ ■ , -/), (34) 

where (pmi{r) describes the basis function for the magnetic quantum number m and spatial mode 
/, and the corresponding annihilation operators that satisfy the canonical commutation 

relations 

[ami,alj] = 6mnSij, [ami,a„j] = [«',„■, a ^] = 0. (35) 
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The basis functions are assumed to satisfy the orthonormality conditions 



and the completeness relation 



(36) 



(37) 



Then, the field operators ( l34l i satisfy the field commutation relations ([T]i. 

In the mean-field approximation, it is assumed that all Bose-condensed bosons occupy a 
single spatial mode, say / = 0, and a single spin state specified by a linear superposition of 
magnetic sublevels. Then, the state vector is given by 



10 = 



/ 

\m=-f 



mO 



I vac), 



(38) 



where |vac) denotes the particle vacuum and ^,„'s are assumed to satisfy the following normal- 
ization condition: 



2 1^-1' = 1- 

m=-/ 



It is straightforward to show that 



<'Am(r)>o = <C(r))o = 0, 
<iAl(r)0'„(r')>o = r„ir)Mr'), 



{il,UrWn(r')Ur")^cir"'))Q = 1 - - C('-)'A:('-')'Ai('-")<AKr"'), 



N 



where <■ ■ ■ ■)() = {(\---- \0 and 



(39) 



(40) 
(41) 

(42) 



(43) 



As shown in Eq. ( l40t . the expectation value of the field operator vanishes as it should in a number- 
conserving theory. Nevertheless, all experimentally observable physical quantities, which are 
expressed in terms of the correlation function of the field operators such as (HTt and (l42T i. can 
have nonzero values. These values agree with those obtained by the U(l) symmetry-breaking 
approach within the factor of 1 /N. 

3.2. Mean-field theory of spin- 1 BECs 

We use Eqs. (HTt and (l42l i to evaluate the expectation value of the Hamiltonian H - Hq -\- V 
over the state (l38T l with / = 1, where Hq and V are given in Eqs. Q and (l24l l. By ignoring the 
terms of the order of 1 IN in Eq. ( l42l l. we obtain 



= (H) 



= N''] X 

^ ^ m=-\ 



' 2M 



H- lJ(r) — pm -H qm 
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-Am + + yl^-n, (44) 



where 



n{r) = <«(r)>o = ^ |<Am(r)|2 



(45) 



is the particle density and F - (F^, Fy, F~) is the spin density vector defined by 



(46) 



m,«— -1 



Because the spin-1 matrices are given by 



1 



0] 

1 

oj 



fv 



I 



ro 


-1 


0) 




'1 





0] 


1 





-1 













lo 


1 


oj 




lo 





-ij 



the components of the spin vector F are exphcitly written as 

Fx = 1^*1^0 + i/'(*)(i/'i + lA-i) + lA-iiAol , 

Fy = [-lAlfAo + |/'S(^/'l - (A-i) + Ci'Ao] , 



(47) 

(48) 

(49) 
(50) 



From the energy functional in Eq. i44i . we find that co must be nonnegative; otherwise, the 
system collapses. This condition also guarantees that the Bogoliubov excitation energy corre- 
sponding to density fluctuations is real, as shown in Sec. |5] 

It is easy to identify the ground-state magnetism of a uniform system when the Zeeman terms 
are negligible {p - q - 0). From the energy functional (l44l l. one can see that the ground state 
is ferromagnetic {\F\ - n) when ci < and antiferromagnetic or polar (\F\ - 0) when ci > 0. 
Here, the order parameter of the antiferromagnetic phase [6] is given by Vn/2(1, 0, 1)^ and that 
of the polar state 0] is ^/n(0, 1,0)^. The former state is obtained by rotating the latter about 
the v-axis by 7r/2, and therefore, these two states are degenerate in the absence of the magnetic 
field. Moreover, all states obtained by rotating these states are degenerate. In the presence of an 
external magnetic field, however, the ground state becomes more complicated due to the linear 
and quadratic Zeeman effects, as shown below. 

The time evolution of the mean field is governed by 



ih 



dt 6rm(r) 
Substituting Eq. (l44l i into the right-hand side gives 



(51) 



ih 



dt 



' 2M 



+ U(r) - pm + qm 



+ CQmf/,„ + ci ^ ^' ■ fm/7<A« (m^ 1, 0, -1), 



(52) 
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which are the multicomponent Gross-Pitaevskii equations (GPEs) that describe the mean-field 
properties of spin-1 Bose-Einstein condensates llalTl]. In a stationary state, we substitute 



i/r,„(r, t) = (A,„(r) e >' 
in Eq. ( |52] |. where /i is the chemical potential, to obtain 



2v72 



' 2M 



+ U(r) — pm + qm 



(Am + conif/,,, -H Ci ^ f ■ f „,„((/„ = iui//„ 



n=-l 



Writing down the three components m - 1,0,-1 explicitly, we obtain 



2v72 



'1m 



+ U(r) - p + q + C()« + ciFr - fj 



V2 



V2 



2M 



+ U(r) + CQti - ^ 



lAo + -^F^i/z^i = 0, 
V2 



2M 



+ U(r) + p + q + cqh - c\F,- fi 



iff-i = 0, 



(53) 
(54) 

(55) 
(56) 
(57) 



where F-t = F^ + iFy. By solving the set of equations (|55]|-(|57]|. we can investigate the properties 
of the ground and excited states of the spin-1 BEC. 

Note that because the system is suspended in a vacuum chamber, the projected spin angular 
momentum in the direction of the external magnetic field is conserved for a long time (> 1 s) ['^. 
Because the atomic interaction given by Eq. (|7]) conserves the total spin of two colliding atoms, 
the spin dynamics obeying the GPEs conserves the total magnetization. Therefore, when we 
search for the ground state for a given magnetization M, - J drF-(r), we need to replace E in 
Eq. (ISTl i with E - AM,, where A is a Lagrange multiplier Then, p in Eqs. (I55ll-(l57]i is replaced 
with p - p + A, which is determined as a function of M,. 

Here, we solve Eqs. (I55ll-(l57li and investigate the ground-state phase diagram Jgl [l^ in a 
parameter space of {q,p). We assume a uniform system for a given number density n. The 
ground state of a uniform system is obtained if we set the kinetic energy and U (r) to be zero. We 
also note that cqu can be absorbed in /i by defining p = fi - con. Furthermore, we may assume 
Fy - without loss of generality, provided that the system is symmetric about the z-axis. Then, 
by performing a gauge transformation to make (^o real, all lAi.o -i become real. (After obtaining 
the real solution, we can recover the general complex order parameter by performing a global 
gauge transformation e'^" and spin rotation e^-^'^~ about the z-axis.) Then, Eqs. (|55] |- (|57] | reduce 
to 



(-p + q + ciF- - p)i//i + ci((Ai + lA-OiAo = 0, 
[/i-Ci((Ai +lA-i)^]l/ro = 0, 
ci((/fi + + (p + q - ciF, - =0. 

The energy per particle is given by 



1 

e - — 

n 



^{-pm + qm^)\i//„ 

m 
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|' + ^co«2 + ici|f|2 



(58) 
(59) 
(60) 



(61) 



From Eq. $59[ . we have either (i) i/tq = or (ii) jl = ci(i/ri + i/f-i)^. In case (i), we have three 
stationary states: 



T 11 

I : (e'-^' yfn, 0, 0) , ei = -p + q + ■:jCQn + ^'^yn, 

T 11 

II : (0, 0, e'-^ ' Vn) , £\\-p + q+ -con + -c\n. 



III 



em = + z-con - 



2cin ' 



(62) 
(63) 

(64) 



where we have recovered the phases of the order parameter with;^'±i = X() + Xz- States I and II 
are fully polarized (F, - +n), and therefore, ferromagnetic, whereas state III has a longitudinal 
magnetization that depends on p because - p/ci. State III is often referred to as an antifer- 
romagnetic state. In case (ii), we may solve Eqs. ( ISSl l and ( l60l l. together with the normalization 
condition 'Zjm=-i llAmP - n, to obtain the following two stationary states; one, referred to as a 
polar state, is given by 



IV: (0,6'^° V^,0)^, eiv = ^co«. 



and it has no magnetization in any direction; the other is given by 
V:(A±i 



2q ^ 



—p^ + q^ + 2c\nq 



Iciq 



{q^ - p^)(-p-^ - q^ + Icinq) 



{-p^- + q^ + 2qcinf 1 
ev = 2 + o'^o"' 



(65) 

(66) 

(67) 
(68) 



and it exists only when the argument of the square root for every iff„, is nonnegative. State V 
occurs for ci < (ferromagnetic case). When 



p^ - q^ +2\ci\nq > and q^ > p^. 



(69) 



all three components and i^o in Eqs. ( |66] | and ( |67] | are nonzero and the magnetization tilts 
against the external magnetic field: 



P(P - q +2q\ci\n) 
2\ci\q^ 



F+ = e'^ 



^(q^-p^){(p^ + 2\c,\nqy-cf 
2\ci\q^ 



(70) 



The polar angle [arctan(|F+|/Fj)] is determined by the interaction and the magnetic field, whereas 
the azimuthal angle {x^ is spontaneously chosen in each realization of the system. Because the 
axial symmetry with respect to the applied magnetic field axis is broken, this phase is referred to 
as the broken-axisymmetry phase [10]. 

Comparing the energies of Eqs. (I62b-(l68b. we obtain the phase diagram in a parameter space 
of (_q,p), as shown in Fig.[T| j^. 
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(a) ci>0 (b) <:i=0 (c) c,<0 




Figure 1: Phase diagram of spin-1 BEC for (a) C[ > 0, (b) C| = 0, and (c) cj < 0. The rotation symmetry about the 
magnetic field is broken in the shaded region. 





phase 


order parameter ifi^ 


F, 


fi^ H- cqh 


e — e - ^cqh 


(I) 


F 


(e'>' Vn, 0, 0) 


n 


—p + q + cin 


-p + q + jC\n 


(11) 


F 


(0, 0, e*-' V") 


-n 


p + q + cin 


p + q + jC\n 


(III) 


AF 




1 n-p/ci p 
\l 2 ^ c, 


q 


^ LC.\n 


(IV) 


P 













(V) 


BA 


Eqs. (HSll, dSlI 


p(-p- +q-+2qcin) 


q + 2ci« - ^ 


{-p^^f\-lqc\nf- 



Table 3: Possible ground-state phases of a spin-1 BEC, where F, AF, P, and BA denote ferromagnetic, antiferromagnetic, 
polar, and broken-axisymmetry phases, respectively. 



Note that the rotation symmetry about the magnetic-field axis is broken in state III, despite 
the fact that the magnetization is parallel to the magnetic field in this phase. To understand the 
underlying physics, we introduce the spin nematic tensor defined by 

1 /f„fv + fvf„\ 

The nematic tensor for the ferromagnetic (I, II) and polar (IV) phases is independent of 
Xz, and given as 





/ 1 

2 





0^ 




'1 












1 

2 










1 







lo 





Ij 




lo 





oj 



respectively. Here, Nxsi - Nyy indicates that the order parameter has rotational symmetry about 



16 



the z-axis. On the other hand, the nematic tensor for state III for;^^, = is given by 



(l+ff) 0^ 
5(1-0') 
1) 



a= 1 1 - 



P 

c\n 



(73) 



Here, Nxx + N^y indicates the anisotropy on the xy plane in the spin space. This is the physical 
origin of the axisymmetry breaking of state III. In fact, it can be shown that off-diagonal elements 
such as Nxy appear when;t', + 0. 

In general, the symmetry of the order parameter can be visualized by drawing the wave 
function in spin space. Because an integer-spin state can be described in terms of the spherical 
harmonics Y'"{s), where s is a unit vector in spin space, the order parameter for a spin-/ system 
can be expressed in terms of a complex wave function given by 



*I'(^) = 2«AmJ7'(^). 



(74) 



Figure |2] shows the order parameter for states I-V. From Fig.|2] one can see that the order pa- 
rameters for states III and V have the same morphology and that the magnetization is parallel to 
the magnetic field in one phase (III) but not in the other (V). Additional symmetry properties are 
discussed in Sec. 17.1.11 



IV |FM V 




phase 27r 



Figure 2: Order parameter for spin-1 stationary states. Surface plots of |'F(i)|- are shown, wliere the gray scale on the 
surface represents arg'P(s). 



3.3. Mean-field theory of spin-2 BECs 

We take the expectation value of the Hamiltonian H - Hq + V over the state in Eq. ( l38T l with 
/ = 2, where Hq and V are given in Eqs. (|2]l and ( l30b : 



EW = {H}o 



' 2M 



+ Uir) - pm + qm 



2 ^l|r-ii2 ^^|^|2 

+ -n^ + -IFP + -|A|- 



(75) 
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Here, the spin-2 matrices are given by 





1 










°1 




fO 


—i 








0^ 


1 







Vf 










i 
















/ 












, fv = 























Vf 





1 















— / 


.0 










1 


0. 




.0 








/ 


0. 


rl 











> 

















1 



























































-1 

















.0 











-2j 















(76) 



As compared to the energy functional (l44l i of the spin-1 BEC, the new term C2|Ap/2 appears in 
Eq. dTST l. where 



A = <Aoo(r)>o 



V5 



(77) 



is the amphtude of the spin-singlet pair. 

For the spin-1 case, the spin-singlet amplitude is uniquely related to the magnetization by 
relation ( fTSI l. and only one of them appears in the energy functional. However, for the spin-2 
case, they can change independently, giving rise to new ground-state phases. This can be best 
illustrated when the system is uniform. In this case, the kinetic-energy term can be ignored and 
the Hartree term con^/2 is constant; therefore, the ground state is determined by the remaining 
terms in Eq. ( iTST l. 

Let us first consider the ground state in the zero magnetic field. The phase diagram in this 
case is determined by the last two terms in Eq. ( iTSl l. as shown in Fig. |3] When c 1 < and C2 > 0, 
the energy of the system is lowered as magnetization increases. Therefore, the ground state is 
ferromagnetic, and the order parameter is given by 1^^''=™ - ( ^n, 0, 0, 0, 0)^ or its rotated state in 
the spin space. 

When c\ > and C2 < 0, the energy of the system is lowered as the amplitude of the 
spin-singlet pair increases. Therefore, the ground state is uniaxial nematic (UN) or biaxial 
nematic (BN), and the order parameters are given by ^""'^^ - (0, 0, yfn, 0, 0)^ (uniaxial) or 
i^*"™ - Vn/2(1, 0, 0, 0, 1)^ (biaxial). While these states are often referred to as polar or antifer- 
romagnetic flllillstl, respectively, here, we adopt the liquid-crystal terminology which implies 
the symmetry of the order parameter (see the insets in Fig. O. In this review, we refer to C4 
phase discussed below as antiferromagnetic, in order to distinguish it from the UN state. The 
order parameter -\/n/2{Q, 1,0,1,0) also represents the BN state, which is obtained by rotating 
iff^"^'^ as /e-'fv'r/2g-!f:ff/4^bmx^ while cannot be obtaind by rotating tfi^^'^'^ since the symmeti-y 
of the order parameters are different. Note that the symmetries of the UN and BN states differ; 
however, the two states are degenerate at the mean-field level. Moreover, the superposition of 
these two states Vn/2(cos ^, 0, V2 sin ^, 0, cos ^) (^ is an arbitrary real) is also degenerate. It has 
been shown that zero-point fluctuations lift the degeneracy lfl2.[l31. 

When both c\ and C2 are positive, neither the ferromagnetic nor UN/BN state is energetically 
favorable and frustration arises, resulting in a new phase referred to as the cyclic phase I.4„.ll.i5ll. 
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The order parameter is given 

^cyclic ^ V^/2(l,0,/>^,0, If, which possesses tetrahedral 
symmetry, as shown in the inset of Fig. [3] In the many -body ground state of the cycHc state, three 
bosons form a spin-singlet trimer and the boson trimers undergo Bose-Einstein condensation, as 
shown in Sec.[8l 

Can 

cyclic 




Figure 3: Phase diagram of the spin-2 BEC at zero magnetic field. In each phase, the profile of the order parameter 
'!'(.«) = 2,„ lAmF^C^Ois shown. 

The time-dependent GPEs for the spin-2 case can be obtained by substituting Eq. (iTSl l in 
Eq. (EB: 



±2 



dt 



_d\lJ±] 



ih 



#0 

dt 



' 2M 



+ U(r) + 2p +4q + cqh ± 2c\F~ - fi 



C2 



V5 



(78) 



2y72 



2M 

•ci 



+ U(r) + p + q + Con + ciF, - /i 



2M 



V6 

—F^il/Q + F^il/^2 



+ U(r) + Con - lu 



- ^^<^;i> (79) 

i/fO + ^ci{F^ifri+F^ifr_i) + ^A4rl, (80) 
^ V5 
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where 

= Fl = 2(^;iAi +'A*i'A-2) + V6(iAl<Ao + (81) 
F, = 2 (|iA2l' - l^-if) + l^il' - (82) 

The time-independent GPEs can be obtained by substituting i/',„(r, f) - ^„{r)e^'^'^^ in Eqs. (l78ll- 
(l80l l. If we assume that the system is uniform (i.e., U{r) - 0), then Eqs. (fTSll-dSOll reduce to 

(4q + 2yo - ft) tp2 + aiA* 2 = -T-'Ai . (83) 
a*iA2 + (4q - 2yo - P) fA'i = -T-iA-p (84) 
7+1/^2 + (^ + 70 - /i) lAi - aiA-i = — Y "^-'^'^^ ^^^^ 

-a \pi+{q-yQ- p) tff_i + y+^_2 = — Y ^""^O' ^^^'^ 
. V6 

/ii/^o - fliAo = — ('"+i/'i + y-^-\}, (87) 

where fi = - cqii, y± = ciF±, 70 = ciF,, 70 = 70 - p, and a = C2A/ yfS. We may use the degree 
of gauge transformation (i.e., the U(l) phase) to make i/'o real. Furthermore, because the physical 
properties of the system are invariant under a rotation about the quantization axis (i.e., the z-axis), 
we may choose the coordinate system to make Fy - such that 7+ = 7_ = 7. A mean-field 
ground state can be obtained as a solution of Eqs. (|83]|-(|87]| under such a simplification. 

Here, we consider the case of 7 = 0, i.e., the case of no transverse magnetization. Then, 
Eqs. (I83]l-(l83 reduce to 



{4-q + 2yo -Jl)4i2 + atff*_2 


= 0, 


(88) 


0*^2 + m - 2-70 - //) ^'^2 


= 0, 


(89) 


(q + yo-p)i//i -ai^*_^ 


= 0, 


(90) 


-a*ipi + (q - yo - p) 


= 0, 


(91) 


(fi - a) ipa 


= 0. 


(92) 



In this case, the free energy per particle is given by 

, 2 

e = - y i-pm + qm^) + ^ + ^fI + ^\A\\ (93) 
n ^—'^ 1 In In 

m=-2 

Because Eqs. (I88l l- (l92l l are decoupled into three parts, the solutions can be classified accord- 
ing to the determinant of the coeflicient matrix of Eqs. (l88b and ( l89b . 

D2={4q-iif-Ayl-\a\\ (94) 

and the determinant of the coefficient matrix of Eqs. ( |90l l and ( |9l1 l. 

D,={q-jlf-yl-\a\\ (95) 
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• If Di and D2 + 0, then we have i^ii - ip^i - and tlii - >^-2 - 0, so that the solution 
is the UN state, that is, the order parameter is given by 



UN: (0,0, Vn, 0,0)^. 



(96) 



The chemical potential is determined from ( |92] l to be /i = ^C2n and F, = 0. The energy 
per particle is found from Eq. ( l93T l to be e = j(co + ^C2)n. 

If Di =0 and D2 0, 4^2 - lA-a = 0. From the assumption that there is no transverse 
magnetization, iA±i 9^ leads to iJ/q = 0. Clearly, 



F1+ : (0, e'^' Vn, 0, 0, 0)^ 



(97) 



and 



Fl_ : (0, 0, 0, e'^-' VH, 0)^ 



(98) 



are the solutions of Eqs. ( l90b and (IgTT i with - n,fi - -p + q + citi and e - -p + q + 
|(co + ci)«, and F~ - -n,ji - p + q + c\n and e - p + q + ^(co + ci)n, respectively. When 
\p\ < \c\ - \c\ \n, a solution with nonzero ili-tx exists, which is given by 

^7■ 



C2 : 



0,e'^ 



n - F^ 



,0 



(99) 



where F 

\ 

5 



C2). The chemical potential and the energy per particle are Ji 



P/(ci - ^ 

q + ^Qn and e = q + i(co + ^02)11 - p^/[2(ci - jC2)n]. At p = 0, this state continuously 
becoms the BN state. 



• If Di and D2 = 0, i/'i = i/'-i = 0. Moreover, if jl a, we have 1//0 = 0. In a manner 
similar to the previous case, Eqs. dSST l and (|89] l have three solutions. Two of them are 
ferromagnetic and they are given by 



F2+ : (e*2 V", 0, 0, 0, 0)^ 



(100) 



with F~ - 2n,jj. - -2p + 4q + Ac\n and e = -2p + 4q + j(co + 4ci)n, and by 



F2_ : (0, 0, 0, 0, e'^-- yfnf 



(101) 



with - -2n,fi - 2p+Aq+Ac\n and e = 2p+Aq + ^{cQ+'Xc\)n. When \p\ < |2ci - ^C2\n 
the other solution is given by 



C4 : 



'.-V^(-f),o,o,o,.-7i(«-^) 



(102) 



where F, - p/(ci - ^q),/* - 4q + \c2n and e-Aq+ i(co + \c2)n - p^/[2(ci - ^C2)«]. 
This state becomes the BN state at p = 0. 
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• If Di and D2 - and // = a, we have i/^i = ij/^i = 0; however, i/'±2 and i/'o can be 
nonzero. They obey \iff2\^ + \4'-2?' + |iAol^ = « and \ip2\^ - |(A-2p = j-F;,, and hence. 



C2' 



« + FJ2 - 



0, I/'O, 0, e'-^ 



« - FJ2 - ifrl 



(103) 



Because a = (c2/5) 



,'ur2+^-2) 



= /i is real, we have ;i'2 + X-2 - 



^(„_^0)2_^2/4 + ^2 

or n. Here, /i and F, are determined from Eqs. ( l88l l and ( [89l ) to be /< = 2^ - f^/{2q) and 
= (yo + p)lc\, respectively, where yo is a real solution of the following equation: 



To + PTo + 4?[? + 2ci(n - i/'o)]ro + 4^^^ = 0, 
and i/'o is determined so as to minimize the energy per particle given by 



€ = 4-q 



1 - 

n 



1 fo^ - /^^ Q 
H — Con H 1 

2 2cin lOn 



F^ 



(104) 



(105) 



In particular, when p = and |<7| > 2ci(n - i/'y), yo is determined to be zero, leading to 
F- = 0. Moreover, when C2 > 0, Eq. (11051 ) is minimized ?Ax2+X-2 - ^ and i^^ = jn+5q/c2 
to be e = 2q' + jCqh - lOq^/ (c2n). The corresponding order parameter is given by 



D2' : 



V" - 10^/c2 ^ /n - 10^/c2 ^ . V" - 10^/c2 



-,o, 



, 0, /e 



(106) 



where we set x±2 - \ ± X- This state continuously becomes the cyclic state at ^ = 0. On 
the other hand when C2 < 0, Eq. (1105b is minimized at;^'2 +X-2 - and the order parameter 
is obtained to be UN for ^ > and BN for ^ < 0. 



If D] = D2 = 0, we solve the simultaneous equations D\ -0 and D2 = to obtain 



2q 



(107) 



When q 0, p a, we find from Eq. ^ that i/'o = 0. Then, y = 2 (i/'ji/'i + lA* jiA-a)- We 
may use Eqs. ([88]l-(l9TJ and (1107b to show that 



y = 2cii/'2i/'i 1 



70 + 3g 



2cii/'*ii/'_2 1 



yo - 3g 

To + 3^ 



= 0. 



(108) 



Thus, we must, in general, have i/'2i/'i - i/'-ii/'-2 - 0. To be consistent with Eqs. (l88ll-(l9Tb. 
we find that either i/'i - i/'_2 = or i/'2 = i/'-i = should hold. In the former case, the 
order parameter is given by 



C34 



n + F- 



, 0, 0, e'^ 



2fi-F, 



(109) 
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with F, - (p - q)/ci and e - 2q - |co« - (p - qf'lilcxn). In the latter case, the solution 
is given by 



C3- : 



2n + F, 



'-,0,0, e'^ 



n — Fy 



(110) 



with Fj. - {p + q)lc\ and e - 2q - ^cqii - {p + qf' l(2c\n). At zero magnetic field 
(i.e., p - q - 0), these two states become the cyclic state, which are related to ^if'^>''='"= = 
V^/2(l,0,iV2,0, l)^as 



1 ^ 





V2 
) 

\ 
V2 



1 ) 



- -I exp 



(■■'4) 



f V - fy 1 

exp I -/ — — arccos — = I il/ 
V2 V3, 



cyclic 



3;r\ 



4 / 



f , + f. 



1 



= -/ exp if J. — exp -/ - — — arccos — — ifi' 



V2 



V3 



f cyclic 



(111) 



(112) 



The above results are summarized in Table|4] By comparing the energy of the obtained state, 
we find the ground state phase diagram of a spin-2 BEC, as shown in Fig.|4]for the special case 
of p = and ^ < 0. 
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state order parameter ifi^ 



fJ- CQll 



7 Con 



F2+ (e'>- V^,0,0,0,0) 
F2_ (0, 0, 0, 0, e*-^ V") 



2n -2p + Aq + Ac\n -2p + 4q + 2cin 
—2n 2p + Aq + Ac\n 2p+Aq + 2c\n 



F1+ (0, e'^' Vn- 0, 0, 0) 
Fl- (0, 0, 0, e'^-' Vn> 0) 



n —p + q + c\n 
-n p + q + c\n 



-p + q + ^c\n 
p + q + jCin 



UN (0, 0, e'>» V", 0, 0) 



C4 ie'x^- yl'^, 0, 0, 0, 



4? + ^C2n - 



2(C]-C2/20)« 



C3+ 0, 0, e'>-> yf^, 

C3 (o, e'>' -y/^^, 0, 0, e'>-2 



C2 (o, .^2^, 0, e'>-' ^/s^, 



2^ 



Cl 



2q ■ 
2q- 



(p-lf 
2c\n 



P 



c,-o/5 



q + T^C2« - 



2(ci-e2/5)n 



C2' 



4^(1 - ^) + 



10« 



7(«-<Ao)'-f^' 
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Table 4: Possible states for the ground state of a spin-2 BEC. Each states is associated with a particular symmetry, as identified in Sec. 17.1.1] In the C2' state, yo is a real solution 
of Eq. )lQ4k and ^2 +X-2 ~ 2^0 = ^ for Q > and ^2 +X-2 ~ 2^0 = for C2 < 0. At p = the C4 and C2 states continuously become the BN state, and the C3± state becomes 
the cyclic state. For C2 > 0, the C2' state becomes the D2' state at p = and the cyclic state at p = ij' = 0. 




Figure 4: Ground-state phase diagram of a spin-2 BEC for p = and ^ < 0, where F2+ states and C3+ states are 
degenerate in the region of F2 and C3, respectively. The M state and C state in Ref. [141 correspond to C3 and D2' in 
this phase diagram, respectively. 



4. Dipolar BEC 



In this section, we consider the dipole-dipole interaction (DDI) that, unlike the i-wave con- 
tact interaction, is long-range and anisotropic. The DDI between alkali atoms is a thousand times 
smaller than the short-range interaction for the background scattering length. However, recent 
experimental developments, such as the realization of a BEC of ^tQjjjg |fl5i[l6i[l7i[l8i[l9i 
l2ofl . a tec hniq ue for creating cold molecular gases 12 ll l22i [2311 . and high-accuracy measure- 
ments I24, 25, 2^, have enabled us to study several consequences of the DDI. Here, we briefly 
introduce the general properties of the DDI in cold gases (Sec. 14.1b . and discuss two impor- 
tant cases: spin-polarized dipolar gases (Sec. 14.2b and spinor dipolar gases (Sec. 14.3b . Excellent 
reviews on dipolar gases have been published by Baranov 1271 12811 and by Lahaye et al. I29ll . 



4.1. Dipole-dipole interaction 

4.1.1. Scattering properties of dipole interaction 

The dipole-dipole interaction (DDI) between two dipole moments with relative position r 
[see Fig.|5](a)] is given by 



Vdd('-) = Cdd 



Cdd'^di^Qnv(r)d2y, 



(113) 



where di^2 are unit vectors indicating the directions of dipole moments, r = |r|, r = r/r, and 
Quvif) (M,y - x,y,z) is a rank-2 traceless symmetric tensor that can be expressed in terms of 

25 



rank-2 spherical harmonics F™(r) as 



(114) 



l6n]_ 



' ^Y^(r)-Y^(r)-Y2Hr) 
iY^ir) - iY^Hr) 
Yhr)-Y-\r) 



iYlir) - iY^Hr) Yl(r) - Y^Hr) 

^Y^(r) + Y^ir) + Y-^r) -iY\ (r) - iY^' (r) 

-2^F0(r) , 



-iYhr) - iY-\r) 



(115) 



In a special case in which the dipole moments are polarized under an external field, the DDI 
takes a simple form: 



1 - 3 cos^ e 



I6n Y^if) 

— Cdd — 
5 



(116) 



where 6 is the angle between the polarization direction and the relative position r [see Fig.|3b)]. 
Because Vdd is negative (positive) for - Q (n/l), the DDI interaction favors a "head-to-tail" 
configuration [Fig.|5](c)] as opposed to a "side-by-side" one [Fig.|5](d)]. 



(a) 



(c) 




(d) ^ f 



t ; 



Figure 5: Interaction between two dipole moments for (a) an unpolarized case and (b) a polarized case. The dipole- 
dipole interaction is most attractive for (c) the head-to-tail configuration and most repulsive for (d) the side-by-side 
configuration. The bottom figures below (c) and (d) show characteristic spin textures. 

The long-range character (~ 1/r^) of the DDI becomes prominent in low-energy scattering. 
When a potential decreases as 1/r" for large r, the pha se shift di behaves in the zero-energy 
limit as A:^'^' if Z < (n - 3)/2 and as k"^^ otherwise f3(J\. For a potential with n > 3, such as 
van der Waals potential (« - 6), the main contribution at A: —> is an i-wave channel So ~ k, 
and therefore, the potential can be described with a single parameter, namely, with the i-wave 
scattering length. On the other hand, when n - 3, all partial waves contribute to the scattering 
as 6i ~ k, and therefore, the scattering process is no longer described with a single parameter 
such as the .s-wave scattering length. Moreover, the anisotropy of the DDI induces coupling 
between partial waves. Because the dipole interaction is described by rank-2 spherical harmonics 
and it therefore has t/-wave symmetry, a partial wave of angular momentum [l, m) is coupled 
to (l,m + Anif) and {I + 2,m + Am/), where Am/ is the change in the projected total spin of 
colliding particles and Am/ = for the polarized case. Although there is no diagonal term for 
the I = channel, this coupling induces the effective potential for I = that behaves as r^^ at 
large distances. Hence, the dipole interaction generates a short-range interaction as well. Yi and 
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You 1 31 , m have shown that away from the shape resonance the effective interaction 

between two dipoles can be described with the following pseudo-potential (see also 1,35.., 36] ): 



VefF(r) = TT-^^ir-) + ^dd(r), (117) 

M 

where a{d) depends effectively on the strength of the dipole moment. 

4.1.2. Dipolar systems 

There are several systems that undergo long-range and anisotropic DDI. The coefficient 
in Eq. ( II 13b is given by c^d = d^/i^-J^eo) and Cdd - i-iod^/iAn) for electric and magnetic dipole 
moments, respectively, where eo and /io are the dielectric constant and magnetic permeability 
of vacuum, respectively, and d is the magnitude of the dipole moment. The typical order of 
magnitude of the electric dipole moment is a product of the elementary charge e and the Bohr 
radius oq, whereas the magnetic dipole moment of an atom is of the order of the Bohr magneton 
fiB ■ Therefore, the ratio of the energy of the magnetic DDI to that of the electric DDI is given by 



2 

a^~10-\ (118) 



(eflo)/e() 

where a - HAnefic) ^ 1/137 is the fine structure constant. 

Due to the strong electric dipole moments, gases of polar molecules are ideal candidates for 
dipole-dominant systems. The permanent dipole moments of the lowest ''^2^ states of heteronu- 



clear alkali dimers are calculated to be of the order of 1 Debye (^ 3.335 x 10"^" Cm) Il37ll . 
However, because the ground state of the molecule is rotationally symmetric and the dipole mo- 
ment averages to zero, we need to apply a strong electric field, typically of the order of 10"^ 
V/cm, to polarize the dipole moments. Several resear ch g roups are now trying to cool molecules 



using techniques such as buffer-gas cooling 112 IL l38l 13911 and Stark deceleration (see 114011 for a 



review). Another way to achieve quantum degeneracy of heteronuclear molecules is to create 
weakly bound Feshbach molecules from cold atoms and to transfer the molecules to the tightly 



bound ro vibrational ground state 11221 |23L 14111 



The effects of the DDI have recently been observed in several systems. The Stuttgart group 
have achieved a BEC of chromium (^^Cr) L15.1 . The ^^Cr atom has a magnetic dipole moment of 
6/iB, which is six times larger than that of alkali atoms, and therefore, the DDI is 36 times larger 
than that of alkali atoms. Moreover, the s-wave scattering length of ^^Cr can be decreased by 
means of Feshbach resonance |42]. Thus, the Stuttgart group realized a dipole-dominant BEC 
and observed magnetostriction [16,117, 18] and anisotropic collapse lfl9l,l20ll . 



A dipole-dominant BEC, albeit a very weak one, has also been realized in a BEC of ^Li 
by almost quenching the 5- wave contact interaction using Feshbach resonance \23}- Fattori et 
al. observed the decoherence of the interferometer of ^^K due to the dipole interaction ll24ll . 
The Berkeley group observed an evidence of the crystallization of spin domains in a BEC of 
8^Rb [25J. 

4.1.3. Tuning of dipole-dipole interaction 

Tunability of the DDI is crucial for systematically investigating the properties of the dipolar 
BEC. The strength of the electric dipole moments can be manipulated by an external electric 
field with the magnetic dipole moment kept constant. Giovanazzi et al. [43] proposed a method 
to control the sign of the dipole interaction as well as its strength by using a rotating field. This 
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method is applicable to both electric and magnetic dipole interactions. Consider a magnetic 
dipole moment under a rotating magnetic field described by 



B{t) = B 



'sin cos Q.t^ 
sin if sin Q.t 
cos ip 



(119) 



where if is the angle between the rotating axis and the dipole moments (see Fig.|6l). The angular 
frequency Q. is set to be much smaller than the Larmor frequency so that the magnetic dipole 
moment adiabatically follows the direction of B. When Q is much larger than the trap frequency, 
the effective dipole interaction is time-averaged over InjQ.. Substituting d - BI\B\ in Eq. ( I113l l 
and taking its time-average, we obtain the following effective interaction: 



<Vdd> - Cdd 



1-3 cos- 9 3 cos^ 



(120) 



The last factor in the right-hand side of ( I120l i can take a value between -1/2 to 1. Hence, one 
can change the sign of the dipole interaction as well as its strength by varying the angle ip. 




Figure 6: Two dipoles adiabatically following the external magnetic field )119) that rotates with angular frequency n. 
The sign and strength of the dipole-dipole interaction jl20t can be controlled by varying the angle ip. 



4.1.4. Numerical method 

In the mean-field treatment, we need to calculate the effective potential induced by the distri- 
bution of the dipole moments d{r), which includes an integral of the form 

Yj r d'r'Q,,{r-r')d,{r'). (121) 

While this integral converges, the 1 /r divergence of the integrand makes it difficult to numerically 
carry out the integration in real space. We therefore use the convolution theorem to rewrite ( I121l i 
as H 

e^.v(^)5;(^)e*'-, (122) 

V k 
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where Q^v{k) = f dre^''^ '' Qi^y(r) and dy(k) = J dre^''^'''dy(r). Here, the Fourier components 
of the kernel Q/jy can be found as follows. First, we note that the spherical harmonics is the 
eigenstate of the three-dimensional Fourier transformation, i.e.. 



■'''■''Ypir)dQr = 4-m'j,(kr)Y'i"ik), (123) 



where k = \k\,k - k/k, ji(x) is the Ith spherical Bessel function and dD.r denotes the integration 
of the angular part of r. Because the radial integration is then carried out as 



Jo 



— r— r dr = -, (124) 
3 



we obtain 



Q^yik) = j dre'''--Q^y{r) = -y (V - ^K^v) ■ (125) 
For spin-polarized gases, we only require 

Ayr 

ez.(^) = -y(l-3cos2^?,), (126) 

where 0k is the angle between the polarization direction and the momentum k. 

In the numerical calculation, we evaluate dy{k) from dy{r) by means of a standard fast Fourier 
transform algorithm — we multiply dy(k) by Qijy{k), sum the result over v, and perform the in- 
verse Fourier transform. The result is not sensitive to the grid size Ar in the coordinate space as 
long as dy{l/Ar) is negligible. However, a small numerical error arises, depending on the system 
size R, because the discrete Fourier transform presupposes a 3D periodic lattice with a unit cell 
of linear dimension R. If the Fourier transform is performed in a cubic region ofR\ the accuracy 
improves if we introduce an infrared cutoff to calculate the Fourier transform of the integration 
kernel as Qpk) = J^^^^^^^drQ.yirye-''-'- M- 

4.2. Spin-polarized Dipolar BEC 

When the dipole moment of every atom is polarized in the direction of an external field, say 
in the z-direction, the condensate is described with a single component order parameter tff(r) that 
obeys the GP equation 



at 



+ y,„p(r) + ||(A(r, tt + 1>dd(r, f) 



Hr,t), (127) 



where ytrap('') is a trapping potential, g = Anti^a/M with a being the scattering length for this 
component, and 

Odd(r,f) = Cdd r dr'^^^^^^\^{r',t)f (128) 

is a mean-field effective potential due to the DDI. Here, we assume a repulsive short-range inter- 
action fl > so that the condensate is always stable in the absence of the DDI, whereas Cdd can 
be either positive or negative due to the tunabihty discussed in Sec. 14.1.31 
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4.2.7. Equilibrium shape and instability 

Due to the anisotropy of the interaction, the equiUbrium shape of dipolar BECs is highly 
nontrivial 47, 45]. The dipole interaction works attractively along the direction of dipole 



moments and repulsively in the perpendicular direction (see Fig. |5j- The attractive part makes 
a dipole-dominant BEC unstable in a homogeneous system. This instability can be qualitatively 
understood from the Bogoliubov spectrum in a homogeneous system [46.1 : 

w = yjet [ek + 2ng {1 + fddO cos2 0^ - 1)}], (129) 
where n is the number density of atoms and 

fdd = — (130) 
3^ 

is a parameter that expresses the relative strength of the DDI against the short-range interaction. 



Equation ( I129l l shows that for edd > 1, the BEC becomes dynamically unstable because o) is 
imaginary for wave numbers below — '\j4-gnM(edti - l)/h. 

When the system is confined in a trapping potential, the quantum pressure due to the confine- 
ment can stabilize the condensate against the attractive interaction. Santos et al. [48] numerically 
investigated the equilibrium shapes of pure dipolar gases with a = 0. They considered an ax- 
isymmetric trap with trap frequencies ll),- and in the radial and axial directions, respectively. 
When the trapping potential is prolate, i.e., when the aspect ratio / - ^jcorloj? > 1, the dipole 
interaction is always attractive. Therefore, as in the case of attractive s-wave interaction, there 
exists a critical number of atoms above which the BEC becomes unstable. On the other hand, if 
the trap potential is sufliciently oblate so that I < I* ^ 0.43 and if the DDI is not too strong (see 
the last paragraph of this subsection), the dipole interaction is always positive and the BEC is 
stable. In the region of /* < / < 1, the aspect ratio of the condensate increases with the number of 
atoms, and finally, it becomes unstable. Yi and You |31, 32] investigated the stability including 
the effect of repulsive short-range interaction using Gaussian variational ansatz, and obtained 



results consistent with those of Ref. H48I1 . 

In Refs. 1 32, 44], the low-energy excitation spectrum is calculated using time-dependent 
variational methods. It is shown that the instabilities of dipolar BECs in prolate (/ > 1 .29) and 
oblate (/*</< 0.75) traps originate from dynamical instabilities of the breathing and quadrupole 
modes, respectively. For the intermediate region, the lowest excitation frequency corresponds to 
a breathing mode far below the above mentioned critical number of atoms, and to a quadrupolar 
one near the critical value. 

In the Thomas-Fermi limit, ^« » 1, the GP equation (1127b has an exact solution ll49ll . Using 
the relationship (6^v - ^rftry)/r^ - -V^Vv(l/r) - 47r5^v5(r)/3, the mean-field dipolar potential in 
Eq. (1128b can be rewritten as 



'l'dd('') = -Cdd 



. . 47r , ■ 
— 0(r) + —n(r) 



(131) 



where <f)(r) = J dr'nir')/\r - r'\, and therefore, (p(r) obeys Poisson's equation, V^^ - -Ann{r). 
Because n{r) is parabolic in the Thomas-Fermi limit, Poisson's equation is satisfied hy (p - 
ao+aix^ +a2y^ +a-iZ?' +a4X^y^+aiy^z?' +a(,z?' +a-ix^ +a%y^ +a()Z' , and hence. Odd is also quadratic. 



Therefore, for a harmonically trapped BEC, all terms appearing in the GP equation ( 1127b are 
quadratic, and as in the simple s-wave case, an inverted parabola for the density profile gives 
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a self-consistent solution of the dipolar hydrodynamic equations. (For an analytic form of the 
solution, see Refs. 1,49. ,50.1 .) 

Eberlein et al. 05011 pointed out that the potential seen by atoms exhibits a local minimum 
immediately outside the condensate, i.e., the sum of the trapping and dipole potentials is locally 
smaller than the chemical potential, causing an instability that brings atoms out from the con- 
densate to fill the dip in the potential. The presence of this instability is in agreement with the 
biconcave shape of a condensate pointed out by Ronen et al. [sT] . By developing a high-precision 
code for dipolar BECs, Ronen et al. L51,1 have shown that the same instability causes a dipolar 
condensate to become unstable in the limit of strong dipole interaction, even though the trap is 
strongly oblate, in disagreement with the result of Refs. 1148113 1113211 . Interestingly, in a strongly 
oblate trap, the wave function of a BEC immediately before collapse becomes biconcave, with 
its maximum density away from the center of the gas. They have also shown that the instability 
of a strongly oblate dipolar BEC is analogous to the "roton-maxon" instability reported for 2D 
dipolar gases (see Sec. 14.2.3b . 



4.2.2. Dipolar collapse 

The Stuttgart group observed the magnetostriction of a ^^Cr BEC, where edd - 0.16 for 
the background scattering length. By reducing the s-wave scattering length using a Feshbach 
resonance, Lahaye et al. observed a change in the aspect ratio of the condensate jlSll . 

The same group also observed a collapse dynamics due to the DDI ll20ll . In the experiment, 
the i-wave scattering length is adiabatically decreased to 30% of the background value, and it is 
then suddenly decreased below the critical value of the collapse. After a certain hold time, the 
BEC is released from the trap and the TOE image is taken after expansion. 

The collapse dynamics proceeds as follows. After a sudden decrease in the scattering length, 
atoms begin to gather at the center of the trap due to the attractive interaction; with an increase in 
the atomic density at the trap center, the three-body recombination losses become predominant; 
as a consequence of the atom loss, the attractive interaction weakens, and the atoms are ejected 
outward due to the quantum pressure. 

Figure [T] (a) shows the measured TOE image (upper panels) together with the results of nu- 
merical simulation (lower ones). We have numerically solved the GP equation that includes the 
three-body loss using the loss rate coefficient that is determined to best fit the measured loss 
curve (Fig. 3 of Ref. [20]). The excellent agreement between the experiment and the theory 
demonstrates the validity of the mean-field description even for dipole-dominant BECs. More- 
over, the numerical simulation revealed that the cloverleaf pattern in the TOF image arises due to 
the creation of a pair of vortex rings, as indicated in Fig.|7](b). When collapse occurs, the atoms 
were ejected in the xy plane (vertical direction to the dipole moments), whereas atoms still flow 
inward along the z direction, giving rise to the circulation. The velocity field shown in Fig. |2](c) 
clearly shows the d-wave nature of the collapse dynamics. 

4.2.3. Roton-maxon excitation 

Santos et al. [52,] have shown that the Bogoliubov excitation spectrum exhibits roton-maxon 
behavior, i.e., the excitation energy has a local maximum and minimum as a function of the 
momentum q, in a system that is harmonically confined in the direction of the dipole moments 
(i.e., the z direction) and free in the x and y directions. If in-plane momenta q are much smaller 
than the inverse size L of the condensate in the z direction, excitations have a 2D character. 
Because the dipoles are perpendicular to the plane of the trap, particles efficiently repel each 
other and the in-plane excitations are phonons. Then, the DDI increases the sound velocity. For 
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Figure 7: (a) Series of absorption images of tlie collapsing condensates for different values of hold time (top) and the 
coiTesponding results of the numerical simulations obtained without adjustable parameters (bottom), (b) Iso-density 
surface of an in-trap condensate. The locations of the topological defects are indicated by the red rings, (c) Velocity field 
of the atomic flow in the A" = plane. Color represents the velocity (red is faster). Reprinted from Ref. l20h . 



q » 1/L, excitations acquire a 3D character and the interparticle repulsion is reduced due to 
the attractive force in the z direction. This decreases the excitation energy with an increase in q. 
When Edd > 1, the excitation energy reaches a minimum (roton) and then begins increasing, and 
the nature of the excitations continuously become single-particle like. As the dipole interaction 
becomes stronger, the energy at the roton minimum decreases and then it reaches zero, leading 
to an instability. The 3D character is essential for the appearance of the roton minimum, which 
does not appear in the quasi-2D system where the confining potential in the z direction is strong 
and the BEC has no degrees of freedom in this direction [53]. The roton-maxon spectrum is also 
predicted for a ID system with laser-induced DDI 1,54,1 . 

4.2.4. Two-dimensional solitons 

It is known that a nonlinear Schrodinger equation with short-range interactions admits soliton 
solutions in one dimension, but it does not in higher dimensions. However, Pedri and Santos llssll 
showed that with a nonlocal DDI, a two-dimensional BEC can have a stable soliton. 

We consider a dipolar BEC polarized in the z direction, and assume that the system is confined 
only in the z direction by a harmonic potential ytrap('') = jMcj^z^, and that it is free in the xy 
plane. Then, the DDI is isotropic in the xy plane. We consider the following Gaussian variational 
ansatz 
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+ z^ 



(132) 



where /q = -\/h/(M(Aj^) and Lp and Ly are dimensionless variational parameters that characterize 
the widths in the xy plane and the z direction, respectively. Using this ansatz, the mean-field 
energy is evaluated to give 

1 1 



ho), 

E{Lp,L,)^-^ 



^ ^ 'l+WW], (133) 



where k = Lp/L^ is the aspect ratio and 



2/c2 + 1 3^2 j arctan V/f^ - 1, k> I 

^W-— - (,._l)^ {iln(^). .<!■ 

For fixed in the absence of the DDI (edd - 0), E(Lp) monotonically increases or decreases 
as a function of Lp, resulting in collapse or expansion. However, in the presence of the DDI, 
because /(k) is a monotonic function of k with /(O) = -1 and /(/c ^ oo) = 2, E(Lp) may have a 
minimum. When the trapping potential is strong and = 1, for simplicity, the energy minimum 
appears for 

fdd < 1 + -x/^-^ < -2edd, (135) 
V 2 aN 



where a is the s-wave scattering length and A^, the total number of atoms (see Ref. II55I1 ). This 
condition holds only for e^^ < 0, which can be achieved using a rotational field. On the other 
hand, Tikhonenkov et al. [56,1 have shown that when dipole moments are polarized in the 2D 
plane, stable 2D soUton waves can be generated without tuning the dipole interaction. In this 
case, the soliton is anisotropic and elongated along the direction of polarization. 

Because a 2D bright soliton in a dipolar gas is stable, the roton instability discussed in the 



preceding subsection does not cause a collapse of the BEC but creates 2D solitons 115711 . If the 
dipole moments are polarized perpendicular to the 2D plane, these solitons are stable as long as 
the gas remains 2D. However, if the dipoles are parallel to the 2D plane, the (anisotropic) solitons 
may become unstable even in 2D if the number of particles per soUton exceeds a critical value. 

4.2.5. Supersolid 

The Bose-Hubbard model with long-range interactions exhibits a rich variety of phases such 
as the density wave (DW), supersoHd (SS), and superfluid (SF) phases mSlialMl. While SF 
and SS phases have a nonzero superfluid density, DW and SS phases have a nontrivial crystalline 
order in which the particle density modulates with a periodicity that is different from that of 
the external potential. Thus, in the SS phase, both diagonal and off-diagonal long-range orders 
coexist |62, 63, 64]. 

A dipolar gas in an optical lattice is an ideal system to realize such exotic phases, because the 
interaction parameters can be experimentally controlled. Goral et al. |65] showed that DW and 
SS phases, as well as SF and Mott insulator (MI) phases, can be accomplished in dipolar gases in 
a 2D optical lattice, and Yi et al. ll66ll investigated detailed phase diagrams in 2D and 3D optical 
lattices. A dipolar gas in an optical lattice can be described with an extended Bose-Hubbard 
model given by 

//ebh = -tY^bftj + hytd + Yj «'<«'■ - 1) + 5 Z ^iM^i - 1) + 5 Z U^hihj, (136) 
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where b, is the annihilation operator of a particle at the lattice site i; hi - bft),, the correspond- 
ing particle-number operator; and t, the hopping matrix element between the nearest neigh- 
bors, and Uq > is the on-site Hubbard repulsion due to the i-wave scattering. Here, t and 
U can be expressed in terms of the Wannier functions w(r - r,) of the lowest energy band as 
U = Anah^lM /dr|w(r)|^and t = jdrw*{r - r/)[-^V2 + Voir)]w{r - rj), where Vq is the 
optical lattice potential ll67ll . The last two terms in Eq. ( I136l l describe the on-side and inter-site 
DDIs with coupling parameters given by 

Uu^CM I dr I dr'\w(r-rir^-^^^^\w(r' -r:r, (137) 



ijdrj dr'\w(r-ri)f^-^-^^^\w(r' -rj)\\ 



where is the angle between the dipole moment and vector r-r,-. Here, Cdd can take both positive 
and negative values by using a fast rotating field (see Sec. 14.1.3b . By changing the depth of an 
optical lattice, the transverse confinement, and the orientation of the dipoles, f, Uo, and {U'jj can 
be tuned independently. 

Let us begin by reviewing the case of U'J^ = 0, i.e., the on-site Bose-Hubbard model. At 
f = and commensurate filling, i.e., the average number of particles per site is an integer «, the 
interaction energy is minimized by populating every lattice site with exactly n atoms; therefore, 
the Mott insulator (MI) phase is realized. The energy cost required to create a particle-hole 
excitation in the MI phase is Uq (on-site interaction energy), and therefore, the MI state is gapped 
and incompressible. On the other hand, for nonzero t, the kinetic energy favors particle hopping. 
When t is sufficiently large to overcome the interaction-energy cost (~ Uq), the system undergoes 
a phase transition to the superfluid (SF) phase that is gapless and compressible. 

In the presence of the long-range interaction, the ground state of the extended Bose-Hubbard 



Hamiltonian (I136l l in a 2D optical lattice is investigated in Refs. ||65[ 16611 using a variational 
approach based on the Gutzwiller ansatz I*?) = H, 2)^0 /;!!«)/, where denotes the state with 
n particles at lattice site /. The coefficients /,| can be found by minimizing the expectation value 
(^|//ebh - /U E; «/l*l') under the constraint of a fixed chemical potential fi. In this calculation, a 
cutoff length is introduced and the interaction between lattice sites with distance larger than 
Ac is ignored. 

When dipole moments are polarized perpendicular to the 2D lattice, the dipole interaction is 
repulsive (attractive) for Cdd > (cdd < 0). If site / is occupied, it is energetically favorable that 
its neighboring sites are equally populated for negative Cdd and less populated for positive Cdd. 
Therefore, for negative Cdd, only MI and SF phases appear as in the case of the on-site Bose- 
Hubbard model, and local collapse occurs for large |cddl due to the strong attractive interaction. 
On the other hand, for positive Cdd, the DW phase emerges as an insulating state. According to 
the Landau theory of the second-order phase transition, phases with distinct symmetry breaking 
patterns cannot be continuously connected with each other Therefore, as the hopping increases, 
the DW state first changes to the SS state and then to the SF phase. 

When dipole moments are polarized along the 2D lattice, say in the y-axis, the DDI is 
anisotropic. The inter-site interaction for Cdd > (cdd < 0) is attractive (repulsive) in the y 
direction and repulsive (attractive) in the x direction. In this case, the particle density along the 
repulsive direction becomes periodically modulated, whereas it remains constant along the at- 
tractive direction, resulting in striped DW and striped SS states. Local collapse occurs for large 
|cddl for both positive and negative Cdd. 

In Ref. i66t] . the phase diagram in a 3D lattice is also investigated, where dipole moments are 
polarized along the z-axis. For Cdd > 0, the DDI is attractive in the z direction and repulsive in the 
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xy direction. Hence, the DW phase forms a checkerboard-type pattern on the xy plane and it is 
uniform in the z direction. The density pattern in the SS phase reflects the same symmetry as the 
DW phase. On the other hand, for Cdd < 0, layers of high-density sites align along the z direction 
due to the repulsive interaction. Hence, layered DW and SS phases emerge perpendicular to the 
z-axis. 

Recently, Menotti et al. IHsIl studied a polarized dipolar gas in a 2D optical lattice beyond the 
ground state. They have shown that there exist many metastable states in the insulator region de- 
pending on the cutoff length A^. The number of metastable states and pattern variations increases 
rapidly with the number of lattice sites. Metastable states also appear in magnetic domains in 
solid-state ferromagnets, classical ferrofulids, and quantum ferrofluids (see the next subsection), 
and they appear to be a characteristic consequence of long-range interactions. 

4.2.6. Ferrofluid 

Let us consider a repulsively interacting two-component BEC, in which atoms in component 
1 have a magnetic dipole moment, whereas those in component 2 are nonmagnetic. The m - -3 
and states of ^^Cr can be a candidate for such a system. We apply a magnetic-field gradient in 
the z direction, so that the two components phase-separate as shown in Fig. [8] 




Figure 8: Schematic illustration of a two-component BEC system in which atoms in component 1 have a magnetic dipole 
moment and those in component 2 are nonmagnetic. The magnetic dipole is polarized in the 7, direction by an external 
magnetic field. The two components are separated due to the field gradient B'{z) < 0. Reprinted from Ref. | 69]. 

The equilibrium configuration of this system can be found by solving the two-component GP 
equation in imaginary time li69ll . We found a stable hexagonal pattern of the density distribution 
of component 1 as shown in Figs.|9](a) and (b). We also found several metastable states start- 
ing from different initial states, such as stripe, concentric, and deformed hexagonal patterns, as 
shown in Figs. |9](c)-(e). The energies for these patterns are almost degenerate; however, each 
pattern is robust against small perturbations. 

A striking difference between the present system and magnetic liquids is that the present sys- 
tem is a superfluid that supports a persistent current. We found a metastable state with a finite 
circulation in component 1 and no circulation in component 2 in a stationary trap. Remarkably, 
the Rosensweig pattern emerges even in the presence of a superflow, indicating that the sys- 
tem simultaneously exhibits a diagonal and off-diagonal long-range order, although the latter is 
imbedded in the GP theory (see Ref. [(59J for more detail). 

4.3. Spinor-dipolar BEC 

Let us now consider a dipolar BEC with spin degrees of freedom. By "spinor-dipolar," we 
imply that the the direction of the spin of the order parameter is not polarized by an external field 
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Figure 9: (a) Isodensity surface of a hexagonal pattern formed in component 1 at B' = -1 G/cm. (b) Column densities 
for component 1 (left), component 2 (center), and the total system (right), (c)-(d) Metastable configurations found at 
B' = -0.3 G/cm, where the isodensity surface (top) and column density (bottom) for component 1 are shown. Reprinted 
from Ref ffill. 



but can vary in space. In such a situation, the system develops a nontrivial spin texture due to 
the DDI and self-organized magnetic patters may form. Possible quantum phases that arise due 
to the DDI have been discussed by using the single mode approximation |70] in which all spin 
components are assumed to share the same spatial mode. The spin ordering and spin waves in 
Mott states in ID and 2D optical lattices have also been discussed in Refs. 1711172117311 . 
In the second-quantized form, the magnetic DDI for a spinor BEC is given by 



Vdd = ^ jdrj dr' : F^(r)Q^y{r - r')F,{r') :, 



(138) 



where /"^(r) is defined by Eq. ( fTTI ) and Cdd - fioigFl-iB)^ /(^^), with gf being the Lande g-factor 
for the atom. The dipole interaction yields a non-local term in the GP equations. In fact, the GP 
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equations for a spin-1 BEC [Eq. (l52T l1 are given by 



in—— 

dt 



- + f/(r) — + 

1 1 



+ conip,,, + ci ■ f„„tff„ + Cdd ^ ft ■ f„miA«' (139) 



and those for spin-2 BEC [Eqs. (l78Tl-(l80ll1 are given by 



+ U{r)- pm + qm 



c 

V5" 



+ Conif/„, + Ci ^ f ■ f,„„(/r„ + -^Al/r*„, + Cdd ^ ft ■ fmnlAn (140) 



n--2 ^ n-—2 

Here, ft is the effective dipole field defined by 

b,= ^ dr'Y,Q,v{r-r')FAr'), (141) 

which is nonlocal, and Bdd = Ci^b/igffiB) works as an effective magnetic field. 

4.3.1. Einstein-de Haas effect 

We first consider the dynamics of a spinor-dipolar BEC at zero external magnetic field. Sup- 
pose that we prepare a spin polarized BEC in an external field. We consider what happens if we 
suddenly turn off the external field. In the absence of the DDI, a spin polarized BEC is stable 
because the total spin of the system will then be conserved. However, in the presence of the DDI, 
the situation drastically changes because the DDI does not conserve the spin angular momentum; 
it only conserves the total (i.e., spin plus orbital) angular momentum of the system. This can be 
understood from the fact that the DDI (I113l l is not invariant under a rotation in spin space due to 
the spin-orbit coupling term (d ■ r). On the other hand, because Vidir) is invariant under simul- 
taneous rotations in spin and coordinate spaces, the DDI conserves the total angular momentum 
of the system. Therefore, the DDI causes spin relaxation to occur by transferring the angular 
momentum from the spin to the orbital part; as a consequence, the BEC starts rotating. This is 
the Einstein-de Haas effect in atomic BECs i74ll . 

We numerically simulated such dynamics for the case of a ^^Cr BEC that is confined in an 
axisymmetric trap. Figure[TO](a) shows the time evolution of the spin angular momentum M- - 
J drY,]n=^-im\il/m{r)\^ and the orbital angular momenfiim = J drXl^^^if/l,(r)(^-ij^)il/^(r), 
where ip is the azimuthal angle around the symmetry axis. The initial state is assumed to be spin- 
polarized in the magnetic sublevel m = -3. We find that as time progresses, the spin angular 
momentum decreases with a concomitant increase in the orbital angular momentum, so that the 
total angular momentum of the system is conserved. 

In the dynamics, the order parameter is given by 

.Am(r,^,z) = e'(^=-'">'^77,„(r,z), (142) 

where (r,(p,z) are the cylindrical coordinates, 77,,, is a complex function of r and z, and is 
the projected total angular momentum on the z axis (7^ - -3 in the present case). The order 
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parameter (I142t is the eigenstate of the projected total angular momentum that satisfies 

(f, + = - j 4r„{r) = J,iff„{r). (143) 

The order parameters for magnetic sublevels m = -3, -2, and -1 at wf = 2 are shown in Fig.fTOl 
(b)-(d). Each spin component except the m - -3 state is a vortex state with the winding number 
7- - m, and the BEC as a whole forms a coreless vortex. 




Figure 10: (a) Time evolution of spin and orbital angular momentum due to the Einstein-de Haas effect, (b)-(d) Isopycnic 
surfaces of (b) i^_3, (c) 1//-2, and (d) at ojt = 2, where the color on the surface represents the phase of the order 
parameter. Reprinted from Ref. |74]. 

The physical origin of the spin relaxation is the Larmor precision of atomic spins around the 
dipole field given by Eq. ( I141I I. Figure[TT](a) shows the dipole field induced by a spin-polarized 
BEC in a spherical trap, and Figs.[TT](b) and (c) show the spin textures at cjt = 2. The whirling 
patterns in the upper and lower hemispheres exhibit opposite directions, reflecting the fact that 
the xy components of the dipole field are antisymmetric with respect to the z = plane. 

In the case of a ^^Cr BEC, the DDI compensates for the kinetic energy to form the spin 
texture. As discussed in the next section, a spin-1 '^^Rb BEC is ferromagnetic, and therefore 
a spin-polarized state is rather stable. In such a case, an external magnetic field applied in the 
opposite direction to the polarization causes a resonance between the kinetic energy and the 
linear Zeeman energy, leading to the Einstein-de Haas effect |75]. In the case of a spin-1 ^^Na 
BEC, the spin-exchange interaction favors a non-magnetic state, and therefore, it tends to destroy 
the spin-polarized state. 

4.3.2. Ground-state spin textures at zero magnetic field 

The DDI is known to create magnetic domains in solid-state ferromagnets. Similar structures 
are expected to appear in ferromagnetic BECs. Here, we consider ground-state spin structures. 
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Figure 11: (a) Dipole field induced by a spin-polarized BEC in a spherical trap, where the color of each arrow denotes 
the field strength (see the left gauge) and the solid curve indicates the periphery of the condensate. (b),(c) Spin textures 
on the (b) z = 2 fim plane and (c) z = -2 fim plane at ait = 2. The length of the arrows represents the magnitude of the 
spin vector projected on the xv' plane, and the color indicates |f |aj^^/A' with ai,o being the harmonic oscillator length and 
N, the total number of atoms. Reprinted from Ref. | 74]. 



The mean-field energy for the DDI can be rewritten as 

= ^ Z Z ^'^^^^2;.vW^v(-^) 
k fiv 

= '-^^Y.[3\k-F(kf-\F(k)\'], (144) 
^ k 

where F(k) is the Fourier transformation of F(r) and Qfiv(k) is given by Eq. (1125b . When the 
BEC is fully magnetized, i.e., \F{r)\ = fn(r), with / being the atomic spin angular momentum 
and «(r), the number density, the last term in Eq. ( 1144b makes a constant contribution determined 
by the density distribution. Then, the dipole-favored structure satisfies k ■ F(k) - for all k, or 
in real space, V ■ F(r) - 0. This condition is the same as that for solid-state ferromagnets fv^, 
and the corersponding stable domain structure is known as a flux-closure structure. In solid-state 
ferromagnets, the magnitude of magnetization is constant and the domain structure is constrained 
by boundary conditions such that the magnetization is parallel to the surface. In contrast, a 
spatial variation of \F(r)\ in proportion to the number density yields a spin texture in trapped 
condensates. 

As compared to solid-state ferromagnets, the unique feature of ferromagnetic BECs is that 
they exhibit the spin-gauge symmetry that generate a supercurrent by developing spin textures 
(see. Sec. 16.2. lb . The supercurrent, defined by Eq. (1264b . is generated by a spatial variation of the 
direction of magnetization. Therefore, even in the ground state in a stationary trap, a supercurrent 
can flow due to a dipole-induced spin texture. 
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Figure. [12] (a) shows the numerically determined phase diagram for a spin-1 ferromagnetic 
BEC |77] as a function of Rjp/^sp and Rtf/^m, where Rjp is the Thomas-Fermi radius; ^sp - 
h/ V2M|ci|«, the spin healing length; and ^dd - ^/ V2Mcdd«, the dipole healing length. We find 
three phases — flower (FL), chiral spin-vortex (CSV), and polar-core vortex (PCV) — whose spin 
configurations are shown in Fig. [12] (b), (c), and (d), respectively. The FL and CSV phases are 
eigenstates of the projected total angular momentum - 1, whose order parameter is given by 
Eq. (11421) . The orbital angular momentum for the order parameter (I1421 i is calculated as 

= [dr «^,;(r)(-/|-)-Am(r)= [ dr J] (Jz - m)\i/,mir)\\ (145) 

^ m=0,±l \ fl J m=0,±l 

Therefore, when 7; = 1, the BEC can have a nonzero orbital angular momentum depending on 
the population in the m = and -1 states. When the system size is smaller than ^dd, the kinetic 
energy, which hinders the spin texture from developing, is dominant, and the spin is almost 
polarized. Most of the atoms therefore reside in the magnetic sublevel m = 1 and the orbital 
angular momentum is quite small. As the BEC becomes larger, the spin texture develops to 
form a flux-closure structure. Consequently, the number of atoms in the m - and -1 magnetic 
sublevels increases, leading to an increase in the orbital angular momentum. The difference 
between FL and CSV phases is the symmetry of spin textures: the spin texture has chirality in 
the latter phase. Furthermore, an increase in the size of the BEC results in the phase transition to 
the polar-core vortex state with 7. = 0, because the Jr = I state costs a large kinetic energy when 
the population in the m - -1 state increases. In the J- = phase, a polar-core vortex appears at 
the symmetry axis. 

All three phases in Fig. [T2]can be realized in spin-1 ^^Rb BECs, which traces on the dotted 
line shown in Fig. [12] (a). The orbital angular momentum along this line is plotted in Fig. [13] as 
a function of the total number of atoms and trap frequency oj. The orbital angular momentum 
in the ground state has a significant value, and it increases up to approximately OAH in the CSV 
phase. 

For higher-spin BECs, similar spin textures are expected to appear in the ferromagnetic 
phase, although the mass current depends on the value of the atomic spin and the core structure 
in the polar-core vortex region is nontrivial. Spin textures in the limit of strong ferromagnetic 
interaction are discussed in Ref. [78J. On the other hand, even in the case of a non-ferromagnetic 
phase such as a spin-1 ^^Na BEC, a spin texture may appear if the dipole interaction is effectively 



enhanced in a pancake trap with a large aspect ratio 07911 . 



4.3.3. Dipole-dipole interaction under a magnetic field 

Now, we consider the effect of an external magnetic field on the DDI. Suppose that a ho- 
mogeneous external magnetic field B is applied in the z direction. Because the linear Zeeman 
term gF~B rotates the atomic spin around the z-axis at the Larmor frequency - gFf^BB/h, 
it is convenient to describe the system in the rotating frame of reference in spin space. The 
spin vector operators in the rotating frame are related to those in the laboratory frame by 

* = F^™" + iF^'°'^ = and F^"" = P.. Then, the dipolar interaction (TM can be 
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Figure 12: (a) Phase diagram of a ferromagnetic BEC in a spherical trap. A BEC of spin-1 ^^Rb atoms traces on the 
dotted line, (b)-(d) Cross sections of spin textures in the xy (upper) and xz (bottom) planes in the (b) flower phase, (c) 
chiral spin- vortex phase, and (d) polar-core vortex phase, where the color and length of each an'ow indicate \F\a^^/N and 
the spin projection onto the planes. Reprinted from Ref. liT^ . 



rewritten in terms of F*™*^ as 



V6 



(146) 



wheree = {r-r')/\r-r'\. When the Zeeman energy is much larger than the DDI, i.e., S-Wl » Cdan, 
the spin dynamics due to the dipolar interaction is much slower than the Larmor precession. In 
a large magnetic field, we may therefore use an effective DDI that is time-averaged over the 
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Figure 13: Orbital angular momentum for a spin-1 Rb BEC in a spherical trap with ti'ap frequency d). Reprinted from 
Ref. Wk. 



Larmor precession period Inlo)]^: 

7y CM r, , \r-r'f-3(: 

X : {f (™"(r) ■ F<""(r') - 3F?°''(r)F?°"(r')} 



'|5 



(147) 



More generally, an effective DDI under an external field applied in the B direction is given by 

^dd = - ^ ^ J^r Jdr' : F^'°'>{r)Q,Ar - r')f *""(r') :, (148) 

where the time-averaged integration kernel 

1 - 3(r ■ Bf 



-{6uv - 3BuBy) 



(149) 



corresponds to Eq. (1120b . In contrast to Eq. ( 1138b . the time-averaged interaction (1147b separately 
conserves the projected total spin angular momentum and the projected relative orbital angular 
momentum. However, the long-range and anisotropic nature of the dipolar interaction is still 
maintained in Eq. (11471 ). While the total longitudinal magnetization is conserved in the presence 
of a magnetic field, the anisotropic dipole interaction can induce spin textures in the transverse 
magnetization |0] . 

Equation (11491 ) provides information about how the dipole moments tend to align in the 
ground state. Longitudinal magnetization tends to become parallel in the direction of the mag- 
netic field, whereas it tends to be anti -parallel in the direction perpendicular to the magnetic field 
[Figs.[T4](a) and (b)]. On the other hand, transverse magnetization tends to be anti-parallel along 
the magnetic field and parallel in the perpendicular direction [Figs. [I4](c) and (d)]. Compar- 
ing the dipole energy for longitudinal and transverse magnetizations, the energy for longitudinal 
magnetization is twice as large (in negative sign) as that for transverse magnetization. 

The quadratic Zeeman effect also contributes to the ground state. As discussed in Sec. |3] 
the quadratic Zeeman effect in a spin-1 **^Rb BEC induces transverse magnetization. A stable 
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Figure 14: Alignments of (a)(b) longitudinal and (c)(d) transverse magnetization that are favored by the time-averaged 
dipole interaction. Due to the Larmor precession, transverse magnetization tends to be anti-parallel along the magnetic 
field and parallel perpendicular to the magnetic field. 



spin texture is determined by the balance between the quadratic Zeeman effect and the dipolar 
effect il mill. 

For simplicity, we consider a quasi-2D trap where the cloud size in the strongly confined di- 
rection is smaller than the spin healing length so that spin texture is uniform in this direction. We 
also assume that the magnetic field is applied parallel to the 2D plane and the total magnetization 
along the magnetic field is set to zero |25]. 

When the quadratic Zeeman energy is stronger than the DDI, the magnetization is restricted 
perpendicular to the external field. Then, the transverse magnetization forms a helix structure, as 
in the case of Fig. [T4](c). On the other hand, when the DDI is stronger than the quadratic Zeeman 
energy, the longitudinal magnetization forms a domain structure to minimize the DDI [Fig. [14] 
(b)]. However, these ground states differ from the magnetic patterns observed by the Berkeley 
group 1 25I 84 1. The equilibrium pattern for this system has yet to be understood. 

A further topic is the rotating spinor dipolar BEC in the presence of an external field. When 
the rotating frequency is close to the Larmor precession, the third or fourth line in Eq. ( I146I I 
remains in the time-averaged dipolar interaction instead of the second line in Eq. ( 1146b . For ex- 
ample, when a spin-polarized BEC is rotated, the spin relaxation is expected to occur resonantly 
with the Larmor frequency. 



5. Bogoliubov Theory 

Quantum and thermal fluctuations as well as external perturbations induce excitations from 
the mean-field ground state. When the excitations are weak, they can be described by the Bo- 
goliubov theory. We express the field operator as a sum of its mean-field value ij/m and the 
deviation from it 6\j/m. 

lAm = lAm + S^m (m = /, / - 1 , ■ ■ ■ , -/). (150) 

Here, the mean-field part can be calculated from the Gross-Pitaevskii theory described in Sec.[3j 
The basic concept of the Bogoliubov theory is to substitute Eq. (1150b in the second-quantized 
Hamiltonian and retain the terms up to the second order in Skfi,,,. The resulting Hamiltonian can 
be diagonalized by means of the Bogoliubov transformations. The Bogoliubov transformations 
are canonical transformations that couple particle-like excitations with hole-like ones via a Bose- 
Einstein condensate and they are non-perturbative. The resultant dispersion relation therefore 
depends on the coupHng constant in a non-analytic manner. 
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We consider a spatially uniform system (U(r) - 0). In this case, it is convenient to analyze 
the system in Fourier space. We obtain the Fourier expansion of the field operator: 

4',n^-^Yak,ne""', (151) 

where Q is the volume of the system. Then, the noninteracting part of the Hamiltonian is written 
as 

Hq = ^(e* - pm + qm^)hkm, (152) 

km 

where - lT?k^/2M and «fo„ = a^j^Jikm- Because the interaction Hamiltonian comprises the 
particle density, spin density, and spin-singlet pair operators, it is convenient to introduce the 
corresponding Fourier-transformed operators: 



Pk 



\ drh{r)e ^^al^Uq+km (153) 

^ qm 

f drF{r)e-"'' = ^ f„„fl;„,fl,+fc, (154) 

qmn 

V27TT J d/rAoo(r)e-'*- = ^(-l)-^-"'fl„„fl,-t,-;„, (155) 



where it is understood that m runs over ;« = /,/- 1, ■■■ , -/. 

We assume that the BEC occurs in the ^ = state, and decompose the terms in Eqs. il52i - 
( 11551 ) into the ^ = component and the rest. Then, Eq. ( 1152b is written as 

//() = '^(-pm + qm^)ho„, + ^ (ft - pm + qm^)hk,n- (156) 

As shown below, the interaction Hamiltonians involves quadratic forms of the Fourier-transformed 
terms in Eqs. ( I153l l-( ll55l l. They can be decomposed similarly as 



Y,-p\Pk- = N{N-\) + Y^:p[pk:, (157) 

k ki^O 

}]--FlFk: = :FIf,:+Yj--PIP'c--^ (158) 

k ki=0 

J]-^l^k: = ■.AlAo-.+Y.-AlAk-. . (159) 



k ktO 



In Eq. ( 1157b . the dominant terms in pk^^o - Yjqm o\mOq+km are the terms with q - and 
q = -k: 



Pk^O 



- Xi*-^*'"^*'" "-fc«"Om) ^ Dk + d[^, (160) 



where - a'^Jiikm describes the density fluctuations from the condensate. Substituting 
Eq. ( 1160b into Eq. (1157b . we obtain 

^ : plpk :- N(N-l) + Y^: {iblbk + DkD k + b\D[^) : . (161) 
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In Eq. (I158I I. the term Fq can be further decomposed into 

Fo^F+ ^ imna\„flq+kn, (162) 

where F - T^mn fmn^l^^On IS the zeroth-order term and the last term is of the second order. On 
the other hand, the term F^^o can be approximated as 

+ alj^^^aon)- (163) 

mn 

This operator describes the spin fluctuations from the condensate. Substituting Eqs. < ll62l t and 
( 11631 ) into Eq. ( 1158b and retaining the terms up to the second order, we obtain 

ft: qtd^inn ktO,mnm' n' 

Similarly, Eq. ( 11591 ) can be approximated as 



^ : A^At : = 45+5_ + 2 5+ ^ (-l)^""afo„fl_;t,_„ + H.c. 

k V i#0,m 

+4 ^ {-\r^%^a^,„al_Jk.-n, (165) 

where H.c. denotes the Hermitian conjugate of the preceding term and 

S- = (5+)"'' = i 2(-l)-^"'"flomao,-„,. (166) 

m 

In the Bogoliubov approximation, we replace oom by V^^m, which satisfies the following 
normahzation condition: 

By subtracting the last term, we can construct a number-conserving Bogoliubov theory without 
introducing the chemical potential. In keeping with the Bogoliubov theory, we retain the last 
terms in operator form. 

5.7. Spin-1 BECs 

The interaction Hamiltonian for the spin-1 BEC is given by Eq. (|23] ) or ( l24l ). In Fourier space, 
it is expressed as 

k 
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where cq = (^o + 2g2)/3 and ci = (g2 - goM'i. Using Eqs. (I156I I, ( 1161b . and (1159b . we obtain the 
effective Hamihonian for the spin-1 BEC: 

+ ^ (ei - pm + ^m^) hkm 

Co + Ci 



2Q 



^ : {2dId, + DkD k + d[d[^) : 



5+ ^ (-1/ "Clkma-k-m + H.c. 



(169) 



kmn 



The first three terms on the right-hand side were studied in Refs. 11851 U, l86|]; the third and 
fourth hnes describe the density and spin fluctuations, respectively, and the last line describes 
the coupling between them. In the following discussions, we shall discuss the case for ^ = 
unless otherwise stated. The case for ^ is discussed in Ref. ifioll . 



5.7.7. Ferromagnetic phase 

When c\ < 0, the ground- state phase is ferromagnetic, and we have 



(170) 



Substituting this into Eq. ( 1169b . we obtain 

1 



= -{co + ci)n{N-l)-pN 
1 

ekfiki + -j{cQ + ci)n{ak\a-k\ + fikx + H.c.) 



+{ek + p)hkQ + (e* + 2p- 2cin)hk-i] ., 



(171) 



where n = N/Q.. The Hamihonian (1171b shows that the m = and m - -1 modes describe 
the massive transverse spin and quadrupole spin excitations, respectively, with the dispersion 
relations given by 



"k,-l 



= ek+P, 



£k + 2p — 2cin. 



(172) 
(173) 



To diagonaUze the Hamihonian of the m - 1 mode, we perform the Bogoliubov transforma- 



tion: 



ak = bk cosh 6 - b , sinh 9. 



(174) 
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Substituting this into Eq. ( 1171b and requiring that the coefficients of the off-diagonal terms such 
as bkb-k vanish, we obtain 



sinh(20) = 



gn 



cosh(26') = 



(175) 



where 



-02 



(176) 



and eic = + (cq + ci)n- The diagonalized Hamiltonian then becomes 



k*0 



where 



Ell ^ ^Jek[ek+2(co + Cl)n] 



(177) 



(178) 



describes the coupled density and longitudinal spin excitations. 

The sum over k in the first line of Eq. (1177b diverges for large k — an artifact that arises due 
to the use of a delta function that does not correctly describe the short-range behavior of the real 
potential. For large k, the summand becomes 



^Jekiek +2gn) - ek - gn 
Subtracting this term, the sum converges, giving 



{gnf 
2ek 



(179) 



Iv^/ I ignf\ 128 ria' 

-^.Je,(e, + 2gn)-e,-gn+—j^—gn^—N. 



(180) 



Thus, the ground-state energy of the system is 



gnN 



1 + 



128 



■pN. 



(181) 



It follows from this that the sound velocity is given by 



c = 



MN 




1 +1 



(182) 



where the leading term is the Bogoliubov phonon velocity and the last term is called the Lee- 
Yang correction to it. 
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The Bogoliubov approximation takes into account virtual pair excitations from the mean- 
field, and therefore, the properties of the ground state are very different from those of the mean 
field. A mean-field ground state in the ferromagnetic phase is 



|vac). 



(183) 



On the other hand, the Bogoliubov ground state is given by 



lO) ~ exp 



01 



|vac). 



(184) 



where cp^ - A^[l - (8/3) yJnaf/Ti], Vk - 1 -H q - ^fciJc^+2) with q = k^l{8nain), and the 
summation is performed over the entire range of ky and k-, but the summation over k^ is restricted 
to the positive side to avoid double counting. In contrast to the single-particle state in Eq. (1183b . 
the Bogoliubov ground state is a pair-correlated state. In momentum space, the Bogoliubov 
ground state exhibits pair correlation having opposite momenta; in real space, it describes the 
pairwise repulsive interaction. To show this, let us write down the many-body wave function in 
the coordinate representation: 

"Viriai,--- ,rA,Q'A,) = (vac|(A„,(ri)---0'„,(ri)|<D). (185) 

Substituting Eq. (1184b into this and ignoring the terms of the order of 1 IN, we obtain iIstIi 



,rNl) 



exp^-|(37r - 8)A^ ^/«fl^) 



■ exp 



(186) 



where r,y = |r, - r j\ and ^ = {^nainY'^l^ is the correlation length. Thus, the many-body wave 
function decays exponentially when any two particles approach closer than the scattering length 
or the correlation (or healing) length. 

5.7.2. Antiferromagnetic phase 

When c\ > 0, the ground state is antiferromagnetic, and we have 



^0 = 0, Ifiil 



2 2 2N 



fikj 



(187) 



ki^Om 



where y = p/{cin) and the last term is retained in order to construct a number-conserving Bo- 



goliubov theory. Substituting this into Eq. ( 1169b . we obtain 

H^^ = ^oHN - 1) + ici« [y^N - l) - pyN 
+ ^ (e* + cin6,„fl)hkm + ^ ['^oiDkD-k + ^^Dt) 



kn) 



+ciiSkS-k + SlSk + n|^i^_i|a*:,oa-i,o) + H.c.} , 



(188) 
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where 



Dk 



z 

m 

z 



y 

1 + 2^ti 



1 - 2'^k.-\ 



(189) 
(190) 



We note that when y i^O, the density and spin excitations are coupled because [Dj, 5[] - yll. 

The m - mode is decoupled and its dispersion relation can be calculated in a manner similar 
to that for the m - 1 mode of the ferromagnetic phase, with the result 



(191) 



The m - ±1 modes are coupled, and their dispersion relations can be obtained by solving the 
Heisenberg equations of motion for the m - +\ modes [6]: 



Cl 



(192) 



In the absence of an external magnetic field {p - y - 0), Eq. (11921) reduces to the excitation 
spectra of the uncoupled density and spin waves: 



(193) 



5.7.5. Domain formation 

The MIT group observed the formation of metastable spin domains llssll . They first prepared 
all atoms in the m - 1 state and transferred half of them to the m = state by irradiating the 
rf field. Then, letting the system evolve freely in time while applying a uniform magnetic field 
to prevent the m = -1 component from appearing due to the quadratic Zeeman effect, they 
found that spin domains develop with the two components being aligned alternatively. It was 
discussed L88i.l89i.i9Q] that this phenomenon is due to the imaginary frequencies of the excitation 
modes. Here, we use the Bogoliubov theory to confirm this hypothesis and to derive a general 
expression for the dispersion relation. The experimental conditions in effect amount to setting 



1^1 1' = l^ol' = 5 - ^ J](nki + «w), i-i = 0. 



(194) 



Substituting this in Eq. (I169l l, we obtain 

Co + Cl 



H 



n{N -\)- —nN 



Zl Con 3ci«\^ / Con cin\^ 



ki=0 

Co , , , 

+ yn(flW)fl-JK) + h-c.) 

Co + Cl 



n{ak\a-k\ + StiS-M + cij^^a-ko + h-c) 
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(195) 



The eigenspectrum can be obtained by writing down the equations of motion for a±ic i and a±k,Q- 
Rewriting the resulting equations in terms of A - 0^,1 + a^j^ j and B - cikfi + a^_i^ y, we obtain 



Assuming that A oc e r^*' and B oc e r^*', we obtain 



(196) 



+2v (2u + v)ek - V M + -V 



(197) 



where u = CQn/2 and v = cin/2. For parameters of the MIT experiment, n ~ 10 cm 



flt ~ 29 A, and ~ 26A 119111 . we find that v/m 1. In Eq. ( I197l i. ignoring higher-order powers 
of v/u, we obtain 

(Ekf ^4 + (2m + v)ek + 2(u + v)ek, 

where the plus (minus) sign corresponds to the density (spin) wave. We see that the energy of 
the spin wave becomes pure imaginary for < v, implying the formation of spin domains. The 
corresponding wavelength defines the characteristic length scale of the spin domains as 



A, = 



37T 



^ (a2 - ao)n ' 



(198) 



This result agrees with that of Ref. 190t], except for a numerical factor Using the above pa- 
rameters, we obtain A ~ 18_um; this is in reasonable agreement with the observed value of 
approximately about 40 fim Km . 

5.2. Spin-2BECs 

We again assume that the total number of particles is conserved: 



Then, the one-body part of the Hamiltonian becomes 

Ha = -pNf, + Z " P^"" - ^^^J 



i^O in 



where 



(199) 



(200) 



(201) 



50 



The interaction Hamiltonian of a spin-2 BEC is given by Eq. ( [30b and its Fourier transform 
becomes 



V 



hi X ■ • ^""^ ■ '^^^^ • +y^*^*) ■ 



(202) 



Using Eqs. ( 11561 ). ( 11611 ). ( I164l l. and ( 11651 ). we obtain the effective Hamihonian for the spin-2 
BEC; 



where 



nk,, 



ki:0 mn 



C2« 



S-* ^ ^(-l)'"flft,ma-A-m + H.C. 



V k^O tn 



ki^Q mn 



mn 



(203) 

(204) 
(205) 



5.2.7. Ferromagnetic phase 

We assume that only the m - 2 magnetic sublevel is macroscopically occupied. Then, we 
have F, - 2no2 = 2{N — 2i#0m «tm)^ ™d the one-body part becomes 



i7o = -2pN + ^ [et + (2 - m)p] 



(206) 



ki^O,m 



Because the spinor that is normalized to unity is = 6m^2, only the z-component is nonzero in / 
and i_ = and Eq. ( 12031) becomes 



V 



= ^84nN + n ^ [coni,2 + 2ci(2hk,2 - 2%,o - 2>hk-\ - 4% ,-.2) 



+ ^C2hk-2 + ^■(,ak,2a-k,2 + «I,2^-t,2) 
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(207) 



where g4 = cq + 4ci. Combining Eqs. (12061 ) and ilOli . the total Hamiltonian is given by 



(ek + g4n)hk,2 + -g4n(ak,2a-k,2 + 2^-4,2) 



+ ^ [(e* + p)nk,i + {ek + 2p - 4cin)%o + (e* + 3;? - 6ci«)nt _i 



+ (et + 4p - 8cin + 2c2n/5) nt,_2] ■ (208) 
It follows from this Hamiltonian that the m - 2 mode gives the Bogoliubov spectrum 

El.2,k = ^el+2g4nek. (209) 

For the system to be stable, the Bogoliubov excitation energy ( 12091 ) must be positive. This 
implies that the s-wave scattering length for the total spin-4 channel must be positive; 

g4 = ^fl4 > 0. (210) 

This condition is the same as the one required for the ferromagnetic mean field to be stable, that 
is, the first term on the rhs of Eq. ( 12081 ) being positive. We also note that the Bogoliubov spectrum 
(12091 ) is independent of the applied magnetic field and remains massless in its presence. This 
Goldstone mode is a consequence of the global U(l) gauge invariance due to the conservation of 
the total number of bosons. 

On the other hand, all the other modes give single-particle spectra: 

ELuk = ^x + P' (211) 

ELo,k = ek + 2p-4cin, (212) 

^L-i,* = ek + 3p-6cin, (213) 

El=-2,k = ffc + 4p - 8ci« + 2c2n/5. (214) 

For the system to be stable, the single-particle excitation energies must also be positive. This 
implies that p > 2c\n and p > (2ci - C2/10)n; these conditions are the same as the ones for the 
ferromagnetic phase to be the lowest-energy mean field. 

5.2.2. Antiferromagnetic phase 

We consider the situation in which only the m - ±2 modes are macroscopically occupied. 
By requiring that the total magnetization be the normalized spinor is given by 



^2 = e'^^ V ^ + f ' ^2 = e'*-- y ^ - f , ^0 = ^.i = 0. (215) 

Because of the quantum depletion of the condensate, F. has to be replaced by {N-Yjki=Om ^km)fz, 
and the one-body part of the Hamiltonian is therefore given by 



Hq = -pF, + ^ (ejt - pm)hkm 

= -pNfz + J][ek- pirn - /,)] hkm. (216) 



ki^O,m 
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The interaction Hamiltonian can be evaluated by noting that p ^ f-{c\ - c2/20), and the total 
Hamiltonian is given by 



H^^ = ^(co + C2/5)nN-^pNf, 



+ 



^ {(^k + gAn\^2?-)hk,2 + (£k + g4n\(-2f')hk-2 



+ I — (^2 ak,2a-k,2 + ^_2«i,-2fl-it,-2) 



+(co - 4ci + 2c2/5)n(f2^*2«*,2«-it -2 + ^2C2^k,2^k,-2) + H.C.]} 



e;i + (ci - C2/5)n + -(ci + C2/10)n/, 



nt,i- 



ejt + (ci - C2/5)n - -(ci + C2/10)n/2 



+2(ci - C2l5)n{^*2C2^k,\a-k-\ + H.c.)} 
+ ^ - C2«/5)njt,o + '^{ClC-2^k,oa-kfi + H.c.)| . 



(217) 



We thus find that the spectrum is classified into three cases: the coupled m - ±2 modes, coupled 
m - +\ modes, and single m - mode. 

We first investigate the eigenspectra of the m - ±2 modes. The relevant part of the Hamilto- 
nian is 



h2 ^ ^ {skafika + ek-2hk-2 



+ — (^2 ^k,2a-k,2 + C2«i,-2«-it,-2 + H.C.) 



+ (co - 4ci + -C2 I ni^2C2^k,2a-k,-2 + (2C-2H,2^k,~2 + H.C.) ^ , 



(218) 



where ek,±2 - ^k + ^4n|^±2p- Here, the complex numbers (±2 can be made real by means of 
unitary transformations fl±jt,±2 — > a±k,±2e"^-- ■ The resulting equations can be diagonalized by 
writing down the Heisenberg equations of motion for ak,±2 + a ^ ^j, and seeking solutions of the 

form ak.±2 + a' 1^ oc txpi-iE^^^^tlh). The resultant eigenspectra read 



/-pAF <2 _ ^ 



eu + g4n + g4n - 



gA 



4 



(219) 



For the ground state to be stable, these energies must be positive; this implies the conditions 
c\ - C2/2O > and co + C2/5 > 0. The former condition is met if the antiferromagnetic phase is 
the lowest-energy state, whereas the latter condition ensures that the antiferromagnetic phase is 
mechanically stable, that is, the first term on the rhs of Eq. ( I217l i is positive. The spectra ( 12191 ) are 
massless regardless of the presence of the magnetic field. This is because they are the Goldstone 
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modes associated with the U(l) gauge symmetry and the relative gauge symmetry, which is the 
rotational symmetry about the direction of the applied magnetic field. 

Next, we investigate the eigenspectra for the coupled m = + 1 modes. The relevant part of the 
Hamiltonian is 

= ^ ki.ini.i + £k,-in-k,-i + (5flt,ifl-*:,-i + H.c.)) , (220) 
where 6 = 2(ci - C2/5)n^2^*2 and 

et,±i s£t+(ci-j^n + ^ (ci + Y^jn/z- (221) 

The eigenspectrum can be obtained by writing down the Heisenberg equations of motion for 3^1 
and S-i-i, and seeking solutions of the form exp{+iEk,±it/h). The eigen matrix is 2 x 2 and it 
can be solved straightforwardly, giving 




+ ^e2+2(c,-|)«e, + i(c,-|)'«2/:2. (222) 

The spectra become massive in the presence of a magnetic field. For the ground state to be stable, 
the spectra must be positive definite. Hence, we obtain the conditions ci - C2/2O > p/2n > 
and C2 < 0; these are in agreement with the conditions for the antiferromagnetic phase to be the 
lowest-lying one. 

The spectrum for the m = can be obtained by means of the Bogoliubov transformation: 




(223) 



This mode becomes massive in the presence of a magnetic field. 

At zero magnetic field, all five dispersion relations in Eqs. (I219l l. (1222b . and ( I223l l become 
gapless and linear, reflecting the fact that in the absence of the magnetic field, the ground state is 
degenerate with respect to five continuous variables. Note, however, that only the first four are 
Goldstone modes because the continuous symmetry concerning the last mode does not reflect 
the symmetry of the Hamiltonian but the hidden SO(5) symmetry 



5.2.3. Cyclic phase 

There are several mean-field solutions for the cyclic phase. Here, we consider the Bogoliubov 
excitations from the following mean field; 




(224) 
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where 4>2 + ip-2 - '2<po = ±7r. Because - and - 0, the interaction Hamiltonian becomes 



V 



^ + ^pNfz - cmf^ Yj Yj 

Co 

2V 



+ S Z ■ ('^'^'^-'^ + ^kDk + H.c.) : +c,nf, ^ J] 



+ ^ YYfufmn{a\hQ,„a[^^ak,j + a^.a'^Jikja-k,n + H.c.) 



0,m «;,-m 



(225) 



^^0 nin 



The interaction Hamihonian shows that the spectra are divided into two parts: the m - +1 
coupled modes and the m - ±2 and m - Q coupled modes. The eigenspectra for the former 
modes can be found by writing down the Heisenberg equations of motion for fl±A:,±i with the 
result 



(226) 



These excitation energies are always positive semidefinite and massive in the presence of a mag- 
netic field. 

The eigenspectra for the m = ±2 and m = coupled modes can be obtained by writing down 
the equations of motion forX± = 0^,2 ± ^1^2' ^± ^ - ^-to' ^ ^k-2 ± ^Lk-2- 

desired eigenspectra are given by 



"±2,0 



2c2n 



and 



4 + 

(^±2,0)^ = Q et+ Con + — z^cjn ± -v|(con)2 - (4 - 3/?) coCi«2 + 



(4 + /?)2 



(ci«)2 



(227) 



(228) 



The solutions in Eq. ( I228I I are always positive semidefinite and massless even in the presence 
of an external magnetic field. This is a consequence of the fact that the mean-field solution is 



degenerate with respect to at least two continuous variables. The solution in Eq. ( 1227b is massive 
and independent of the external magnetic field. 

5.3. Dipolar BEC 

In the presence of a magnetic dipole-dipole interaction (DDI), the Hamiltonian is 
^"/'^'' }]rjr)^-^V^ + U(r)-pm + qm'^yUr) + jn(rf + cmr)f 

+ ^-f Jdrdr' F,{r)Q,Ar-r')F,{r'\ (229) 
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where Cdd = yUo(gF/^B)^/(47r) and 



(230) 



In general, the long-range nature of the DDI yields rich physics in a trapped system, such as 

nontrivial equilibrium shape and spin textures. The Bogoliubov spectrum in the presence of the 

DDI has been studied using a Gaussian ansatz and numerical analysis iilsllil 111. Here, we 

discuss the Bogoliubov spectrum for a spinor dipolar EEC under an external field 

We assume that the system is strongly confined in the z-direction and free in other directions 

so that it is effectively two-dimensional. Consequently, the order parameter is decomposed as 

follows: 



ifr{x,y,z) = h{z)i//{x,y), 
where h(z) is assumed to be normalized to unity: 

h\z)dz = 1. 

oo 

Substituting Eq. (123 li t into Eq. (I229l t, we obtain an effective 2D Hamiltonian: 



H: 



2D 



where p = (x, y). 



- - r 

J — CK 



h\z)dz. 



(231) 



(232) 



(233) 



(234) 



and /fdd is the two-dimensional version of the time-averaged dipolar interaction ( 11491 



2D 



Cdd 



2jhHz)dz 

Here, the integration kernel 0j^^\k2D) is given by 



fdpf dp'J]F,{p)FAp')J]e-''^°-^''-P'^Q^^\k2D). (235) 



'C3d) 



(236) 



Qf^\k2u) = / § / dzhHz) J dz'hHz')e''''''-''^Q,Ak: 

= J dzh\z) J t/z'/z2(z')e*^'"*y[l-3(B-^3D)'](V-3Mv)> (237) 

where Q/ivik^^^^) is the Fourier transform of Eq. ( 11491 ) and B is a unit vector along the external 
field. 

To proceed with the calculations, we assume that 



1 

(2;rt/2)i 



(238) 
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After some algebraic manipulations, we obtain 



Q(k) =1-3 cos^ 9b + 3 yfn&rfc{kd)kde- 



,k-S 



X l^cos^ 6b - sin^ Osik.x cos <pB + ky sin ., 



(239) 



(240) 



where 0b and (pB specify the direction of the magnetic field: B - (cos (pB sin 6b, sin ^b sin 0^, cos 0^). 
With this kernel, the Gross-Pitaevskii equation in the presence of the magnetic DDI can be ex- 
pressed as follows: 



ot 



V + U{r) + am + cnn 

2M 



(241) 



where the linear Zeeman term is eliminated by moving into the rotating frame of reference in 
spin space with the Larmor frequency, and b is an effective dipole field given by 

W = \ dp' ^ e'*=« <^-'''>aA:2D) ^(V - 3B^By)Fy(p'). (242) 



Note that because the spin and orbital parts are decoupled in the dipole kernel (12371 ). it is not 
necessary to choose the coordinate axes in real space that coincide with those in the spin space. In 
what follows, we take the spin quantization axis along the magnetic field, which is not necessarily 
parallel to the direction of the strongly confining potential, and their relative direction is given 
by B = (cos ^B sin 6b, sin (pB sin 6b, cos Ob)- 

In the following discussions, we discuss the ferromagnetic phase (ci < 0) such as the f - I 
hyperfine manifold of ^^Rb. As in this case, we assume that the coefficient of the quadratic 
Zeeman effect is positive: ^ > 0. When the total longitudinal magnetization is set to zero, the 
magnetization is transverse due to the quadratic Zeeman effect. In fact, it can be shown that 



f<1 










V(l+§)/2 









(243) 



(244) 



gives a stationary solution of Eq. ( 1241b with U (r) - 0, where 

- ^ q 

^ 2«(|ci|+CddQ(0))' 

The expectation value of the spin vector for this spinor is nonvanishing only for the x-component: 

F, = Yj^^>pTntMn^ = ^«(l - §2)<5^,,. (245) 

inn 

To investigate the Bogoliubov excitation, we expand the order parameter from the mean-field 



stationary state ( 12431) : 



Zl 



km 



-i(k-r-ci)ktfh) 



(246) 



57 



Substituting Eq. (1246b into Eq. (I2411 l. we obtain the 6x6 eigenvalue matrix equation: 

-)(:). 

where and vt are three-component vectors and //q i is the 3x3 matrix. The spin wave mode 
M, V oc (1, 0, -1)^ decouples from the other two modes and the eigenfrequency is obtained as 

col ^[ek + q- Cdd«(l - q)l2Q(k) + QiO)]) [e^ + Cdd«(l + q)[Q(k) - Q{0)]) . (248) 

Note here that due to the anisotropy of Q(k), the right-hand side of Eq. (12481 ) becomes neg- 
ative depending on the direction of k. This result implies that the initial state given by ( I243I I 
is dynamically unstable. Figures [T5] (a)-(c) show the distributions of |Imwt| (red) and -RecJk 
(blue) in momentum space which correspond to the dynamical instability and Landau instability, 
respectively. For > 0, the sign of a>k is determined so that the corresponding eigenmode 
satisfies the normalization condition m^m^ - v[vt — 1. 







(a) 


(b) 


(c) 




d 


80mG. no helix 




160mG. no helix 




2nOmG. no lielbt 


"a 


d 










• 
• 




La 
o 

=? 


IlULJfr] Ol^^Ml.4s"' 




|Imwj:| 0^™0.44s'' 




ni — — ru ^-i 



-0.05 0.0 0.05 
fe/27r [;jm"^] 



Figure 15: Real and imaginary parts of Ci)i^ in Eq. 12481 for the case of f/ = 1 /im and n = 2.3 X lO'* cm"'. The region with 
nonzero imaginary part shows the dynamical instability, while the region with Retj;t < shows the Landau instability. 
The external magnetic field is (a) 80 mG, (b) 1 60 mG, and (c) 200 mG. 

The anisotropic distribution of the instability suggests the existence of the periodic pattern in 
real space. Recently, the Berkeley group observed the crystallization of transverse magnetization 
in a quasi two-dimensional ^^Rb BEC [25]. In one of the experiment, they started from the 
initial condition in which the magnetization form a helix and let the system evolve with time. 
The Bogoliubov spectrum (12481) obtained above does exhibit dynamical instabilities that exhibit 
various patterns such as a helix and stripes. However, it fails to explain the Berkeley experiments 
in terms of the time and length scale. 

6. Vortices in Spinor BECs 

A scalar BEC can host only one type of vortex, that is, a U(l) vortex or a gauge vortex. 
However, a spinor BEC can host a much richer variety of vortices. We begin by discussing why 
this is the case. The properties of a vortex can be characterized by looking at how the order 
parameter changes along a loop that encircles the vortex. Because the order parameter must be 
single-valued, it should satisfy 
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(249) 



where is a 2/ + 1 -component order parameter; ip, the azimuthal angle which specifies the 
location of the loop; and ^, a set of "coordinates" which, together with (f, completely specify 
the order parameter. For the case of a scalar BEC whose order parameter is a single complex 
function i/r - |i^|e"*, the vortex solution to Eq. ( |249t is 

ilf{ip,^)^W"-^, (250) 

where n„, - 0, ±1, +2, ■ ■ ■ is the winding number If «,v + 0, the density at the vortex core must 
vanish. The corresponding superfluid velocity is given by 

^ [^*V<A-(V^-)^], (251) 



2Mni 

where M is the atomic mass and n - Substituting Eq. (12501) in Eq. (12511) . we find that 

h 

Vs = n,,—'V(f>. (252) 
M 

Therefore, the circulation of the superfluid velocity is quantized in units of k - h/M. 

Due to the internal degrees of freedom, spinor BECs allow several solutions other than 
Eq. ( 1250b . depending on the symmetry of the order-parameter manifold. A general order pa- 
rameter of a spinor BEC is described as 

VH^= VHe'^f/(ff,Ar)^o, (253) 

where n is the particle density; (p, the global gauge; U(a,p,y) - e^'^" e^'^y^ e^'^'"^ , an SO(3) 
rotation with spin matrices f^: ,. ^ and Euler angles a, /3, and y; and ^o, a representative normalized 
spinor for the order parameter. Hence, the spacial variation of a,/3, and y as well as characterize 
the structure of a vortex. In the following, we investigate what types of vortices are allowed in 
each phase. A general classification of topological excitations based on homotopy theory is given 
in Sec.O 



6.1. Mass and spin supercurrents in spinor BECs 

The mass current density jt^^**) for a spinor BEC is defined as the sum of currents in all spin 
components: 

= 5^ Z ['^^(^'A,,,) - (V^:>„,] . (254) 

m 

It can be shown from the GPEs that satisfies the continuity equation: 

^ +V-/™^^^> = 0. (255) 
at 

In a similar manner, we define the spin current density jj,^'""' as 

fT' = ^ YS^,),nn [C(V<^„) - (VC)<^«] ■ (256) 

mn 
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At zero magnetic field, satisfies the continuity equation for the component of the spin: 



dt 



(257) 



not conserved, and Eq. ( 1257b is replaced with the equation of motion of spins given by 



In the presence of an external magnetic field, the spin angular momentum perpendicular to it is 



(258) 



2qn 



where B is a unit vector in the direction of the magnetic field and N is the 3x3 symmetric 
nematic tensor defined in Eq. (fTTT i. The dipole interaction induces an effective magnetic field and 
therefore modifies the first term in Eq. ( 12581 ) 



The superfluid velocity field for the particle and spins, defined as v^™"*^^ = /'"^^'*'/« and 



(spin) _ .(spin) 



/n, respectively, are described in terms of a normalized spinor ^ as 



„(spin) 



(259) 
(260) 



In contrast to the case of a scalar BEC, the circulations of Eqs. ( |259t and ( 12601 ) are not always 
quantized, as shown below. 

6.2. Spin- 1 BEC 

The rotation matrix for the order parameter of a spin-1 BEC is given by 
/'e-'X<*+r) cos2 ^ _£± sin/3 e-"'"-^^ sin^ 



U(a,/3,j). 



V2 

COSyS 



i(a+y) 



^sin/3 e'l'^+y' cos-' ^ 



(261) 



The spin-1 BEC in the absence of an external magnetic field has two phases: ferromagnetic and 
polar Below, we discuss vortices in each phase. 

6.2.1. Ferromagnetic phase 

A representative order parameter of a spin-1 ferromagnetic BEC is = (1,0, 0)^, where 
T denotes the matrix transpose. Substituting this and Eq. (1261b in Eq. ( 12531 ). we obtain a general 
order parameter of a spin-1 ferromagnetic BEC: 



feno 



,'(0-r) 



^sin^ 

2 P 



e sm 



(262) 



2 ^ 



The direction of the magnetization for this state is given by 



s = — = (sin j6 cos a, siuyS sin a, cos /3) . 
fn 
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(263) 



The linear combination - y in Eq. ( 1262b reflects the spin-gauge symmetry which impHes that 
the rotation in spin space through angle y is equivalent to the gauge transformation by 4 



For simplicity of notation, we set cp' - (/> -y. We substitute Eq. ( I262l i in Eq. ( 12591 ), obtaining 



s 



^(mass,F) ^ ^ [ y^' _ p y^,] (264) 

Thus, the spatial variation of the spin gives rise to a supercurrent, and vice versa. 

In the case of a scalar BEC, the superfluid velocity is given by Vj - (h/M)'V(p; it is therefore 
irrotational, i.e., V x Vs = 0. However, this is not the case with the ferromagnetic BEC. In fact, 
using 

V X (fa) = (V/) X fl + / (V X a), (265) 

we find that 

V X v^™'''^^ = sinyS ( VyS X Va) . (266) 

Introducing unit vectors to = (cos /3 cos a, cos yS sin a, - sin /3) and hq = (- sin a, cos a, 0) which 
satisfy s - toX Hq, we obtain Vy6 = 2^ to^^s^ and sin j6Va - 2^ u^^Vs^. Then, we can express 
Eq. ( 12661 ) in terms of s as 



= — 2^ ^/j/i'<5vv'(fO/iVi;,0 X (mqvVv) 

h -^-^ 

Afivfi' v' 

= ^ [ei/ivi^(V.?^ X Viv) - (fo^Vi^) X (MovViv)] 

/l/JV 

fivA 

This relation is known as the Mermin-Ho relation |95]. V x v*"^'*^'^-' is nonvanishing because of 
the contribution from the Berry phase. To see this, we rewrite Eq. ( 12641 ) as 

^(mass,F) - A(l _ coSyS) Vof = ^ V(0' - ff) 

and integrate both sides along a closed contour C. Due to the single-valuedness of the order 
parameter, the integral of V(0' - a) should be an integral multiple of 2n. Hence, we obtain 

^(mass,F) . _ ^ _ cos/3)Vq' ■ dt = —H,, ill,, : integer). (268) 
TV/ TV/ 

The term ^( 1 - cos y6) Va ■ dt gives the Berry phase enclosed by the contour C. We thus conclude 
that in the ferromagnetic BEC, the circulation alone is not quantized, but the difference between 
the circulation and the Berry phase is quantized. 
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As a special case, let us consider the order parameter ( 12621 ) with a - -n„ip and (p' - n^.tp; 
then, Eq. (|262] | changes from ^ = (e'^""'^, 0, 0)^ to (0,0, 1)^ as yS changes from to n. This 
implies that the In^ vortex is topologically unstable. On the other hand, if we take a = -n„(p 
and (f>' = («„, + 1)^, then Eq. (|262] | changes from ^ = (e'(2n„+iv^ q, 0)^ to (0, 0, e''ff as /? changes 
from to n. Thus, the 2n„ + 1 vortex is unstable against the decay to the singly quantized vortex 
which is stable 101 • 

We can utilize these properties to create a doubly quantized vortex from a vortex-free state 
by changing j0 from to tt as a function of time while keeping a - (f and (p' - -ip [96, 97]. Such 
a scheme is realized if a spin polarized BEC is prepared in the m - 1 state under a quadrupole 
field B(r,(p,z) = (Bj.(r, z)cos(-(/;), Z?x(r, z)sin(-(^), B;(r, z))^, where (r,(p,z) are the cylindrical 
coordinates and > 0. At t - 0, B; is set to be positive and much greater than B±. As B, is 
changed adiabatically from » to B- <g: -B^, the atomic spins follow the direction of the 
local magnetic field asa - (f and /3 = aictan(B±/B^), resulting in the transformation of the order 
parameter from (1,0,0)^ to (0, 0, e"^"^)^; thus, a doubly quantized vortex is imprinted. Using 
this method, the MIT group have observed multiply quantized vortices ll98ll . 

Any axisymmetric configuration of if/ is homotopic to one of the following two types of 
vortices: a coreless vortex [Fig.[l6](c)] and a polar-core vortex [Fig.[l6](b)]. When (f>' - +a, the 
vortex is coreless; in particular, the order parameter for the case of (p' - a = (f is given by 



coreless 



COS 



2 13 



V2 



2 

sinyS 
■sin^^ 



(269) 



in which the singularity can be removed by choosing yS = at r = 0. Let us consider the con- 
figuration in which /3(r - Q,(p,z) = and P(r = ro,(p,z) = tt, where ro is the radius at the 
(cylindrical) boundary of the system. Then, the spinor (1,0, 0)^ at the origin with no vortex sin- 
gularity and (0, 0, e^"^)^ at the boundary with a doubly quantized vortex attached to the inverted 
spin configuration. The corresponding superfluid velocity is nonsingular at the origin; in fact. 



(mass.F) 



= —(1 -C0Sf3)ifi, 

Mr 



(270) 



where ip is the unit vector in the azimuthal direction. The coreless vortex has been observed by 
the MIT group [99J, where a vortex was imprinted using a quadrupolar field. 

On the other hand, when (p' - and a - ip, the vortex core has a singularity in the ferromag- 
netic order, where the order parameter at r — > oo is given by 



polar-core 



c2 



(271) 



2 



Because the m - +1 components have nonzero vorticities, the density in Eq. ( 12711 ) must vanish 
at the center of the vortex. In fact, the free-energy analysis shows that the vortex core is filled 
by a polar state in which only the m - component is occupied; hence, this vortex is called a 
polar-core vortex. As a consequence, only the spin density vanishes at the vortex core and the 
number density remains finite. A variational order parameter that describes both the vortex and 
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its core simultaneously is given by 



1 -/2(r)cos4| -g2(r)sin'* 
e'^g(r) sin^ f 



(272) 



where f(r) and g(r) are the variational functions subject to /(O) = g(0) - and /(oo) - g(oo) - 1 
and it can be obtained so as to minimize the Gross-Pitaevskii energy functional. The spontaneous 
formation of the polar-core vortices has been observed by the Berkeley group lilOOll . 




Figure 16: (a) Order parameter for the spin-1 feiTomagnetic phase, where the aiTow indicates the direction of the local 
magnetization, (b) Polar-core vortex given by Eq. 127 It . where the ferromagnetic order parameter has a singularity at the 
vortex core which is filled by a polar state, (c) Coreless vortex given by Eq. j269t with /3(r = 0) = and fi(r = kq) = n/l. 



To calculate the spin superfluid velocity ( 1260b . it is convenient to use the Cartesian coordinate 
in spin space ^ = {^x,ty,^zV rather than the spherical coordinate ^ = (^i, ^o, ^-i)^- They are 
related to each other by 



1 



-10 r 
-/ -/ 

V2 0, 



(273) 
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In the Cartesian coordinate, spin matrices are given by (ffi)vA - -i^fivA, where ij,v,A - {x,y,z}. 
The ferromagnetic order parameter ( 12621) is then given by 

^feno ^ _f!l(iy + ia^^ ^ __L^i + iu), (274) 

V2 V2 



where t - cos (f>'tQ - sin (f>'uo, u = sin (p'to + cos (f>'uQ, and to and uq are defined below Eq. ( 12661 ). 
Substituting Eq. ( 12741 ) in Eqs. ( 1259b and ( 12601 ) and using some vector calculus formulae, we 
obtain 

^(.ass,F)^|^^^^^^^ (275) 



M 

(spin,F) _ ^ ^(massB^ ^ 
2M 



v,A 



The first term on the right-hand side of Eq. ( 12761 ) is the flow of spins carried by the particle, and 
the second term results from the spin flow induced by the gradient of spins that is perpendicular 
to s. The spin superfluid velocity v^'^p"'^* satisfies the foUowing relations: 

V ■ vJ,'P"'^^ = s^V ■ v^™^^'^^ + v^^^^'P^ ■ Vs^ - ^(s X V^s)^, (277) 

2 ■ vl^^P"'^^ = V ■ v^'"^^^'^), (278) 
f 



Z^4'^x»')<^'"''1=0, (280) 

where we have used the Mermin-Ho relation ( 1267b to derive Eq. ( |279b . 
6.2.2. Polar phase 

We next consider the order parameter of a spin- 1 polar BEC, where = (0,1, 0)^. A general 
order parameter of the polar state is given by 



coSjS 
4=sinfl 



(281) 



The fact that i^p"'^"^ does not depend on y reflects the SO(2) symmetry around the quantization 
axis d - (sinyScosa, sinjesina, cosyS); exp(-/f ■ d)tl/P°^'^'^ - i/rP"'"'". Substituting Eq. (12811) into 
Eq. ( 12591 ). we obtain 

^(mass.P) ^ ^y^^ (282) 



Thus, unlike the case of the ferromagnetic phase, the circulation is quantized; however, the unit 
of quantization is half the usual h/M. To understand this, let us consider a loop that encircles a 
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vortex with a fixed radius. Then, each point on the loop is specified by the azimuthal angle if. 
We note that the single-valuedness of the order parameter (1281b is met if we take, for example, 
a - ^ - n^yipjl and fi - nl2. Then, the order parameter at r — > oo is given by 



lalf- vortex 



-1 ^ 





(283) 



where n„ is an integer It follows then that the circulation is quantized in units of h/(2M) rather 
than the usual h/M: 



,(mass,P) 



de : 



h 

2m' 



(284) 



The underlying physics for the half quantum number is the Z2 symmetry of the order parameter 
of the polar phase; Eq. ( 12811 ) is invariant under n gauge transformation {(p ^ <l> + ^) combined 
with a 7T spin rotation around an axis perpendicular to d; d — > -d. More generally, the single- 
valuedness of the order parameter is satisfied if (p changes by n„7T along the closed loop, where 
n„ is even if d{(p - In) = d((p = 0), and odd if d((p - 2n) - -d{tf = 0). Thus, the polar 
phase of a spin-1 BEC can host a half-quantum vortex for the case of - 1 jlOlll . which is 
also referred to as an Alice vortex or Ahce string 110211 . The dynamic creation of half-quantum 
vortices was discussed in Ref. i l03[l . As in the case of the polar-core vortex in the ferromagnetic 
phase, the atomic density can be finite at the vortex core of the half-quantum vortex by filling the 
ferromagnetic state [(1,0, 0)^ for the case of Eq. ( 12831 )1 at the core. 




Figure 17: (a) Order parameter for the polar state, which is invariant under n rotation about the Ct axis together with 
the 7T gauge transformation, (b) Half-quantum vortex of the polar phase of a spin-1 BEC. The order parameter rotates 
through an angle rr about the Ct axis defined in (a) as one circumnavigates the vortex. To satisfy the single-valuedness 
of the order parameter, this round trip must be accompanied by a phase change in ir (phase dift'erence between A and B), 
resulting in a half-quantized vortex. 

To calculate the spin current, we again use the Cartesian coordinate in which the polar order 
parameter is given by - e"^d. Substituting ^po'""" in Eq. ( 1260b . we obtain 



(spin.P) 



M 



Z 



(285) 
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If d is restricted in a two-dimensional plane perpendicular to a unit vector a (in spin space), the 
circulation of 2^ a^vjj^'""'^* is quantized: 



,^^(spin,P) ^ ^ ~^"»' ("i^' • integer). (286) 



Here, the unit of the quantization is h/2M due to the spin-gauge Z2 symmetry. The rotation and 
divergence of the v^^"""'^* are calculated to be 

V.v« = -A(rfxV%. (288) 



Note that at zero magnetic field, the continuity equation (I257l l becomes 

-V.(«v«) 

M 



dt 



^ Z Z f/jy^ [(5,n)(?v(5,(?i) H- ndyidididj)] 

i-x,y,z. vA 



M ■ 

i-x,y,z vA 

dF III 

57 = M Z ^ 5,(«5,rf). (289) 

!=.v,y,z 

This result implies that 2, diindid) has to be either zero or parallel to in a stationary d texture. 
Otherwise, the local magnetization would develop with time and the order parameter would no 
longer belong to a manifold for the polar state. The rf-field around the half-quantum vortex of 
Eq. ( 12831 ) is d - (cos«,vi^/2, sinn,^.(^/2,0), which satisfies Yjidiindid) = -n\,dl{Ap-) if n(r,(p) is 
axisymmetric, i.e., ■ Vn = 0, where is the unit vector in the azimuthal direction. Therefore, 
the half-quantum vortex can be a stationary state. 



6.2.3. Stability of vortex states 

The stability of a vortex state can be studied by using the multicomponent GP equations. 
An important feature of the spinor BEC is that each of the components can accommodate dif- 
ferent vortices (see, for example, Eq. ( 128 11 1): however, this gives rise to relative phases between 
spinor components, and from the energetics point of view, certain conditions must be met. From 
Eq. ( l44l l. we find that only the term ci \F\^/2 depends explicitly on the relative phases. Expressing 
the term in terms of components, we obtain 

yl/l' = y[(I<AiI'-I'A-iI¥+2W'(I«AiI' + I!A-i|') 

+4|.Ao'Ai'A-il COS Ofi +X-1 - 2;ro)] , (290) 

where Xm = ^^gi^m) is the phase of the m-th component. To minimize the energy, +X-\ - '^Xo 
must be nn, where n must be even if ci < and odd if ci > 0. Substituting ;t:-„, - x'm +Xmf^ where 
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if is the azimuthal angle, we must have;tfj + x'-\ ~ "^x'q - ^« and;^^" + x"-\ - '^x'o - 0- Ass umin g 



that the maximum vorticity is 1, we find that the following three vortex states are allowed 1110411 : 



(xlXo,X-i)-(hl,l), (1,0,-1), (1,1/2,0), (291) 

where the last state is allowed only if i/tq - 0, and therefore, it is described as (1, none, 0). 

The (1,1,1) state is the usual singly quantized vortex state with an empty vortex core, where 
the order parameter is given by 

e'^diAile'^i, |i/ro|e'^», |(A-i|e*'i)^- (292) 

The (1,0,-1) state has a vortex in the m = 1 component and an anti-vortex in the m - -1 
component with the vortex core filled by the m -Q state, where the order parameter of this state 
is given by 

me''^'^''\mM-i\e-'"'''^''V. (293) 

This is the polar-core vortex in the ferromagnetic phase or an integer spin vortex in the polar 
phase. The (1, none, 0) state corresponds to the Alice vortex [Eq. ( I283l ll whose order parameter 
is given by 

(|iAi|e'^'+'^\0, |(A-i|e'^'')^, (294) 

where the vortex core is filled by the m - -1 component. Except for the (1,1,1) state, the vortex 
core is filled by a non-vortex component which plays the role of a pinning potential, stabilizing 
the vortex state in the absence of an external potential [.1041 . 

6.3. Spin-2BEC 

The rotation matrix is 

f/(«,y8,r) = 

g-2i(ff+y)(-4 -2e-'<2"+y>C^5 y/6e-^'"C^S^ -2e-'^^"-^^CS^ g-2i(a-y)g4 

2g-/(<^+2y)(-3^ _2g-/(a+y)c252 _ ^g-'" gin 2/3 26-'''"-^^ S ^ -le-'^^-^y^CS 

V6e-2'Vc252 ^e-'>sin2yS i(l+3cos2/3) - ^e'^ sm2/3 ^e^'yC^S 

2e'("-^y)CS^ 2e'''"-y^C^S^ ^e'" sin 2/3 -2e'^"^y^C^S^ -2e'^"^^T'^C^S 
giKa-y)gA 2e'^^-y^CS^ ^/6e^'"C^S^ 2e'(2«+y)C^5 g2/(a+7)(-4 

(295) 

where C = cos | and S = sin|. The spin-2 BEC has four mean-field ground-state phases in 
the absence of the magnetic field: ferromagnetic, uniaxial nematic, biaxial nematic, and cyclic. 
The uniaxial and biaxial nematic phases are degenerate at the mean-field level ( fTST i jlOSll . The 
uniaxial and biaxial nematic phases are also referred to as the antiferromagnetic phase in Refs. 
lit] and the polar phase in Ref . | IJ | . 
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6.3.1. Ferromagnetic phase 

The representative state of the ferromagnetic phase is - (1, 0, 0, 0, 0)^. Substituting this 
and Eq. (1295b in Eq. ( 1253b . we obtain 



-ferro 



,i(fl>-2y) 



V6cos2f sin^f 
2e'" cos I sin^"* | 

2ia ■ 4 yS 



(296) 



Similar to the spin-1 case [see Eq. ( I262b l. the order parameter has the spin-gauge symmetry 
with y now replaced by 2y because the magnitude of the spin is 2 rather than 1. The superfluid 
velocity is similarly obtained as 



(mass.F) 



M 



(297) 



whose circulation is not quantized as in the case of the spin-1 ferromagnetic BEC. Applying the 
discussions about the stability of the vortex states for the ferromagnetic spin-1 BEC, it can be 
shown that the 4n,v + m,v vortex (m„, = 0, 1, 2, 3) is unstable against the decay to the m,v vortex. 

The relationships for the spin superfluid velocity corresponding to Eqs. ( I277b -( I280I ) change 
to the following ones: 



Vv 



{spin,F) 
(spin,F) 



= 2 



. ^(mass.F) ^ j,(mass,F) . y _ ^ y2j)^ 



2V- v' 



(mass,F) 



vAcr 



J]S,[VX v'r''''] = V X v(-^-P> = ^ J] e,.,s,iVsy X VS.) 



(298) 
(299) 

(300) 

(301) 



6.3.2. Uniaxial nematic phase 

The representative state of the spin-2 uniaxial nematic phase is ^™""' - (0,0, 1,0,0)^ and 
substituting this and Eq. ( 12951 ) in Eq. ( 1253b . we obtain the general order parameter of this phase: 



-2e"'" sinyScosyS 
^(3cos2/3-l) 
2e"^ sin/? cos y6 
e^'^sin^jS 



(302) 



Similar to the case of the spin-1 polar phase, the spin degrees of freedom do not contribute 
to the superfluid velocity that only depends on the spatial variation of the gauge angle: 



(mass,uniax) 



= — V0. 



M 



(303) 
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This is a general property for the non-magnetized phase in spinor BECs. 

The shape of the uniaxial order parameter in the spin space is shown in Fig. [18] (a). As one 
can see from Fig.[T8](a), the uniaxial nematic phase has a SO(2) symmetry around the direction 
of d - (cos a sin sin a sin cos /I). The uniaxial nematic phase also has Z2 symmetry as in the 
case of the spin-1 polar phase. However, in the present case, Z2 symmetry is independent of the 
gauge angle (p: is invariant under a tt rotation around an axis perpendicular to d. Hence, 

although the uniaxial nematic phase can accommodate a vortex, as shown in Fig.[l8](b), the mass 
circulation for this vortex is zero. We call it a 0-1/2 vortex, where and 1/2 indicate the gauge 
transformation and rotation angle in spin space around the vortex, respectively. An example of 
the order parameter around the 1/2-spin vortex is given as 



spin-vortex 



V6n 





J2J3 




(304) 



where (p is an azimuthal angle around the vortex. 

The spin superfluid velocity for the uniaxial nematic phase is calculated to be 



(spin,uniax) 



M 



(305) 



If d is restricted in a plane perpendicular to a, the circulation of the spin superfluid velocity 
2// fl/iV^^'""'""""'' is quantized in terms of a unit of 3h/M. 

Because the spin and gauge is completely decoupled in this phase, the mass circulation is 
quantized in terms of a unit of h/M. An example of such a vortex is (0, 0, e"^, 0, 0)^. 



(a) 



U(I) 



(b) 









Figure 18: (a) Order parameter for the spin-2 uniaxial nematic phase, which has SO(2) symmetry corresponding to 
rotations about d. The uniaxial nematic phase is also invariant under a n rotation about the C2 axis, i.e., this phase has a 
Z2 symmetry, (b) The order parameter configuration for the 0- 1/2 vortex, around which the order parameter rotates by rr 
about the C2 axis in (a). 



6.3.3. Biaxial nematic phase 

The representative state of the spin-2 biaxial nematic phase is ^Jj'" = (1/ V2, 0, 0, 0, 1/ 
and substituting this and Eq. (|295b in Eq. ( 1253b . we obtain the general order parameter of this 
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phase: 



biax 



- J sin^yS^ cos 2y - i cos p sin 2yj 
e"'" sin y8(cos y6 cos 2y - i sin 27) 

-^sin^y6cos2')' 
-e'" sin y6(cos y8 cos 2y + / sin 2y) 
e^'" [( 1 - i sin^ yS) cos 2y + / cos /? sin 2y] ^ 



The superfluid velocity only depends on the spatial variation of the gauge angle: 



(306) 



(mass,biax) 



M 



(307) 



The profile of ^i'™ is depicted in Fig.[T9](a). 

When we rotate the order parameter by 7r/2 about the C4 axis in Fig. [19] (a) as we circum- 
navigate a vortex, the order parameter changes its sign, as shown in Fig. [19] (b) (see, the phase 
difference between A and B). Therefore, this vortex has a mass circulation of h/(2M). We call 
this vortex the 1/2-1/4 vortex, where 1/2 and 1/4 denote the gauge transformation and the spin- 
rotation angle, respectively, around the vortex. If we rotate the order parameter by n about the 
C2 axis around a vortex, then it goes back to the initial state, as shown in Fig. [19] (c). In this 
case, the vortex has no mass circulation and it is referred to as a 0-1/2 vortex. On the other 
hand, if we rotate the order parameter by n about the axis, the sign of the order parameter is 
reversed [Fig. [19] (d)]. This vortex has a mass circulation of h/(2M), and it is called a 1/2-1/2 
vortex. Because the roations about C4, C2, and are non-commutable, these vortices obey the 
non-Abelian algebra, which is discussed in Sec. 16.61 

The vortices shown in Fig.[T9]are the simplest examples, and in general, the order parameter 
can vary in a more complex manner. However, in order to meet the single-valuedness, the mass 
circulation in the biaxial nematic phase has to be quantized in terms of a unit of h/(2M). The 
spin superfluid velocity is very complicated, and its circulation is not quantized. 



6.3.4. Cyclic phase 

A representative state of the spin-2 cyclic phase is fg^'^'"^ = (1/2, 0, // V2, 0, 1 /2). Figure l20l 
(a) shows the profile of the order parameter of the cyclic phase *!* = F™^^^^'"^. One can clearly 
see that it has a three-fold symmetry about the (1,1,1) axis [C3 axis in Fig.|20](a)]. This implies 
that the system has a one-third vortex [d06ll . The rotation matrix about the (1,1,1) axis through 
angle 6 is given by 

(308) 

For 6 = 27r/3, this reduces to e'*'''^^e"*^g^'^'"^; the single-valuedness of the order parameter is met 
by the gauge transformation of(p — > (p + 27T/3. Because this is one-third of the usual 27t, the cyclic 
phase can possess a one-third vortex. Similarly, for 9 = 47r/3, Eq. ( 13081 ) reduces to e^'''^-'e"*^g*''^'"^; 
the single-valuedness of the order parameter is met by the choice of A(p - 47r/3. Thus, the cyclic 
phase also has a two-third vortex. The cyclic phase can also accommodate the 1/2-spin vortex, 
around which the order parameter rotates by tt about the direction of one of the lobes [Fig. [20] 
(c)]. The vortices in the cyclic phase also obey the non-Abelian algebra (see Sec. 16.61 ). 
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Figure 19: (a) Order parameter for the spin-2 biaxial nematic phase, which has one symmetry axis of the fourth order 
(C4) and two axes of the second order (C2 and C^). (b)-(d) The configuration of the order parameter around (b) a 1/2-1/4 
vortex, (c) a 0-1/2 vortex, and (d) a 1/2-1/2 vortex. 



6.4. Rotating spinor B EC 

The stationary state of a rotating BEC can be obtained by minimizing F - E - f2L,, where 
the rotation axis is assumed to be in the z-direction. In the case of a spin-1 BEC, we have 



r dr Y,r,r, 



2v72 



" 2M 



The stationary state is obtained by requiring that 

6F_ _ 



(309) 



(310) 



where //,„ is the chemical potential for the m-th component that is determined so as to conserve 
the total number of particles - J dr^jm liAmP' total magnetization - J dr(\il/i\^ — 
and total projected angular momentum L, - J dr 2, t(/*{!^if/i, where - -ih dldw . The vortex 
lattice formation in the presence of an external rotation i s disc ussed in Ref . llOTIl . The vortex 
phases of spin-2 BECs under rotation is discussed in Ref. lIlOSll . 



6.5. Spin-polarized dipolar BEC 

The dipole interaction is long-rang e and anisotrop ic, and it has significant effects on the 
ground state of rotating BECs 1 109, 1 10. .1 1 ll[TT2lll 13ll . Let us assume that all the dipoles of the 
atoms are polarized in the same direction. Then, the interaction is described by 



y(r) = CQ6{r) + CiA- 
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1 - 3 cos^ ( 



(311) 




Figure 20: Order parameter of the cyclic phase (1,0, i'/ V2, 0, 1)^, where the blue, yellow, and green colors denote 
2rg[Zm ^'"'fm*^'^'"^] = 0. 2;r/3, and 4^/3, respectively. The order parameter has a two-fold symmetry about the x, y, and z 
axes and a three-fold symmetry about the (±1, ±1, ±1) axes, (b) The 1/3-1/3 vortex, where the order parameter is rotated 
by 2n'/3 about the C3 axis from A to B in a clockwise direction, results in the gauge transformation of -2n/3. (c) The 
0-1/2 vortex, around which the order parameter is rotated by n about the y-axis. 



where the first and second terms on the right-hand side are the i-wave contact interaction and the 
dipole interaction, respectively, and 6 is the angle between the dipole moment and the relative 
coordinate r between two atoms. Thus, depending on the direction of the polarization, the dipole 
interaction can be either attractive or repulsive. This anisotropy of the interaction leads to the 
deformation of the vortex core and the structure of the vortex lattice. 



Yi and Pu 111 1211 considered a situation in which a system in an axisymmetric harmonic trap 
is rotated about the z-direction, and the dipoles are polarized in the x-direction in the frame of 
reference corotating with the system, and found that the vortex core is deformed into an elliptic 
shape. This is because the dipoles undergo repulsive and attractive interaction in the x- and 
y-directions, respectively. 



Cooper et al. 1 109,] and Zhang and Zhai 111 1 111 considered a situation in which the interparticle 
interaction is so weak that the lowest Landau level approximation holds, and found that with an 
increase in the strength of the dipole interaction, the vortex lattice undergoes structural changes 
from triangular to square and to stripe configurations. Such cha nges were not found in Ref. [1 1211 . 
and possible reasons for the discrepancy are discussed in Ref. 11 1 1 311 . 

In the aforementioned work, the dipole interaction is treated as a perturbation to the s-wave 
contact interaction. Very recently, the Rice group [26.1 used a shallow zero crossing in the wing 
of a Feshbach resonance of ^Li in the \F - l,mf = 1) state to decrease the scattering length to 
as less as 0.0 Iab- This opens up several new possibilities on the dipole-interaction-dominated 
BEC. An interesting question is how the degeneracy of a fast rotating BEC is lifted by a weak 
but genuine dipole interaction and what is the corrsponding many-body ground state. 
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6.6. Non-Abelian vortices 

The topological charge of a vortex determines how the order parameter changes as one cir- 
cumnavigates the vortex. For example, if the phase of the order parameter of a U(l) vortex 
changes as e'"'^, the topological charge of this vortex is said to be n. One important feature of 
a spinor BEC is that the order parameter has a geometrical shape in spin space, as illustrated in 
Figs. [T6l420l Take, for example, the polar phase. This phase has two generators e^^ ''^ and e"'C2, 
where the former describes the rotation about the symmetry axis and the latter describes the n- 
rotation about the perpendicular axis followed by the gauge transformation by n [see Fig.fTTl'a)]. 
These two operations commute with each other; thus, vortices of the polar phase of a spin-1 BEC 
are Abelian. Similarly, the vortices in the ferromagnetic phase are also Abelian because there 
is only one symmetry axis set by the magnetization. Therefore, the spin-1 BEC can host only 
Abelian vortices. 

It is predicted ill 14ll that the cyclic phase of a spin-2 BEC can host non-Abelian vortices. 
To understand this, we note that the order paramet er of the cyclic phase is invariant under the 
following 12 elements of the tetrahedral group T il06ll : 1, Ix - e'^'", A = e'^'", h — e'^-" , 
C = e2m/3g-2;r/(f,+t;+f,)/3V3^ j j j q2 ^ j qI^ ^nd LC^ , where f^ is the ^l component 

of the spin-2 matrices (// = x,y,z) [see, Fig. l20ll . Furthermore, these elements are classified into 
four conjugacy classes: 

• (I) integer vortices: {!); 

• (II) 1/2 - spin vortices: {/j, /,,, /,); 

• (III) 1/3 vortices: {C,I,CJyCJ,C]; 

• (IV) 2/3 vortices: {C^, l,C^,IyC^, LC^}. 

Topological charges in the same conjugacy class transform into one another under the global 
gauge and/or spin rotation. The order parameter of a 1/2-spin vortex can be constructed from 
that of a non-vortex state ^^non-voitex 

^l/2voitex _ ^('f^(/3^non-voitex (312) 

where if is an azimuthal angle around a vortex line. 

The non-Abelian characteristics of the vortices manifest themselves most dramatically in 
the collision dynamics. In general, when two vortices collide, they reconnect themselves, pass 
through [Fig. |2Tld)], or form a rung that bridges the two vortices [Fig. l2TT b).(c)l. When two 
Abelian vortices collide, all these three cases are possible, and one of them occurs depending 
on the kinematic parameters and initial conditions. However, when two non-Abelian vortices 
collide, only a rung can be formed; however, reconnection and passing through are topologically 
forbidden because the corresponding generators do not commute with each other. In fact, the 
nonzero commutator of the two generators gives the generator of the rung vortex. Figure |22] 
illustrates a typical rung formation. When a core of 1/3 vortex in the cyclic phase is filled with 
the ferromagnetic state, which is possible for a certain parameter set, it is possible to observe 
such dynamics of vortex lines by using a phase-contrast imaging technique that can detect local 
magnetization. 
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Figure 21: Collision dynamics of two vortices, (a) Initial configuration, where A or B represents an operation that 
characterizes the corresponding vortex (a set of spin rotations and gauge transformations for the case of a cyclic BEC). 
The vortex on the right end, which is connected to B, is identified as ABA^K The configuration in (a) is topologically 
equivalent to (b) and (c), where a rung is formed. If the vortices are a pair of vortex and anti-vortex, i.e., A = B"', the 
rung in (b) disappears, giving rise to reconnection, whereas the rung in (c) corresponds to a doubly quantized vortex that 
costs a large kinetic energy. If A and B are commutative, passing through is also possible because the configurations 
of (a) and (d) will then be topologically equivalent. However, when A and B are not commutative, the collision always 
results in the formation of a rung. 




Figure 22: Numerical simulation for the collision dynamics of non-Abelian vortices. 



7. Topological Excitations 

When a system undergoes Bose-Einstein condensation, some symmetries of the original 
Hamiltonian are spontaneously broken. The classification of the symmetry breaking can be 
carried out systematically using a group-theoretic method. The symmetry of the order parame- 
ter determines the types of possible topological excitations. In this section, we present a brief 
overview of some basic notions and theoretical tools for investigating symmetries of the order 
parameter and topological excitations. 

7.1. Symmetry classification of broken symmetry states 
7.1.1. Order parameter manifold 

The order parameter for a spin-/ BEC is described with a set of 2/ + 1 complex amplitudes 
= {(Am), and the order parameter manifold M is denoted as M = C^^^' . In this section, we treat 
a normalized spinor ( - II 1^ II with its order parameter manifold given by 

„ = U(2/.l,= Hl>i^HVli>. ,313, 

Consider a group G of transformations that act on M while leaving the mean-field energy 
functional invariant. We assume that the total spin and number of particles of the system are 
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conserved. Then, the energy functional is invariant under SO(3) rotations in spin space and U(l) 
gauge transformations. The full symmetry group G for a spinor BEC is therefore given by 

G = SO(3)f X U(1V, (314) 

where subscripts F and (p denote the spin and gauge, respectively. Any element g & G can be 
represented as g = e'^e^'^-"e^'^'^e^'^''>' , where a, /3, and y are Euler angles [see Fig.|23](a)] and 
ffj. j; are (2/ + 1) x (2/ + 1) spin matrices. When an external magnetic field is applied, e.g. in the 
z-direction, G reduces to 



Gb =U(lk xU(lV 



(315) 



because the system will remain invariant only under rotations about the z-axis apart from the 
gauge transformation. Any element g' G Gb is represented as g' = e"^e^'^-''. 

For each ^ e M, the isotropy group (or the little group) Gf is defined as the subgroup of G 
that leaves ( invariant: G^ = {g e G \ g^ = For example, spin-1 polar phase ^p"''"" = (0, 1,0)^ 
is invariant under a spin rotation about the z-axis, g-'t>'^poiar _ ^poiar^ therefore, G^poiar - 
U(1)k. 

Next, we define the orbit G(^) of f by G(^) - {g^ \ g e G), i.e., G(Q is the set of all points 
obtained as a result of letting all transformations in G act on ^. Because elements of G do not 
change the mean-field energy, the orbit G(^) constitutes a degenerate space in M. For the case of 
a spin-1 polar BEC, a general element of the orbit G(^p°'''0 is given by 



■polar 



(0] 




'-^sinyS^ 


1 




coSjS 






l$sin^J 



(316) 



The orbit is thus parametrized by (p and d = (cos a sin sin a sin yS, cos /?) that specifies the 
quantization axis. The manifold of d (a three-dimensional unit vector) is a unit sphere which is 
referred to as S^, the two-sphere [see Fig.|23](b)], whereas the manifold of (p is U(l). Therefore, 
we conclude that the order parameter manifold for the spin-1 polar state is given by Tjpo''" s 
5^xU(l)^. 





Figure 23: (a) Euler rotation with Euler angles (ff,/?, y). (b) Two-sphere S" which is the manifold of a topological space 
of a three-dimensional unit vector d. 

Mathematically, such a manifold is defined as a left coset of G(. For a subgroup H c G, the 
left cosets of H are defined as G/H - {gH \g e G). The coset space is identical to the orbit of 
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H, and it forms a group if i/ is a normal subgroup of G, i.e., gHg - H for e G. The coset 
space G/H is equivalent to G/H' if and only if H and H' are conjugate, i.e., H' - gHg^^. On the 
other hand, one can prove that the isotropy groups of all points on the same orbit are conjugate, 
i.e., Ggg - gG^g^^ for 3g e G. Therefore, the coset space does not depend on the location in the 
orbit. We use H for a representation of the conjugate class of the isotropy groups gG^g^\ and 
define the order parameter manifold as 

R^G/H. (317) 



Here, H represents the remaining symmetry of the ordered state, and R describes the broken 
symmetry, or a manifold of the degenerate space. 

We next discuss the discrete symmetry of the order parameter manifold. The spin-1 polar 
phase has spin-gauge coupled discrete symmetry. In fact, ^p°'^' - (0,1, 0)^ is invariant under 
a :7r-rotation about an axis perpendicular to the z-axis followed by the gauge transformation of 
e'". This can also be understood from the fact that Eq. ( I316l l is invariant under the simultaneous 
transformations of {(p, rf) — > (0 + n, -d). A group for such an operation { 1 , e'^e"'''") is isomorphic 
to Z2 = {0, 1) because (e"^e"'*^^''^)^ - 1. Therefore, the correct isotropy group of is the 
dihedral group (Doo)f-4 - ^(^)f^ ^ (^2)^,0, and the order parameter manifold is given by 

SO(3)f xU(l), 



R' 



,polai' 



U(1U X (Z2)f,, 



^^xU(l)^ 

(^2)^,0 



(318) 



where \J(1)f xi ('L2)F,<j> implies that when the nontrivial element of Z2 acts on an element g e 
\5{\), g changes to g"'. 

In the case of a spin- 1 ferromagnetic phase, the order parameter - ( 1 , 0, 0)^ is invariant 
under a spin rotation about the z-axis followed by a gauge transformation of the same amount as 
the rotation angle; g'^g-'f-^^fefo _ ^fen-o -pj^jg ^ spin-gauge coupled U(l) symmetry denoted 
as // = U(l)f_+0, where the subscript F, + (f)\s shown to indicate the nature of the coupling. A 
general form of the spin-1 ferromagnetic order parameter is given by 



■ferro 







V^cos^f' 





- g'(0-r) 


^sin^ 


loj 




, e'" sin^ f , 



(319) 



Note that the linear combination - y in Eq. ( |319t manifestly exhibits a spin-gauge symmetry, 
i.e., the equivalence between phase change and spin rotation. Individual configurations of ^f'^™ 
are completely specified by the entire range of Euler angles {a,l3, j - (p), and therefore, the order 
pai^ameter manifold is given by 



_ SO(3),xU(l), ^ 



(320) 



7.7.2. Symmetry of order parameter and stationary states 

There exists a close relationship between the symmetry of an order parameter and a stationary 
state of the free-energy functional. To illustrate this, let us consider a smooth real function 
f(x^ + y^) on manifold M^. The set of operations that leave / invariant constitute the group 
G - SO(2) of rotations about the z-axis, where the origin x = y = is the only point that has 
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SO(2) symmetry as an isotropy group, whereas any other point traces a circle, which is called 
an orbit, because every element of G acts on that point (see Fig. l24l l: evidently, the isotropy 
group of each point on a circle consists only of the identity element of G. The union of all orbits 
that have mutually conjugate isotropy groups constitutes what is referred to as a stratum. In the 
present example, all orbits of the concentric circles constitute one stratum, and the origin forms 
the other isolated stratum because the isotropy group of the origin, which is SO(2), is different 
from those of other neighboring points. Now, the fact that the origin is an extremum (maximum, 
minimum, or saddle point) of the function f(x^ + y^) suggests that an orbit that constitutes an 
isolated stratum gives a stationary orbit of a real smooth function, as illustrated in Fig.l24l 

W (b) /A 




Figure 24: (a) Orbits of f(x^ + y^) on a manifold K^. The union of all concentric circles constitute one stratum, and the 
origin forms an isolated stratum, (b) The origin x = y = is isolated in its stratum, and it is an extremum of any real 
smooth function f{x- + V"). 



Michel generalized this idea by subs tituting M? with an arbitrary manifold M, and SO(2) with 



an arbitrary compact Lie group G 111 1511 . Michel showed that if the orbit G(^) is isolated in its 
stratum S {(), i.e., if there exists a neighborhood U of G(^) such that U r\S(() - G(^), then G(^) 
is a stationary orbit of a smooth real function on M that is invariant under G. Here, stratum S (f ) 
is the union of all orbits that share the same type of isotropy group, that is, ^ and ^' belong to the 
same stratum if and only if G^ and G^' are conjugate. 

M ichel ' s the orem has bee n uti lized to find the ground states of su perfluid H e- 3 jl 16ll . p- 
wave 1 117 , 118 1 and li-wave 1 105 1 superconductors, and spinor BECs 1 1 19 , 120ll . In the case 



of a spin-1 BEC, we can uniquely determine the orbits of polar and ferromagnetic phases from 
their isotropy groups. For example, = (1,0, 0)^ and ^o'""^ - (0, 1,0)^ are the only solutions 
of e''^e-''^^'*^ = ^ and e"'^^"^^ = ^, respectively, indicating that stratum 5 (^f™-°.P°i=^) consists of 
(-^(•^feno.poiar-) ^jone. Therefore, the ferromagnetic and polar states are always extrema of the 
mean-field energy functional. The order parameter for such states does not depend on parameters 
that appear in the energy functional, and these states are called inert states. An inert state is robust 
against a change in the interaction parameters, implying that it is a stationary point of each term 
in the energy functional. In superfluid He-3, A, Ai and B phases are inert states. 

It is possible, however, that the energy functional has other orbits of extrema that depend 
explicitly on the coefficients involved in the functional. The corresponding states are called 
non-inert states. Michel showed that the gradient of a G-invariant function is tangential to the 
stratum, implying that if we choose one stratum and find an orbit in it that minimizes the energy 
functional, then the orbit is stationary in the entire manifold. In contrast to inert states, the order 
parameters of non-inert states are determined by the competition among several terms in the 
energy functional. 

In the case of spin-1 spinor BEC, non-inert states appear in the presence of an external mag- 
netic field, where the entire symmetry is reduced to Gb = U(l);r^ x 11(1)^,. The mean-field energy 
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per particle for a spin-1 BEC in a uniform system is given by 



= ^co« + ^cin|/|2-p^m|^j2^<?^m2|^„|2, (321) 

m m 

where / = Cn^mn^n IS the Spin expectation value per particle. Consider an order parameter 

^^P = (cos;r,0,sin;t-)^ (322) 

where < x < ^ and x ^ ■'r/2. Different values of x lead to different orbits, whereas the 
relative phase between and the global phase may vary in each orbit. The rotational symmetry 
about the magnetic field is broken, and the isotropy group of is given by H^^ - (C2)f,0 - 
{1,6'" 6^'^="} c Gb- Substituting in Eq. ( I321l i and minimizing e with respect to x, we obtain 
the following stationary state: 

and the corresponding energy functional 

e[(^^]^^con--f^+q. (324) 

2 4-C]!! 

The AF state corresponds to the spin part of the A2 phase in superfluid He-3 where the magnetic 
field changes the population of up-spin and down-spin Cooper pairs. 
The energies for polar and ferromagnetic states are given by 

= ^con, (325) 

e[^''n^\(co+ci)n-p + q. (326) 

Comparing Eqs. ( 1324b . (I325I I. and (1326b . we can obtain the phase diagram shown in Fig.[T] 

While symmetry considerations are helpful to find stationary states, they elude those states 
that have no remaining symmetry. As an example, let us consider a broken-axisymmetry state ifioll 
whose C2 symmetry axis is not parallel to the external field. This state cannot be obtained from 
symmetry considerations because Gb is completely broken. Nonetheless, we can obtain a state 
that has the same symmetry (as a subgroup of G) as the real ground state from the symmetry 
consideration. 

7.1.3. Procedure to find ground states 

According to the discussions in the preceding subsection, the ground state that has remaining 
symmetries can be found by the following procedure 1.121,1 : 

1 . List all subgroups H of G; 

2. for each H, find an order parameter invariant under h e H; 

3. if a non-inert state is obtained, then find the parameters that appear in the order parameter 
so as to minimize the mean-field energy; 

4. find the lowest-energy one among the obtained states. 
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P 

c\n 



< 1, 



(323) 



Here, we present a list of the subgroups of SO(3) and show how to carry out step 2 of the 
above procedure. The subgroups of SO(3) are well known and they are classified into continuous 
and discrete groups, as listed below. 



• Continuous symmetries 

The only continuous subgroup of SO(3) is U(l). Let us choose the z-axis as the symmetry 
axis. An infinitesimal transformation h e H can be expressed as a combination of spin 
rotation and gauge transformation: 



where 6r and 6^ are infinitesimal real values. The order parameter that is invariant under h 
is an eigenstate of f,. Conversely, the eigenstate of f^ with eigenvalue m has U(l) symme- 
try and it is invariant under a spin rotation e^'^-'f' accompanied with gauge transformation 
e'""^. For a spin-/ system, there are 2/ H- 1 states that have the U(l) symmetry. We denote 
this symmetry as U(l)/r +,„0, where the subscript is shown to explicitly indicate the sym- 
metry between the spin and the gauge degrees of freedom. The combined symmetry for 
;« states is called continuous spin-gauge symmetry, and it relates spin textures to the 
superfluid current, as discussed in Sec.|6] 

• Discrete symmetries 

The elements of the point groups of SO(3) can be described with a combination of the 
following generators: 

C„ : the cyclic group of rotations about the z-axis through angle InUfn with k - 1 , ■ ■ • ,n — 
1. The group is isomorphic to Z„ and it has generators {C,^). 

Dn : the dihedral group generated by C„, and an additional rotation through n about an 
orthogonal axis. The group is isomorphic to Z„ x: Z2 and it has generators {C„z, Cix). 

T : the point group of the tetrahedron whose generators are given by {C2z, Cj,^x+y+z)- 

O : the point group of the octahedron whose generators are given by {C4,, C^^^+y+z^ C2,x+y}- 

Y : the point group of the icosahedro n com posed of 12 fivefold axes, 20 threefold axes, 
and 30 twofold axes. It is shown jl 19ll that a state with icosahedron symmetry does 
not exist for / < 4. 

Here, we specify the symmetry axis for convenience. The generator C„n denotes a 2n/n 
rotation around the direction of Q. In the present case, C„si can be expressed in terms of 
the spin operator f as 



which is described by a (2f+ 1) x (2f+ 1) matrix. Combined with a gauge transformation, 
the state invariant under the transformation h = C„sie'^ can be obtained by solving the 
eigenvalue equation 



h = e'*^e-'f='*- ^ 1 + i6^ - if,6z, 



(327) 




(328) 



(329) 



If the eigenvalue is 1, the invariance of the eigenstate is satisfied by a spin rotation alone, 
and therefore, the gauge symmetry is completely broken. On the other hand, the eigenstate 
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with an eigenvalue e"'** 1 is invariant under a simultaneous discrete transformation in 
spin and gauge. Such a symmetry is called discrete spin-gauge symmetry. What we need 
to do is to find simultaneous eigenstates of a set of generators for all discrete symmetry 
groups. 

The results for the spin-1 BEC are summarized in Table |5] and the structure of each order 
parameter is shown in Fig.l25](a). 



phase 


H 


generators 






e 


F+ 


U(l)f_+^ 




(1,0,0) 




\{cq + ci)n - p + q 


F_ 






(0,0,1) 




j(co + ci)n + p + q 


P 


(Z)oo)f.,0 


{e-'f^te"'C2.] 


(0,1,0) 




jcon 


AF 








Ci/i' ' y c,n j 





Table 5: Symmetry and order parameter structure for the spin-1 spinor BEC, where H is the isotropy group, ^ is a 
representative order parameter, and e is the corresponding energy of the system. 



7.7.4. Symmetry and order parameter structure of spin-2 spinor BECs 

It is straightforward to carry out the procedure for spin-2 spinor BEC. The mean-field energy 
for a spin-2 spinor BEC per particle is given by 

e[(] = \can + ]^cMf? + \cin\A\!- - pY^m\U^ + qY^m\j\ (330) 

m m 

where A - (2^2^-2 - + (1)1 VS is the spin-singlet pair amplitude. We first consider 

stationary states that have isotropy groups H e G. Among them, the state with H c Gb can 
survive as a stationary state in the presence of an external field {p + Q and q + 

• Continuous group U(l): 

There are five states that have continuous symmetries: 



F2+ 


(1,0,0,0,0)^ 


77f2+ 


= U(l)^^,2^ : 


{e>^^e-''^% eF2, -- 


= —2p + 4q + 2cin, 


(331) 


F2_ 


(0,0,0,0, 1)^ 


7/f2- 


= m)F,-2^ ■■ 


{e-'^-^g-'t*), eF2_ 


— 2p + Aq + 2cin, 


(332) 


FU 


(0, 1,0,0, 0)^ 




= U(l)^^^^ : 


e'^e-'"^% fFi, = 


1 

-P + 1+ -jCxn, 


(333) 


Fl_ 


(0,0,0, 1,0)^ 




= m)F,-^ ■■ 


e-'^e-''^% ~epi_ = 


1 

p + q+ 2*^1"' 


(334) 


UN 


(0,0, 1,0, 0)^ 




= (D^)f, : {e 


-''-'te'^Ci.}, euN 


1 

= C2«, 

10 


(335) 



where we indicate the generators of the isotropy group in curly brackets { ) and e - e - 
C()n/2. These states are stationary in the presence of an external magnetic field, where the 
isotropy group //un is reduced to \5{\)p . 



80 




Vtttr -f"" 




Figure 25: Order parameters for (a) / = 1 and (b) / = 2 spinor BECs. Shown are the surface plots of the ampUtude 
I Y,m i/'mi7.m(s)l, where >7,m(s) is a rank-/ spherical harmonics in spin space. The gray scale on the surface represents 
the phase. The continuous and discrete symmetry axes are indicated along with the gauge transformation factors. 



• Tetrahedron group T: 
We solve the simultaneous eigenstate of C2z and C^x+y+z' whose matrix forms are given by 



n o^ 

0-1000 

10 

0-10 

10 1. 
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.x+y+z 
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-2i 


2 
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-li 


V6 


2i 


-1 . 



(336) 



The only solution is = ^(1,0, / V2, 0, 1)^ (up to the phase factor). Here, satisfies 
Ciz^c = and Cxx+y+z(c - e'^^^^^c- Therefore, the generators for the tetrahedron group 
as a subgroup of G are C2z and e''"^" ^^Cz, x+y+z'- 

C :i(l, 0, / V2, 0, 1)^ He = Tf,^ : {C2,, e-'^'^'^Cxx+y^z), ec = 2q. (337) 

This state is called a cyclic state. Cyclic states need not be stationary in the presence of an 
external field because T ^G-q. 
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• Dihedral groups D„: 

In a manner similar to the tetrahedron group, we find 

BN : J^(l, 0,0,0, 1)^ //bn = (D4)f„^ : {e"'C4„C2.}, ^bn = 4^ + (338) 



D2 



V2 
V2 



, 0, sin;^', 0, 



cos;t' 



, Hu2 = (-02)^ : {C2z, CiA, ~em = 4^cos2;^r + 
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(339) 



D2' : - 
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1 i-,0,i\ 2 1 + — ,0, ,/l - 



C2n 



2 



Hi 



D2' 



(Z)2)f : {C2,-, C2x), em' = 2^ 



lOg 



(340) 



There are two energy minima D2 and D2' in the order parameter manifold with an isotropy 
group £>2- When p - q - Q, UN, BN, and D2 states are all degenerate, and;^' in D2 state 
can take any arbitrary real value, while D2' phase goes to the cyclic state. In the presence of 
an external field, however, BN, D2 and D2' states need not be stationary because £)„ (J: Gb- 

• Cyclic groups C„: All states that have cyclic symmetries are non-inert states. Here, we 
have chosen the order parameter to minimize the mean-field energy. 
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C2' : ^ (fl, 0, b, 0, c)^ , Hc2' = (C2)f. : {C2,). 
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(342) 

(343) 

(344) 
(345) 



The order parameter and energy for C2' state are the same as those of C2' state in Table |4] 
The obtained results are the same as those listed in Table |4] The structure of each order 



parameter is shown in Fig.l25](b). 
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7.2. Homotopy theory 

7.2.1. Classification of topological excitations 

Bose-Einstein condensates can accommodate topological excitations such as vortices, monopoles, 
and Skyrmions. These topological excitations are diverse in their physical properties but have 
one thing in common; they can move freely in space and time without changing their charac- 
teristics that are distinguished by topological charges. The topological charges take on discrete 
values and have very distinct characteristics independently of the material properties. It is these 
material-independent universal characteristics and robustness to external perturbations that make 
topological excitations unique. 

The classification of topological excitations is best made by homotopy theory il22l Il23ll . 
This theory describes what types of topological excitations are allowed in what order-parameter 
manifold. For example, n2{U{l)) - implies that monopoles, which are characterized by the 
second homotopy group n2, are not topologically stable in systems described by scalar order 
parameters. This theory also describes what happens if two defects coexist and how they coalesce 
or disintegrate. 

By way of introduction, we first consider a BEC described by a scalar order parameter \li{r) — 
\il/ir)\e"^^''\ To classify line defects such as vortices, we take a loop in the condensate and consider 
a map from every point r on the loop to the phase 0(r) of the order parameter i}/: 

r i-> (j)(r) = argi^(r). (346) 

If the loop encircles no vortex (e.g., loop A in Fig.l26]l, the image of the map covers only a part 
of the unit circle, as shown in the right-hand side of Fig.|26l If the loop encircles a vortex (e.g., 
loop B), the image covers the entire unit circle at least once. If it covers the circle n times, it is 
said that the winding number of the vortex is n. The crucial observation here is that the winding 
number does not change if the loop deforms or moves in space as long as it does not cross the 
vortex. This property can be used to classify loops according to their winding number. 



vortex 




Figure 26: Mapping from a loop in real space onto the order-parameter manifold S (unit circle) according to the 
correspondence 0(r) : r — > argi^(/-). This mapping defines the first homotopy group ni(S^). 

If two loops {a and /'/, are continuously transformable without crossing singularities of the or- 
der parameter, they are said to be homotopic to each other and written as {„ ~ (b- The homotopic 
relationship is an equivalent one in that it is symmetric (i.e., ta ~ (a), reflexive (i.e., if {„ ~ (t. 
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then (h ~ ig), and transitive (i.e., if ~ (b and lb ~ (c^ then la ~ (c)- By this equivalence 
relationship, all loops are classified into equivalent classes called homotopy classes: 

[A], [^2], [^3],--- , (347) 

where (/ = 1, 2, 3, ■ ■ ■ ) is an arbitrarily chosen loop from the homotopy class [/",] because all 
loops belonging to the same class are equivalent and continuously deformable to each other. 

Mathematically, a loop I is defined as a mapping from / = [0, 1] to a topological space such 
that £iQ) - {(\) = xq, where jcq is called a base point. If two loops £„ and £h share the same 
base point xq and they are continuously deformable to each other, they are said to be homotopic 
at If the base point is not shared, they are said to be freely homotopic. The constant loop 
c is defined as the one such that c(f) = xq for e [0, 1]. The inverse of loop I is defined as 
ri(0 = ^(l-f)for''f6 [0,1]. 

The product of two homotopy classes {{\\ and {(2} is defined as 

[tx]-[h]^[(i-ti\, (348) 

where li ■ I2 denotes the product of two loops in which i\ is first traversed and then I2 is traversed. 
Members of \i\ ■ £2] that are homotopic to £1 ■ I2 need not return to the base point xq en route to 
the terminal point (see the dashed curve in Fig.l27b. 




Figure 27: The product of two loops ([ ■ (2 in which ([ is first traversed and then (2 is traversed. Loops homotopic to 
tl ■ t2 need not return to xq en route to the terminal point, as shown for the dashed curve. 

With the definition of the product ( 1348b . a set of homotopy classes form a group. In fact, they 
satisfy the associative law 

iVi] ■ [fi]) ■ [^3] = [A] ■ ([ii] ■ [^3]); (349) 
those loops that are homotopic to the constant loop c constitute the identity element [c] such that 

[c] ■ M = M ■ [c] forany (350) 
finally, the inverse [^"'] of [£] is the homotopy class that consists of inverse loops of [I] so that 

[r'] -m^m ■[£-'] ^[c]. (351) 
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The group defined above is called the fundamental group or the first homotopy group and it is 
defined as 7ri(M, xq), where xq is the base point of all loops. In many situations in physics, the 
choice of the base point is not important, and we shall omit it and write 7:i(M,xq) simply as 
7ri(M). 

The fundamental group is said to be Abelian if any two elements commute (i.e., [m] ■ [n] = 
[n] ■ [m] for ^[m] and ^[n]) and non-Abelian if some elements do not commute. For the case of 
U(l) vortices, the fundamental group is Abelian and ■ is simply addition. For example, if two 
singly quantized t/(l) vortices coalesce, a doubly quantized vortex results. This can be described 
as [1] ■ [1] = [2]. Thus, the fundamental group is the additive group of integers Z: 

7ri(5') = Z. (352) 

Point defects such as monopoles can be characterized by considering a sphere E that covers 
the object (see Fig. l28T l. We consider a map from each point r on E to a point m in the order 
parameter space M: 

m-.I.^ M. (353) 

We classify elements of such a map by regarding as equivalent any two elements that can trans- 
form into each other in a continuous manner For the case of the monopole, we consider M to 
be a two-dimensional sphere or 2-sphere S^, as shown in the right-hand side of Fig. |28] In this 
case, if the map (I353l l wraps n times, we have 712(8^) = [n]. 




Figure 28: Mapping from a spliere in real space onto the order-parameter manifold M according to the correspondence 
'F(r) : r —> M. This mapping defines the second homotopy group n2(M). Here, we take M to be — a two-dimensional 
unit sphere (2-sphere). 

Skyrmions and knots are classified by the third homotopy group. We assume that the order 
parameter becomes uniform at spatial infinity. Then, the three-dimensional space is compacted 
into a three-dimensional sphere . To help envisage it, imagine a two-dimensional plane and 
identify all infinite points. Then, the two-dimensional plane is compacted into a two-dimensional 
sphere 5^. By considering the map from to the order parameter space M and by identifying 
maps that can be continuously transformed into each other without crossing any singularity of 
the order parameter space as belonging to the same equivalent class, we can introduce elements 
of the third homotopy group nT,{M). 

Higher homotopy groups can be defined in a similar manner Moreover, the zeroth homotopy 
7To(M) classifies the domain walls and gives the number of disconnected domains; in particular, 
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7Tq(M) - implies that the order-parameter manifold is connected. Using homotopy theory, 
one can elicit insights into the topology of complicated structures, which is difficult to obtain 
intuitively. Here, we list some useful formulas of homotopy theory: 



7T„iS'") = 



( z 


if m 


= n > 1; 





if m 


> « > 1; 





if m 


= 1 and n > 2; 


z 


if n 


- 3 and n -2 (Hopf charge); 


I Z2 


if n 


- m + \ > 4 or n = m + 2 > 4, 



(354) 



where implies that trivial configurations (i.e., those contractible to a point) exist and Z2 = {0, 1 ) 
is the two-element group. 

Let G be a Lie group and H be its subgroup. If G is connected (i.e., 7ro(G) - 0) and simply 
connected (i.e., n\{G) - 0), the following isomorphism holds: 

niiOIH) - no{H). (355) 

Similarly, if n2{G) - 7Ti(G) = 0, then 

7T2(G/H) s 7Ti(H). (356) 

As a corollary, when G - SU(2) and // = U(l), we obtain 

7r2(SU(2)/U(l)) = 7ri(U(l)) = Z. (357) 

Some order-parameter spaces R and their first, second, and third homotopy groups are sum- 
marized in Table |6] The topological objects classified by homotopy groups are summarized in 
Table|7] For example, in the case of ferromagnetic BECs, tti (S0(3)) = Z2 implies that there exist 
two types of linear defects that are singular and nonsingular; 7r2(SO(3)) = implies the absence 
of point defects and 2D Skyrmions; and 7r3(SO(3)) = Z implies the presence of a nonsingular 
soli ton-like object such as a Shankar monopole (Skyrmion) 1 124, 125]. 

Here, we define a Skyrmion as a nonsingular point object characterized by a 712 charge in two 
dimensions and by a n^, charge in three dimensions. In the polar phase of a spin-1 condensate, a 
7T2 Skyrmion takes a configuration in which the director points in the positive z-direction at the 
center, flares out in the radial direction, and eventually points in the negative z-direction as one 
goes away from the center On the other hand, a 7:3 Skyrmion is exemplified by a knot soliton, 
as discussed in Sec. 17.3.41 In contrast, a monopole is a singular point object in three dimensions 
(a singular point object in two dimensions is called a vortex). The charge of the monopole is 
characterized by a 712 charge; examples include a hedgehog ('t Hooft-Polyakov monopole) in 
a polar condensate (see Sec. I7.3.2l i. Note, however, that the 712 monopole cannot be created 
in a ferromagnetic condensate because its order-parameter manifold is SO(3), which does not 
give a non-trivial 712 charge; in this case, the monopole must be attached by a Dirac string (see 
Sec. 173:21) . 



7.2.2. Relative homotopy groups 

To fully understand broken-symmetry systems, it is important to take into account not only 
the global boundary conditions but also various constraints imposed on the order-parameter 
field. Relative homotopy groups provide a useful tool for characterizing effects when the order- 
parameter manifold changes its character on a surface or when it is constrained by boundary 

86 





R 




^2 ^3 


scalar BEC 


U(l) 


Z 





Heisenberg spin 


52 





z z 


nematics 




Z2 


z z 


biaxial nematics 


SU(2)/Q 


Q 


z 


ferromagnetic BEC 


SO(3)f,^ = RPl,^ 


Z2 


z 


spin-1 polar BEC 


[5|xU(l)^]/(Z2)ir,^ 


z 


z z 


cyclic BEC 


[SO(3)f X U(l)^]/rF,^ 




z 



Table 6: List of homotopy groups. Subscripts F and denote the manifolds for the spin and gauge degrees of freedom, 
respectively. Q and T* respectively denote the quaternion group and the binary tetrahedral group, the latter being a 
subgroup of SU(2)xR. 





defects 


solitons 




domain walls 


dark solitons 




vortices 


nonsingular domain walls 




monopoles 


2D Skyrmions 






Skyrmions, knots 



Table 7; Topological objects (defects and solitons) described by homotopy groups. 



conditions. Relative homotopy groups can also be used to discuss properties of a defect that is 
itself regarded as a singularity in the order-parameter field in ordinary homotopy theory. Finally, 
the relative homotopy group is useful to calculate the absolute homotopy group via the exact 
sequence of homomorphisms. 

As an illustrative example, let us again consider a ferromagnetic BEC. Suppose that the 
system is contained in a cyUnder of radius ro and impose a boundary condition that the spin 
points in the positive z-direction at r = and outward or downward on the wall at r = ro, where 
r is the radial coordinate. Then, the spin texture would look like 



^ sin^r) 



(358) 



Because the direction of spin is s{r, ip, z) - (e^ cos (p + Cy sin (p) sin /3{r) + cos j3(r), the above 
boundary conditions are met if we choose j3{0) = and y6(ro) = 7r/2 or /3(ro) = n. The topological 
objects corresponding to jS(ro) - n/2 and y8(?"o) = ;r are referred to as the Mermin-Ho vortex and 
Anderson-Toulouse vortex, respectively. It can be shown that the circulation at the wall is 1 
for the Mermin-Ho vortex and 2 for the Anderson-Toulouse vortex. They are good examples of 
relative homotopy because both of them are not stable in free space and they can be stabihzed 
only if the boundary conditions on the wall are imposed by some means. 

In general, the nth relative homotopy group nn(R,R) consists of the elements of 7t„(R) from 
which those of n„{R) are subtracted. We consider an n-sphere in R and let it expand into R by 
continuous deformation. Such an expansion may be regarded as a one-to-one mapping of nniR) 
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into a subgroup of 7r„{R). Let Im(7r„(^) 7t„{R)) be the image of this homomorphism. Then, 
by the definition given above, the relative homotopy 7t„(R, R) is a quotient group of Im(7r„(^) — > 
n„(R)) in TTniR): 



JT„(R,R)^ 



n„(R) 



(359) 



Im(7r„(^) ^ TTniR)) 



The relative first homotopy group 7Ti{R,R) can cope with planar defects with nontrivial 
boundaries. A well-known example is a planar soli ton (Maki domain wall) in the A phase of 
superfluid helium-3, where the order-parameter manifold is R = (S^ x SO(3))/Z2 within a dipole 
healing length from the domain wall and R = SO(3) outside this layer due to dipole locking. 
The domain wall may be regarded as a linear defect as seen from the edge of the wall. We 
take the axis perpendicular to the domain wall as the z-axis. Then, the order-parameter field 
defines a mapping from the z-axis into the contour in R whose end points belong to R. Because 
7ri(SO(3)) = Z2, such contours are also classified as Z2, of which the nontrivial one corresponds 
to the Maki domain wall. 

The relative second homotopy group 7Ti{R, R) allows the stabilization of a point object subject 
to a boundary condition on a sphere enclosing it. Another example is the boojum-a point defect 
that appears in the A phase of superfluid helium-3 when it is contained in a hard-wall sphere. 
Because ;7r2(SO(3)) = 0, the bulk superfluid A phase does not support a topologically stable point 
defect. However, because the I vector must be perpendicular to the surface, a point singularity 
for the /-vector appears on the surface that is connected to a nonsingular vortex texture in the 
bulk. 

7.3. Examples 
7.3.1. Line defects 

The line defects are characterized by the first homotopy group (or the fundamental group) 
7[\{M). For the case of i-wave superconductors, liquid helium-4, and spin-polarized or spinless 
gaseous BECs, the order parameter is scalar and the order-parameter manifold is M = U(l). The 
fundamental group is therefore the additive group of integers: 



For the case of a spin-1 BEC, the ground state can be ferromagnetic or polar. The order- 
parameter manifold of the ferromagnetic phase is SO(3) and the fundamental group is ;ri (SO(3)) - 
Z2 - {0, 1 ), the two-element group. The spin configuration that corresponds to class is continu- 
ously transformable to a uniform configuration, whereas the spin configuration that corresponds 
to class 1 describes a singly quantized vortex that is topologically stable. Because 1 -1- 1 = in 
Z2, the coalescence of two singly quantized vortices is homotopic to a uniform configuration. 
Conversely, there is a configuration called the Anderson-Toulouse vortex that is nonsingular at 
the origin but has circulation of two away from the origin. 

Next, let us consider the polar phase of a spin-1 BEC. The order-parameter manifold is 
^poiar _ |-^2 ^ U(l)0]/(Z2)/r,^. As wc havc discusscd in Sec. 17.1.11 the Z2 symmetry arises 
from the fact that the order parameter given by Eq. ( 1316b is invariant under the simultaneous 
transformation of (0, d) ((p + n, -d). Because G = xU{l) is not simply connected, to cal- 
culate ni(RP°^^'), we hft G to G' = 5^ x M and use Eq. ( 1355b . Then, the points that are identical 



;ri(U(l)) = Z. 



(360) 
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to (0, d) in G' are {(f> + rm, d) for even n and (0 + nn, -d) for odd n, and those points constitute 
an additive group of integers n. Hence, we have 

7ri(/;P°'^') - Zf,^. (361) 

As in the case of a scalar BEC, vortices in the polar phase are classified by integers. However, 
due to the Z2 discrete symmetry, the minimum unit of circulation is one-half of the usual value 
of h/M, as discussed in Sec. 16.2.21 

In a similar manner, we may expect that the uniaxial, biaxial, and cyclic phases in spin-2 
spinor BECs can support various types of vortices due to the discrete symmetry, as discussed in 
Sec. |6] Some of them have fractional circulations, and some others have no mass circulation. 

The topological charge can be expressed algebraically as a functional of the order-parameter 
field. The topological charge corresponding to 7ri(5 ') can be expressed as a line integral of the 
velocity field. 



7.3.2. Point defects 

Point defects are characterized by the second homotopy group. Because 7r2(U(l)) = 0, scalar 
BECs have no topologically stable point defects. When the order-parameter manifold is S^, it 
can host monopoles because 7:2(8^) = 1. 

Here, we consider the situation shown in Fig.|28j where m{r) is a mapping from S^XaS^- . The 
topological charge of the point defect can be calculated as follows. By definition, A^2 gives 
the number of times the vector m{r) shown in Fig. |28] wraps 5^. Expressing the components of 
this vector in spherical coordinate as m ^ - sin a cos ji, niy = sin a sin /3, and - cos a, we can 
calculate N2 as 



A^2 = — I dQ. 



0(6, i 



(362) 



where d£l = sin 0d6d<p, with and <p are the spherical coordinates of r, and the last term on 
the right-hand side is the Jacobian of the transformation of the coordinates. The fact that the 
integrand is the Jacobian explicitly indicates that N2 gives the degree of mapping from in real 
space to 5 ^ in order-parameter space. The right-hand side of Eq. ( 13621 ) can be directly expressed 
in terms of m as 

where 

j = ^e,ji'«/(Vmj X Vmt). (364) 
We can use Gauss' law to rewrite the right-hand side of Eq. ( I363l l as a volume integral: 

N2 = j n^dr, (365) 



where 



nUr)^^V-j(r) (366) 
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gives the density distribution of point singularities. 



't-Hooft-Polyakov monopole (hedgehog) 

Next, we consider the polar phase. The second homotopy group is given by 7T2(R^°^'") = Z, 
and therefore, it can accommodate point defects. To investigate their nature, we write the order 
parameter in the following form: 



^1 



polar 



V2J, 



(367) 



By setting 
hedgehog: 



0, we obtain the spherical monopole called 't-Hooft-Polyakov monopole or 



d{r) 



(368) 



It follows from Eq. ( 1363b that the topological charge of the hedgehog is A^2 - 1 ■ Moreover, j 
defined in Eq. ( I364l i is related to the circulation of spin currents given by Eq. (I287l i as 



X V 



(spin,P) 



(369) 



which implies that the surface integral of the spin circulation is quantized in terms of 4h/M; 

[j]d,{Vx vf ■ dS = -^N2. (370) 



AUce ring 

The m onopo le in the spin-1 polar phase is energetically unstable against deformation into an 
Alice ring [126]. The Alice ring is a combined object of two defects characterized by tti and 712 
(see Fig.|29]l. Viewed far from the origin, it appears to be a point defect (A^2 = 1); however, along 
the ring C on which the order parameter is singular, it appears to be a continuous distribution of 
line defects (A^i = 1). This topological object is called the A lice ring, and it has been predicted 
to be realized in an optically trapped spin-1 -^Na BEC fl2i^. 

Dirac monopole 

The magnetic monopole was originally envisaged as a magnetic analogue of the quantized 
electric charge: 

V ■ B ^ 47Tg6(r), (371) 

where g denotes the strength of the magnetic monopole. The solution of this equation is easily 
found to be B = gr/r^ and the corresponding vector potential is given by 

" I m g(l +cos6>) 

(-3',ji:,0) = — — e^, (372) 



r(r - z) r sin Q 
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Figure 29: Alice ring comprised of continuously distributed half-quantum vortices along a contour C. Far from the 
origin, it appears to be a monopole. 



where (r, 6, ip) are the polar coordinates and is the unit vector in the (/^-direction. The vector 
potential (I372I I reproduces the magnetic field (1371b . except on the positive z-axis along which the 
magnetic field exhibits a singularity called the Dirac string: 



rotA = g-r - Ang6{x)6{y)6{z)e„ 



(373) 



where e, is the unit vector along the z-direction. 

The Dirac monopole can be created in the ferromagnetic phase of a spin-1 BEC 112711 . Sub- 
-(f, a ^ (fi, and /J — > 6* in Eq. (1319b . we obtain 



stitutmg I 



^Dirac 



fOirac 



/-V 



fe-^'-f cos2 

V2 

sin^ |. 



(374) 



The corresponding superfluid velocity v^™"'*^' - (ft/M)Im(^' V^) [defined by Eq. (I259b l is the 
same as the vector potential ( 13721) of the Dirac monopole with the identification g = h/M. 

It is interesting to note that the order parameter (1374b exhibits a doubly quantized vortex only 
along the positive z-axis. One can say that because the first homotopy group of SO(3) is Z2, 
which allows only a singly quantized vortex called a polar-core vortex as a stable topological ob- 
ject, the doubly quantized vortex must terminate at the magnetic monopole. In fact, the singular 
Dirac monopole can deform continuously to a nonsingular spin texture if we take the following 
parametrization: 



1-Diiac 



V2 



sin/? 



sm 



/3(r,<p,z) = 6(l-t)+7Tt. 



(375) 



We can see that as t increases from to 1 , the order parameter continuously changes from the 
Dirac monopole at f = to a nonsingular texture ^ = (0, 0, 1)^ at f = 1. 



7.3.3. Action of one type of defect on another 

When two different types of topological objects coexist, there exists an action of one object 
on the other as they move relative to each other As an example, consider a spin-1 polar phase 
in which two monopoles with charge 1 merge in the presence of a half-quantum vortex. As a 
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monopole turns around the vortex, it changes sign. Thus, if they merge along hne 1, the order- 
parameter field must deform, as shown in Fig.|30](b), and the combined object has the topological 
charge of two. On the other hand, when they merge along line 2, the order-parameter manifold 
can deform into the one shown in Fig.[30](c), and the topological charge is zero, indicating that 
the two monopoles annihilate in pairs. This example illustrates that the charge of the monopole 
changes its sign as it moves around the half-quantum vortex along curve 2. 



half-quantum vortex 




Figure 30: (a) Coalescence of two monopoles with charge one in the presence of a half-quantum vortex, (b) If they merge 
along the dashed curve 1 shown in Fig. (a), the topological charge of the combined object is two. (c) If they merge along 
the dashed curve 2, the topological charge of the combined object is zero. 



7.3.4. Skyrmions 
Shankar Skyrmion 

The third homotopy group characterizes topological objects called Skyrmions that extend 
over the e ntire three -dimensional space. A prime example of this is the so-called Shankar 



monopole 11124 112511 . although it is actually not a monopole but a Skyrmion. The topology 



of the Shankar Skyrmion is n-3(SO(3)) = Z, which is supported by the ferromagnetic phase of a 
spin-1 BEC. The order parameter of the ferromagnetic BEC is characterized by the direction of 
the local spin s and the superfluid phase that can be parametrized in terms of two unit vectors 
such as i and u that form a triad with s: s - i xu. Concrete forms of s and t are introduced in 
Sec. ISO 

The class one e lement of 7r3(SO(3) ) = Z may be realized by rotating the triad at position r 
about the direction llT2llT24l[l2ilT29ll 



n=-f(r)N3, N^eZ (376) 
r' 

through an angle /(r)A^3, where /(O) - 2ji and /(oo) = 0. It follows from the last condition that 
the triad is uniform in spatial infinity, and thus, the three-dimensional space is compacted to 5^. 
The case of n = 1 is shown in Fig. [311 

The topological charge of the third homotopy group can be introduced in a manner similar to 
that of the second homotopy group. Now, the mapped vector has a four-component vector that is 
normalized to unity. One such representation is 

r . fir) fir) 

n = - sm , «4 = cos . (377) 

r 2 2 

The invariant of such a mapping is again given by the Jacobian of the transformations: 

A^3 = J dreafjyseijknadinpdjtiydkns, (378) 
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Figure 31: Shankai' monopole with charge one. 



where e^i- and Sa^ys are the completely antisymmetric tensors of the third and fourth orders, re- 
spectively, and the Roman letters run over x, y, z and the Greek letters run over 1,2,3,4. 

Knot soliton 

The third homotopy group also classifies knot solitons. Knots differ from other topological 
excitations such as vortices, monopoles, and Shankar monopoles in that they are classified by 
a linking number while others are classified by winding numbers. Knots are characterized by 
mappings from a three-dimensional sphere to S^. As in the case of Shankar monopoles, the 
S ^ domain is prepared by imposing a boundary condition that the order parameter takes on th e 
same value in every direction at spatial infinity. Here, we consider the spin-1 polar phase il30|] . 
The order parameter manifold of the polar phase is given by /;p°'^' = (Sg x U(1)0)/(Z2)s,0, of 
which neither U{1) nor Z2 symmetry contributes to the homotopy groups in spaces higher than 
one dimension; hence, we find that 7r3(/?''°'") = 7:^(8^) = Z. The associated integer topological 
charge Q is known as the Hopf charge, and it is given by 

e = ^ jd\ EijtTij^^k, (379) 



where 'Fij - diJlj - djJii - d ■ {did x djd) 11 13 111 . Note that the domain (r) is three-dimensional 



while the target space {d) is two-dimensional. Consequently, the preimage of a point on the target 
constitutes a closed loop in 5^, and the Hopf charge is interpreted as the linking number of 
these loops: if the d field has Hopf charge Q, two loops corresponding to the preimages of any 
two distinct points on the target 5^ will be linked Q times [see Fig.[32](a)]. Figure[32](b) shows 
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an example of the d field of a polar BEC with Hopf charge iQ. 

Knots can be created by manipulating an external magnetic field. In the presence of an 
external magnetic field, the linear Zeeman effect causes the Larmor precession of d, whereas d 
tends to become parallel to the magnetic field because of the quadratic Zeeman effect. Suppose 
that we prepare an optically trapped BEC in the m = state, i.e., d - (0,0, 1)^, by applying 
a uniform magnetic field in the z-direction. Then, we suddenly turn off the uniform field and 
switch on a quadrupole field. Because of the linear Zeeman effect, d starts rotating around the 
local magnetic field, and therefore, the d field winds as a function of t, resulting in the formation 
of knots. 

Figure [33] shows the creation dynamics of knots in an optical trap subject to the quadrupole 
field, where the upper panels show the snapshots of the preimages of d = -z and d - z and 
the bottom panels show cross sections of the density for m - - 1 (bottom) components on the xy 
plane. The density pattern in m-- \ components is characteristic of knots; a double-ring pattern 
appears corresponding to one knot. As the d field winds in the dynamics, the number of rings 
increases. This prediction can be tested by the Stern-Gerlach experiment. 

Topological objects are conjectured to play important roles in the formation of our universe. 
Among them, knot solitons have recently attracted the interest of cosmologists since Faddeev 
and Niemi s ugge sted that knots might exist as stable solitons in a three-dimensional classical 
field theory liill]. 




Figure 32: (a) Preimages of two distinct points on S " form a link, (b) Spin configuration of a Icnot witli Hopf charge 1 
in a polar BEC, where arrows inicate the d field of the polar phase. The solid and dashed lines trace the point where d 
points to X and -z, respectively, and form a link. 



8. Many-Body Theory 

In spinor BECs, the mean-field energy is much larger than the spin-exchange interaction. 
Therefore, we first determine the spatial dependence of the condensate wave function indepen- 
dently of the spin degrees of freedom, and then, we consider the spin state by assuming that 
the obtained spatial wave function is shared by all spin states. This approximation is called the 



'Technically speaking, the configuration in Fig. 1321 (b) is an unknot of a pair of rings with linking number 1, because 
the preimage of one point on S " forms a simple unknotted ring. 
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Figure 33: Dynamics of the creation of knots in a splierical optical trap under a quadrupole magnetic field. Snapshots of 
the preimages of rf = (0, 0, - 1 and d = (1,0,0)^ (top), and cross sections of the density for m = - 1 components on the 
xy plane (bottom). 



single-mode approximation. In this section, we examine the many-body spin states by using this 
approximation. 

Suppose that the coordinate part of the ground-state wave function cpoir) is independent of 
the spin state. Then, (f>o{r) is the lowest-energy solution of the GP equation 



' 2M 



+ t/trap +Co(A^- l)|0ol' 



bo - yU0o, 



(380) 



where f/trap is the trapping potential. Let yUo and fii be the first- and second-lowest eigenvalues of 
Eq. ( 13801) . respectively. The criterion for the validity of the single-mode assumption is that the 
first excited energy /ii - //q is much larger than other energy scales: 



fii-fiQ»\p\,\ci\N/V''', \C2\N/V 



/eif 



(381) 



where = (J dr\(f>()f) ' is the effective volume of the system. When this condition ( 13811 ) is 
satisfied, the field operator i^m (m = /, / - 1, ■ ■ ■ , -/) may be approximated as iff„,{r) ^ ai„(po{r), 
where a„, is the annihilation operator of the bosons that have magnetic quantum number m and 
coordinate wave function Sn. 



8.1. Many-body states of spin-1 BECs 
8.1.1. Eigenspectmm and eigenstates 

As explained above, the field operator of a spin-1 BEC in the single-mode approximation can 
be expressed as lIsslHIsell 

(Am(r) = a„,(t>Q{r) (m = 1,0,-1). (382) 
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For the sake of simplicity, let us consider the case of zero magnetic field (i.e., p - q - 0). 
Substituting Eq. ( 1382b in Eq. (l24l i and replacing 2„ al,a,n and 2,„ „ al^alanam with .^V and N(N - 
1), respectively, we obtain 

H^ffjdr cf,; {-^ + U{r) + 'j-(N - l)|0ol'j 4>o + c[:F^:, (383) 

where 

1 

(384) 

m,n=-l 

and c'j = (ci /2) J t/r|0ol"*; 0o is determined by Eq. ( I381l l subject to the normalization condition 

|0ol'fl'r=l. (385) 



Substituting the solution of Eq. ( 1381b back into ( 1383b . we obtain 

H = fiN - c'qN{N - 1) + c\ : :, (386) 
where s (co/2) J (irl^ol"*- From Eq. ( fTSI ), we find that the total spin operator F satisfies 

■.F^:^N{N-l)-3S^'S, (387) 

where 

5 = ^(fli;-2flifl_i) (388) 
V3 

is the spin-singlet pair operator We use Eq. ( 13871 ) to rewrite Eq. ( 1386b as 

H^^iN - c'qN{N - 1) + c'l [a?(A^ - 1) - 3S'S\ . (389) 

Because .^V and S ^5 commute, the eigenvalue problem reduces to finding their simultaneous 
eigenstates. 

Let |0) be the state in which there is no spin-singlet pair, i.e., S |0) - 0. Using the commuta- 
tion relation 

[5,(5')*] = -k{2N + 2k+\){S'^f-\ (390) 

we obtain 

5"'"5(5"V|0) = ^k{2N -2k+\){S^f\Q). (391) 



We then define the A;-pair state |^): 



1.) . , ^^^'^^l»> (392) 



It follows from the commutation relation ( 1390b that 



S'^S\k) ^^k(2N -2k+l)\k). (393) 
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8.1.2. Fragmentation 

Next, we consider a state \k, d) in which 2k atoms form spin-singlet pairs and the remaining 
N -2k atoms reside in the m = state along the d direction (\d\ = 1). Defining 



+ idy d^ — idy 

= — '■ — -^^ay + dM[) H — '■ — r^— 0-1, (394) 

the corresponding eigenstate and eigenenergy are given by 

\k,d}^zJ(alf-^''(S-^f\0}, (395) 
Ek = ijN-c'qN(N-1)+c\ [N(N-l)-2k(2N-2k+l)], (396) 

where Z,„ is the normalization constant. In the absence of an external magnetic field, all magnetic 
sublevels are degenerate. If we take the quantization axis along the z-direction, we can specify the 
eigenstate in terms of the total spin F, magnetic quantum number F-, and number of spin-singlet 
pairs k: 

\k, F, F.) ^Z-HS ^'f(F-f'''' (fl"[)^|vac), (397) 

where F - N - 2k, and 

F- = ^ (f-)mn aifln = ^{u^ai + a[.^ao) (398) 

m.n 

is the lowering operator of the magnetic quantum number 

Next, we consider a situation in which all F atoms points in the z-direction in addition to the 
k spin-singlet pairs, i.e., F. - F in Eq. (13971 ), and seek the number of the atoms that are present 
in the m - state. We introduce the basis state |ni,«o,n-i) in which the magnetic sublevels 
m= 1, 0, -1 are occupied by ni, no> ™d n_i atoms, respectively. Then, |vac) = |0, 0, 0) and 

^ (vac|fl|^5*a|!)flo('§ ')'^(fl[)''|vac) 



<vac|a[i-'^(i'')'<(fl|)'|vac) 
(F.Q.Q\S''aia,^(S^f\F.Q.Q) 

(399) 



<F,0,0|5*(5')*|F,0,0) 
Substituting 0^30 = - 1 and using the fact that Sg commutes with 5 ' , we have 



(aLao) = — 1- (400) 

" <F,0,0|5*(5')*F,0,0) 



We use Eq. ( |390t to rewrite 



= jk(2N + 2k+l)(S')''-^ +(S^fS. (401) 
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Noting that S\F,l,0}^S \F, 0, 0) = and N\F, 1 , 0) = (F + 1 , 0), we have 



{2F + 2k + 3){F, 1 , 0\S '(5 "0 |f , 1 , 0) 
(2F + 2k+ 1){F, 0, 0|5 *- 1(5 'l^', 0, 0) 

(2F + 2k + 3)i2F + 2k+l)----(2F + 5) 
(2F + 2k+ 1)(2F + 2fe - 1) ■ ■ ■ ■ (2F + 3) 

2k _ N-F 
2F + 3 ~ 2F + 3' 



- 1 



Because (ajfii) - (fl_jfl_i) - F and Yjm=-i{^m^m) - N, we obtain 



ni 
ho 
n-i 



= <fl|fli> = 



NF + F^ + N + 2F 



2F + 3 



N-F 



2F + 3 
(N - F)(F + 1) 



2F + 3 



(402) 
(403) 
(404) 



When all atoms form spin-singlet pairs, i.e., 
sublevels are equally populated: 



F = 0, Eqs. ( I402b -( l404t show that all magnetic 



N 



(405) 



That is, the condensate is fragm ented in th e sense that more than one single-particle state is 



macroscopically occupied 111321 14 l86l 113311 . This is in sharp contrast with the mean-field re- 



sult that predicts that for ci > 0, the ground state is either polar (no = A^) or antiferromagnetic 
(«! - h-i - ^). This discrepancy originates from the mean-field assumption that there exists one 
and only one BEC. However, such an assumption cannot be justified when the system possesses 
certain exact symmetries such as rotational symmetry in the present case. Equation (I403I I shows 
that «() decreases rapidly with an increase in F. In fact, when F - 0{ y/N), we have Hq - 0{ y/N). 
This implies that although mean-field results break down at zero magnetic field, the validity of 
mean-field theory quickly recovers as the magnetic field increases. This affords an example of 
why fragmented BECs are fragile against symmetry-breaking perturbations. 

Finally, we note an interesting relationship between a mean-field state and the many-body 
spin-singlet state. A mean-field state is one in which aU particles occupy the same single-particle 
state: 



|vac). 



(406) 



where Ud is defined in Eq. (I394l i. Substituting dx - sin 6 cos 0, dy - sin 6 sin (p, and d, - cos 6 in 
Eq. (I394I I, we obtain 



1 



a. ,. a , 

sin Oe"'' + a' cos 6 + sin de 

V2 V2 



|vac). 



(407) 
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We calculate the equal-weighted superposition state of \6, (f)) over all solid angles dO. = sin OdOdcf). 



|sym) = — I smOde \ dcp \e. 
47r Jo Jo 



if A? is odd; 

p=(aA^ - 2a'fl' )^|vac) if is even. (408) 

This result shows that the spin-singlet state is the superposition state of a ferromagnetic state 
\Q, 0) over all angles of magnetization. Conversely, the mean-field state \Q, (p) may be interpreted 
as a broken symmetry state of the spin-singlet state with respect to the direction of quantization 
{9,4,). 

8.2. Many-body states of spin-2 BECs 

Next, we discuss the many-body states of spin-2 BECs in the single-mode approximation. 
Substituting (j/mif) - a„<p(r) in Eq. (l30l l. we obtain the spin-dependent part of the Hamiltonian 

as 

where 

f = Y,fmnal,a„ (410) 

mn 

is the spin vector operator, 

fz = Yj'nala,„ (411) 

m 

is its z-component, and 

= {s^y = 1 ^(-i)'"a,„a_,„ (4i2) 

m 

are the creation (S+) and annihilation {S-) operators of a spin-singlet pair 

While the operator >S+, when applied to a vacuum, creates a pair of bosons in the spin-singlet 
state, the pair should not be regarded as a single composite boson because S+ does not satisfy 
the commutation relations of bosons. In fact, the operators S± together with S, = (2N -H 5)/4 
satisfy the S U(l,l) commutation relations: 

[S„S^] = +S^, [S^,S-] = -24- (413) 

Here, the minus sign in the last equation is the only difference from the usual spin commutation 
relations. This difference, however, leads to some important consequences. In particular, the 
Casimir operator S^, which commutes with S± and S^, does not take the usual form of the 
squared sum of spin components but instead takes the form 

= -S+S- + Sj - S,. (414) 
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8.2.1. Eigenspectrum and eigenstates 

The eigenstates are classified according to the eigenvalues of the Casimir operator (14141) . 
Because S^S- = S\ — — is positive semidefinite, must have a minimum value ^S™'". 
Recalling that S, = {2N + 5)/4, >S™" can be expressed in terms of a non-negative integer A^o as 
iS™" - {2N() + 5)/4. Let |0) be the eigenvector corresponding to this minimum eigenvalue; then, 
S-\(p) = 0. It follows that 

S^\(p} = Sf^iSf" - 1)10); (415) 

hence, the eigenvalue of is given by v = 5(5' - 1) with S = (2A^o + 5)/4. 

The operation of S+ on |0) increases the eigenvalue of by 1 and that of .^V increases by 



2, as observed from the commutation relations ( 1413b . Therefore, the allowed combinations of 
eigenvalues S{S - 1) and S~ for and S., respectively, are given by 

S = (2A^o + 5)/4 (A^o = 0, 1,2, . . .) (416) 

and 

5, =5 +A^s (A^s =0, 1,2,...), (417) 

where A^o and A^s satisfy 

A^ = 2A^s+M), (418) 

and we may therefore interpret A^s as the number of spin-singlet "pairs" and A^o the number of 
remaining bosons. 



To find the exact energy eigenvalues of Hamiltonian ( 1409b . we first note that S± and 
commute with !F: 

[>S±,n =0, [S^f] =0. (419) 

The energy eigenstates can be classified according to quantum numbers A^o and A^s, total spin 
F, and magnetic quantum number F~. We denote the eigenstates as \No,Ns,F,F.;A}, where 
A = 1,2,..., gNo,F labels orthonormal degenerate states with gNo.F provided in Ref. [5[] . 

The first commutator in Eq. ( 1413b implies that S± plays the role of changing the eigenvalue 
of iS; by ±1 by creating or annihilating a spin-singlet pair. The eigensystem can be constructed 
by first preparing the eigenstates that do not involve spin-singlet pairs, 

\No, 0, F, F-.; A), (420) 

and then, by generating new eigenstates through successive operations of S on it: 

i«,«,,K,.,.,> = 45!^5^:M£iiL, (421, 



where || ||s -\J{i^\i//), and the orthnormality condition is 

{No, Ns, F, F,; A'\No, Ns, F, F,; A) = Sa'j- (422) 

To write down the eigenstates (1421b explicitly, we need to find the matrix elements of 5+, of 
which only the nonvanishing ones are 



(A^o,A^s + l,F,F,-A\S + \Nu,Ns,F,F^-A} = V(^s + l)(A^s + A^o + 5/2). 
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(423) 



Finally, the energy eigenvalue for the state [A'o, Ns, F, F,; A) is given by 



^ = ^[^(^ + 1) - + Jf^^sCA^ + A^o + 3) - pF„ (424) 

where the relationship 2A^s + Nq - N is used. 

The total spin F can, in general, take integer values in the range < F < 2No- However, there 
are some forbidden values [5]. That is, F - 1,2, 5, 2No- 1 are not allowed when A^o - 3k(k e Z), 
and F - 0, l,3,2A^o - 1 are forbidden when A^o = 3A; + 1. We prove this at the end of this 
subsection. 

To gain a physical insight into the nature of the energy eigenstates \N{),Ns, F, F~; A), it is 
useful to express them in terms of some building blocks. In fact, we can express them in terms of 
one-, two-, and three-boson creation operators. Let A^"*' be the operator such that when applied 
to the vacuum state, it creates n bosons that have a total spin F - f and magnetic quantum 
number F, - f. While the choice of a complete set of operators for the building blocks is not 
unique, we choose the following five operators for constructing the eigenstates; 



A^^'"'" = flz"'" (425) 

A|f'"'' = -^[(al/ -2a\al^+2alal^]^ (426) 

^(2)1- ^ _^[2^/2alal- ^/3{a]f] (427) 
Vl4 

Ajf " = --^[ V2(flo)^ - 3 V2fli t alaii + 3 V3(fl|)2fl^2 

-h3 y/3al(a[^f - 6 V^a^ajalj (^28) 

Af "'" = -^[(a"[)^- V6a^a"[aJ-H2(fl^)2fllj. (429) 



Note that Af^^ and A*^^^ vanish identically because of the Bose symmetry, and that the operators 
A^"*' commute with 

We next consider a set S of unnormalized states: 

l«12, «20, «22, «30, «33> ^ (fl|)"-(A[,2)^)«^o(4(2)t)«.z(4(3)ty,,o(^(3)1y,,3 ,^^0), 

(430) 

with n\2, "20, «22, «30 = 0, 1, 2, . . . , 00 and ^33 = 0, 1. The total number of bosons and the total 
spin of the state ( 14301 ) are - n\2 + 2(n20 + "22) + 3(n3o + ^33) and F - F-- 2(ni2 + «22) + 3^33, 
respectively. For given and F, «33 is uniquely determined by the parity of F: 

n33=Fmod2. (431) 

It can be shown jstl that S forms a nonorthogonal complete basis set of the subspace 'H(f_=f) in 
which the magnetic quantum number F^ is equal to the total spin F. It can also be shown [5] that 
the energy eigenstates can be represented as 

(iF_)A^(A[,'>'")"»P(w,=0)(a^)"'HAf^)"^HA|,'^^)"»(Af')"'^|vac), (432) 
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where «i2, nao, n22,«30 = 0, 1, 2, . . . , co, ^33 = 0, 1, and AF = 0, 1, . . . , 2F. These parameters are 
related to {No, Ns,F, F^} as 



A^o = "12 + 2«22 + 3n30 + 3n33, (43 3) 

A^s = "20, (434) 

F - 2«i2 + 2n22 + 3«33, (435) 

F, = F-AF, (436) 

and the corresponding eigenenergy is given by Eq. (14241 ). Note that the states defined in (I432I I 
are unnormahzed, and the states having the same energy (i.e., those belonging to the same set of 
parameter values {A^o, Ns, F, F-}) are nonorthogonal. 

It might be tempting to envisage a physical picture that the system, as in the case of "^He, to 
be made up of n„f composite bosons whose creation operator is given by A^'" . However, this is 

an oversimplification because the operator A^'' does not obey the boson commutation relation. 

Moreover, the projection operator P(Ns=i)) in ( 14321 ) imposes many-body spin correlations such 
that the spin correlation between any two bosons must have a vanishing spin-singlet component. 
Note that two bosons with independently fluctuating spins generally have a nonzero overlap with 
the spin-singlet state. The many-body spin correlations of the energy eigenstates are thus far 
more complicated than what a naive picture of composite bosons suggests. On the other hand, as 
long as quantities such as the number of bosons, magnetization, and energy are concerned, the 
above simplified picture is quite helpful. As an illustration, we provide an explanation for the 
existence of forbidden values for the total spin F. For example, to construct a state with F - 
or F - 3, composite particles with total spin 2 must be avoided, namely, nn - n22 = 0. Then, 
we have A^o = 3(n3o + n33), implying that F - or F - 3 is only possible when A^o - 3k(k e Z). 
For a state with F = 2 01 F - 5,we have nn -H n22 = 1 and A^o = 1 + "22 + 3(«3o + n33), implying 
that A^o ^ ^k{k € Z). The above simplified picture is also helpful when we consider the magnetic 
response discussed below. 



8.2.2. Magnetic response 

Here, we consider how the ground state and the magnetization F. respond to the applied 



magnetic field p. From Eq. (1424b . we see that the minimum energy states always satisfy F. - F 
when p > 0. The problem thus reduces to minimizing the function 



E(F„N,) 



[FziF, + 1) - 6A^] + -^Ns(2N - 2Ns + 3) - pF,. 



2yeff 



5yeff 



(437) 



Here, we only consider the case of C2 < 0. Detailed investigations of this problem can be 
found in Ref. ||5|]. In this case, it is convenient to introduce a new parameter 



C'l = Cl 



El 

20' 



(438) 



and consider the energy as a function of F, and / = 2A^o - F^: 

yeff 



E{F„r) 



2 yeff 

C2_ 

40y< 



y^" / C2 \ 1 

— V'^^}^-2 

-1(1 + 2F + 6) + const. 
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(439) 



While the averaged slope AF-/Ap ~ V /Cj coincides with that in MFT, the offset term C2/(SV ) 
in Eq. ( I439I I makes a qualitative distinction from MFT, that is, the onset of the magnetization 
displaces from p ^ lo p ^ Icil/i^'V^). Note that the slope V'^/c'^ and the offset Ical/CSV^'^) 
are determined by independent parameters. A typical behavior of the magnetic response when 
\c2\ » Cj is shown in Fig. [34] 




Figure 34: Typical dependence of the ground-state magnetization on the applied magnetic-field strength, for q < and 
c'j > 0. 

We now calculate the Zeeman-level populations of the ground states for C2 < 0. In MFT, the 
lowest-energy states for C2 < have vanishing population in the m = 0, +1 levels. In contrast, 
the exact ground states derived in the preceding subsection, (A^^^ ' (a!^)"'- (Aj^' ' )"^- (Af^ ' )"--' |vac) 
with n22 = 0, 1 and n^j, = 0, 1, have nonzero populations in the m -Q,+l levels. The exact forms 
for the averaged population {al„a,„) are calculated as follows. The above ground states have the 
form (Ay^")'*'^|0) oc (>S+)'^^|0), with \(p) being a state with a fixed number {s = n\2 + 2n22 + 3n33) 
of bosons satisfying S-\4>) - 0. The average Zeeman population for the ground states, 

\a„am) = 7 7. , (440) 

is then simply related to the average Zeeman populations for the state \<p) as 

(a^a„) = {a'lamh + ^ ^^ i{ala„,}o + <fli,„fl-,„)o + 1), (441) 

where {al„a„,}o = {(p\a^a,„\(p} / ((p\(p}. This formula implies that when A^s s> s, the Zeeman popu- 
lations of the ground states are sensitive to the form of 

With this formula, it is a straightforward task to calculate the average Zeeman-level popula- 
tions for four types of ground states, 

(Af'^'f'(al)"'HAf')"''-(Af')"''\vsic} (442) 
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with n22 = 0,1 and n^i = 0, 1 . A striking feature appears in the leading terms under the condition 
I ni2 ^ Ns- The results are summarized as 



<fl!fli> ~ (fl' fl_i> ~ A''s(l + n33)/«i2 



(443) 



and 



(a'flo) ~ A^s(l + 2«22)/ni2- 



(444) 



These results show that the population of each magnetic sublevel depends very sensitively on the 
spin correlations. As a consequence, a minor change in the magnetization might lead to a major 
change in the population. Such dramatic changes originate from bosonic stimulations caused by 
the term a^al^ in Aj^" or the term {02)^0' ^ in Af^' . 

8.2.3. Symmetry considerations on possible phases 

The possible ground states of spinor BECs can be inferred based on symmetry arguments 
without solving the GPEs. The spin-2 BEC has five internal degrees of freedom corresponding 
to magnetic quantum numbers m - 2, 1, 0, -1, -2. Let fl„,'s be their amplitudes. Under rotation 
of the coordinate system, they must transform so as to guarantee that the linear combination 



(445) 



m=-2 



is scalar, where Y^iO, (p) is the spherical harmonic function of rank 2 with 6 and (p being the polar 
and azimuthal angles, respectively, in the spherical coordinate system. Let 

k - {kjc,ky,k,) - (sin0cos(/>, sin0sin0, cos^) (446) 

be the unit vector in the (6, (f>) direction. In terms of the unit vector k, we can rewrite Eq. ( I445I I 
as 



d/= J—k^Mk, 



(447) 



where superscript T denotes the transpose, and 

fl2 + fl-2 - sJ^l^O 

M = r /(fl2-fl-2) 



1 



—ai + a_i 

lAxv Ipxy ^xz 
Ipxy fvv lAyz 
lA.c lAyz ^zz , 



i(a2 - fl-2) 

— fl2 — fl-2 — -^flO 
-/(fll + fl_i) 



— fll + fl_i 
— /(fll + fl_i) 

2 -^fl() 



(448) 



The order parameter is thus characterized by a 3 x 3 traceless symmetric matrix TrM - with 
unit normalization 



Tr( M*M) = 1. 



(449) 
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Let us now consider the eigenvalues of this order parameter matrix: 

det(£i- V3M) = £^-35£+ =0, (450) 

where V3 is added in front of the matrix for the sake of convenience, and S and T are the 
ampHtudes of the spin-singlet pair and spin-singlet trio, respectively: 

^ = = ^'Ao-'Ai'A-i +«A2«A-2, (451) 

'J 

= -9{i!j2iIj\ + tAitA-2) + V6.A()(-<Ao + 3<Ai<A-i + 61^21^-2). (452) 
Substituting S - xy and A - 2 ^/3(x^ + y^), Eq. ( 145 Oi l reduces to 

+y^ - 3Exy = 0, (453) 

and hence, the solutions are 

E - -(x + y), -(tox + o?y), -{cj?x + ojy), (454) 

where cj - e^"'^^ and 

x^-^-^(a+ Va2-4853)^ . (455) 
y 485 ^ ' 

We can infer from this that the spin-singlet pair and the spin-singlet trio of atoms form building 
blocks of the ground state of the spin-2 BEC. It is noteworthy that this fact arises solely from the 
fact that the order parameter matrix of the spin-2 BEC is traceless symmetric. 

9. Special Topics 

9.1. Quenched BEC: Kibble-Zurek mechanism 

Topological excitations can be spontaneously created when a system undergoes a phase tran- 
sition in a finite time. The basic concept is that after the phase transition domains of the new 
phase emerge at causally disconnected places, and therefore, the order parameters of the new 
phase at different locations are not correlated. When the order parameters grow to overlap, they 
may be able to adjust dynamically so that they are connected smoothly; if they cannot, singu- 
larities should develop in the order parameter space, giving rise to topolo gical excitations. Such 
a scenario of topological defect formation was first discussed by Kibble in the context of 
cosmic-string and monopole formations in the early Universe, and an experimental t est f or the 
Kibble scenario in condensed matter systems was proposed by Zurek [135] (see Ref. il36ll for a 
review). 

Bose-Einstein condensates of dilute atomic vapor offer an ideal opportunity for studying the 
Kibble-Zurek mechanism because the temperature, strength of interaction, and external param- 
eters such as magnetic field and trapping potentials can be changed in a time s horter than the 



characteristic time scale of topological defect formation. In fact, scalar vortices 111371 113811 and 
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spin vortices lllOOl Esll have been observed to emerge spontaneously upon both temperature and 
magnetic-field quenches. 

Topological defect formation via the Kibble-Zurek mechanism can occur in a ferromagnetic 
phase of a spinor BEC, where th e phase transition is triggered by a sudden change in a magnetic 
field (magnetic field quenching) ilOd 139 . 14d 141 1. The ground- state phase of a spin-1 ferro- 



magnetic BEC is shown in Fig. \T\ where p and q are the coefficients of the linear and quadratic 
Zeeman effects, respectively. When we consider a type of magnetic-field quench for which the 
total magnetization of the system is conserved, the quench occurs along a constant-p line of the 
phase diagram. As in the experiment described in Ref. [ 100], we consider the case of p - that 
implies that the total magnetization of the system is kept to be zero. Then, the state of the system 
can be in one of the following two phases. 

When q is larger than the critical value q'c = 2|ci|«, with n being the number density, the 
ground state is the polar phase ^ifP"'""^ = (0, ■\/n, 0) which is nonmagnetic. When q < q^, the 
ground state is the broken-axisymmetry (BA) phase in which the system develops transverse 
magnetization. The ground state for ^ < is the broken-antisymmetry (BA) phase with the 
order parameter given by 



(456) 



where q = q/qc- The BA phase has transverse magnetization = Fy + iF^. — e'" ^Jl - q^ 
which breaks the U(l) symmetry corresponding to the direction of the transverse magnetization. 
By rapidly quenching the magnetic field from above to below q^, local transverse magnetization 
develops in random directions, and thus, spin vortices are spontaneously created. 

To analyze such a situation, let us assume that the initial state is the polar phase (0, s/n, 0) 
and expand the field operator as 

<Am = ^ ( ^^'"-o + "-t.™) e''^*""*'^""'^ (457) 



k 



Substituting this into the Hamiltonian and keeping the terms up to the second order in ak,m with 
^ 9^ 0, we obtain the Bogoliubov Hamiltonian for the polar phase. We are interested in how the 
modes m-±l grow with time. The Heisenberg equations of motion for the modes are lll39ll 

d 

ih — ak±\ = (ek + q + c\n)ak±\ + c\na-k^\, (458) 
at ' ' ' 

with the solutions given by 

. / Ekt .ek + q + cm . Ekt\^ 
ak,±i(t) = I cos — - I — sm — \ak,±i (0) 

where 
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gives the Bogoliubov spectrum. If Ek has an imaginary part, the corresponding mode is dynam- 
ically unstable and it grows exponentially. Because c\ < and q > Q for spin-1 '^^Rb atoms, 
exponential growth occurs for ^ < ^c, in agreement with the phase boundary. 



9.7.7. Instantaneous quench 

We can use the solution ( 14591 ) to investigate how the transverse magnetization develops with 
time. Because i/'o - ^In, the magnetization operator F+-F'_-Fx + iFy is written as 



= V2;^(0'|(r) + iA-i('-)). 
It follows that the correlation function of the transverse magnetization is given by 
{F,(r,t)F.(r',t)} ^ ,tl^.!-*(-'). 



(461) 



(462) 



The exponential terms describe dynamical instabilities of which the dominant term is the one for 
which the imaginary part of Ek is maximal. Let be the "most unstable" wave number for 
which llmE^I is maximal. It follows from Eq. ( I460l l that = for q^/l < q < q^ and 



j2M (q, \ 



(463) 



for q < qc/2. Thus, for q < qc/2, the typical size of magnetic domains is ~ fc^^; for q^/l < q < 
q^, it is determined by the width of the distribution of ImEk- In any case, the size of the spin 
domains grows as the final value of q increases, as experimentally observed in Ref. jstl- The time 
scale T for the magnetization is given by 



T - 



which is h/ -\/4q(qc - q) for q^jl < q < q^ and hjq^ for q < <7c/2. 



(464) 



q{t) = qc\l 



9.7.2. Finite-time quench 

Next, we consider the case in which the magnetic field is quenched linearly in a finite time 

tq: 

(465) 

In this case, the Bogoliubov excitation energy ( 1460b also varies with time and the spin correlation 
function is estimated to be 

{F+{r, t)F4r', f)) oc J dkexp ^ J \Ek(t')\t'dt' + ik ■ (r - r') 

(X exp[/(f)-|r-rf/4], (466) 

where is the correlation length and /(f) determines the growth rate of magnetization. In two 
dimensions, we have B 13911 



16^^ 

M2 



■f(TQ - t) 
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(467) 



and 



tan ' 



(468) 



The time scale for magnetization is determined by f(t) - 1 . Solving this equation by assuming 
that f <K tq, we obtain 

tQ-(^Y. (469) 



Substituting this into Eq. (1467b . we obtain 



4 



The same power law was obtained in Ref. 014011 . 

9.1.3. Numerical simulations 

The topological defect formation caused by a quantum phase transition can be best illustrated 
by numerical simulations. Figure[35]shows the time evolution following an instantaneous quench 
of the magnetic field from q > to q - Q (Fig. lSSl a)) and q - q^jl. Here, we assume that the 
system is a two-dimensional disk with a hard wall at radius 100 /vm and the potential is assumed 
to be flat inside the wall. The initial state is considered to be a stationary state of the GP equations 
and the following quench dynamics is obtained by using the time-dependent GP equations. To 
trigger the dynamical instabilities that cause defect formation, we assume small fluctuations in 
the initial amplitudes of the m-+ \ components according to 

flt.±i(0) = ff + i^, (471) 

where a and fi are random variables whose amplitudes x follow the normal distribution p(x) - 

Figure [35] shows snapshots of the amplitude |F+(r)| and the phase argF+(r) of the transverse 
magnetization after the quench. We observe that many holes emerge spontaneously after 100 ms. 
These holes represent topological excitations called polar-core vortices in which the m = ±1 
components have vortices with the cores filled by the m - components. 

By counting the number of vortices in the final state, we can determine the number of spin 
vortices. Alternatively, we can calculate the same number — the spin winding number w — from 
the algebraic relation 

where C{R) is a circle of radius R and w is a count of the number of rotations of the spin vector 
in the xy plane along C{R). Figure[36]shows the R dependence of the ensemble average of w^(R). 
We note that (w^(7?))avg is proportional to R for large R in agreement with the KZ theory, whereas 
it is proportional to R^ for small R. 
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Figure 35: Spontaneous magnetization following the quench from q > q^\o q < q^. Shown are time developments of 
the magnetization (upper) and its direction arg |F+| (lower) for (a) ^ = and (b) q = The black circles in (b) 

indicate a topological defect. Reprinted fromRef. H39il . 

If spin vortices can be generated randomly, (w^(/?))avg should be proportional to the area of 
the system and hence to R^. The KZ scaling law oc R that we observe for large R arises from 
the fact that the number of vortices matches that of anti-vortices inside the disk due to spin 
conservation. In other words, for small q, magnetic domains are aligned in such a manner as to 
cancel the local spin when averaged over the spin correlation length. Thus, over a greater length 
scale, the magnetic domains can be created independently. This is the underlying physics that 
makes the spin conservation compatible with the KZ postulate of independent defect creation at 
a long distance and the KZ scaling law valid for 2D systems. 



9.2. Low-dimensional systems 

It had been conjectured from the work of Peierls and Landau in the 1930s 11 14211 14311 that 
two-dimensional systems with continuous order parameters cou ld no t have a conventional long- 
range order, and this was confirmed by Mermin and Wagner il44ll and by Hohenberg 1 145]. 
Long wavelength fluctuations are more important in one and two dimensions than in three di- 
mensions and deviations from the mean-field theory are much greater. However, there still exists 
a type of quasi-long-range order in the low-temperature phase in two-dimensional systems that 
is characterized by a power-law decay of the order parameter correlation function. Such a phase 
transition is known as the Berezinskii-Kosterlitz-Thouless (BKT) transition. In this subsection, 
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Figure 36: R dependence of the variance of the spin winding number w(/?) for instantaneous quench of the magnetic 
field too = and to ^ = q^jl. The dashed and dotted lines are respectively proportional to R and i?^. Reprinted from 
Ref. {vm. 



we first explain the BKT transition for a single-component Bose gas, and then, we extend it to a 
spinor gas. 

9.2.7. Berezinskii-Kosterlitz-Thouless transition in a single-component Bose gas 
We start from a Ginzburg-Landau energy functional 



(473) 



which gives an equilibrium probability for iff as Peq[i//] oc exp[-j^]. Here,y6 > and a changes 
its sign from positive to negative as the temperature decreases from above to below the critical 
temperature T^, so that the equilibrium value of tfr changes from i/^ = at T > to = t/^o = 
VnZs at r < T,. 

In three dimensions, the one-body correlation function (i/'*(r)i^(0)) remains finite at large 
distance below T^, exhibiting a long-range order. In constrast, above Tc, the correlation function 
decays exponentially: 

<iA*(r)i/r(0)) oc for r ^ oo, (474) 



where ^+ - ypTWl oc (T - T^Y^^^ except in the immediate vicinity of Tc (see e.g. Ref. lll46ll ). 
In this case, the order is short-ranged. 

In two dimensions, the effect of fluctuations is much more drastic. We first assume that 
|t/r| takes its equilibrium value everywhere and consider only phase fluctuations. Substituting 
tf/{r) = i/roe''^*''> in Eq. (l473l l. we have 

= ^ J /:o|V0| Vr = i/To kVk'Pk, (475) 
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where = -/oliAol and 4>k is the Fourier component of 4>{r) satisfying 



Because the 



gradient of induces a superfluid current as Vs = ih/M)'V(f>, the coefficient Kq is related to the 
superfluid density as 



M^Ko 



so that H is rewritten as 

Because Eq. ( 14771 ) is quadratic in we can calculate the correlation function as 



liAol exp 



■Z 



|1 - cos(^ ■ r)] 



In a two-dimensional system, the summation of k is replaced by an integral as 



1 k T C dk C 
{1 - cos(A: ■ r)) = — - cos(^rcos0)] 

(2ny Kq J k Jo 

J dk- 



ksT_ 
IttKq 



Jo(kr) 



(476) 

(477) 

(478) 
(479) 

(480) 
(481) 



where J()(z) is the zeroth Bessel function. For k > 1 /r, we can ignore the rapidly oscillating 
function Joikr), whereas the contribution from k < 1/r is small. We also introduce a cutoff 
wavelength ^ corresponding to the coherence length, and neglect spatial structures smaller than 
^. Then, the correlation function is calculated as 



<r('-)^(0)) = |-Aorexp 



far 

2kKo J^-i 



dk 
T 



where 



ksT 
2nKo' 



(482) 



(483) 



Thus, phase fluctuations destroy the true long-range order; however, the quasi long-range order 
characterized by a power-law behavior remains. 

Next, we consider the effect of fluctuations on the amplitude A crucial effect arises when 
\t//\ passes through zero. If \^(r)\ is finite everywhere in a singly connected region, the phase (p of 
the order parameter is defined as a single-valued function. On the other hand, if i/' = at r = r,, 
can be a multi-valued function. The points {r,) correspond to vortices around which (p varies 
by an integral multiple of 2n. 

The phase transition for the appearance of a quasi-long-range order can be captured by a 
simple discussion of vortex nucleation. When a single vortex passes between two distinct points, 
the relative phase between these two points changes by +27r, implying that free vortices destroy 
the phase coherence. Hence, the quasi-long-range order is destroyed when free vortices are 
thermally excited. 

Ill 



The free energy for a single vortex can be calculated as follows. The circulation of the 
superfluid velocity around the vortex is quantized as ■ d( — h/M. Assuming axisymmetry 
around the vortex, we have Vs - fve^HMr), where (r, ip) is the polar coordinate around the vortex 
and e^p is a unit vector in the direction of (f. Then, the energy cost required to create a single 
vortex is given by 

£v = £c + ^J rdr #P°KP (484) 

= £e+^^olnJ|j, (485) 

where R is the radius of the system; ^, the radius of the vortex core; and E^, the vortex core 
energy arising from the region of r < ^. On the other hand, because the vortex can freely move 
in the area of R^, the entropy is given by 5 = ln(^^/^^)- Then, we obtain the free energy of a 
free vortex as 

F^Ey-TS ^ (nKo - Ik^T) In || j + E,. (486) 
One can immediately see from this free energy form that the phase transition occurs at 

Tco = f^, (487) 

above which free vortices are generated and destroy the phase coherence. From Eqs. ( I483I I and 
(1487b . the exponent rj must satisfy tj < 1/4. 

A more careful argument requires the interaction between vortices to be considered. If there 
exists a pair of a vortex (charge 1 ) and an anti-vortex (charge - 1 ) at distance d, the velocity fields 
generated by these vortices at a distance larger than d cancel each other out. As a result, the 
vortex pair has a finite energy of 



d 



Ep^i^id) = 2£c + 2nKo In - , (488) 



which indicates that a vortex and an anti-vortex attract each other and form a pair at low temper- 
ature. Because vortex pairs do not destroy phase coherence (see Fig. [37l i. the superfluid phase 
transition can be interpreted as a dissociation of vortex pairs. 

If there exist many vortex pairs, the attractive interaction between a vortex and an anti-vortex 
is weakened by the screening effect. Kosterlitz and Thouless 1. 147.] analyzed this phase transition 
by mapping the system to a two-dimensional Coulomb gas of unit charge q = V^r^. They 
introduced a "dielectric constant" s that renormalizes K^) to K = Kofe, and constructed the 
renormalization-group flow equations. Their main result is that the superfluid density jumps at 
the transition temperature that satisfies the universal relation: 



This universal behavior in superfluid helium was observed by Bishop and Reppy 1114811 . With 
regard t o col d atom systems, the BKT transition of a single-component Bose gas was observed 
in Ref . im 
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Figure 37: Local phase of the order parameter around (a) a single vortex and (b) a vortex pair, where the direction of 
each aiTow indicates a phase that varies between and 2n. The phase around a single vortex varies, causing a velocity 
field proportional to 1/r, whereas that away from the vortex pair can be uniform. 



9.2.2. 2D spinor gases 

We now consider the case of a spinor. As discussed above, vortices play a crucial role in the 
superfluid phase transition in two-dimensional systems. When a system acquires spin degrees of 
freedom, the BKT phase transition becomes highly nontrivial because the nature of the vortices 
depends on the order parameter manifold, as discussed in Sec. [7] It is well known that though 
the XY spin model exhibits a BKT phase transition where the order parameter manifold U(l) 
is the same as that of a single-component R ose g as, the BKT transition does not occur for the 



Heisenberg spin model at finite temperature 1114711 . This is because the order parameter manifold 



for the Hisenberg model, which is isomorphic to S^, does not accommodate a topologically 
stable vortex because niiS^) = 0. As discussed in Sec. |7] the order parameter manifold for the 
polar phase of a spin-1 EEC can accommodate half-quantum vortices, and therefore, the BKT 
transition occurs due to the binding of half-quantum vortex pairs. On the other hand, if the 
condensate accommodates several types of vortices, the BKT transition is expected to occur in 
several stages. Here, we discuss two BKT transitions of ferromagnetic spinor gases that occur 
under an external magnetic field. 

For simplicity, we neglect the trapping potential and assume an infinite system. The energy 
functional for a spin-1 spinor gas under an external magnetic field B in the z-direction is given 
by 



{ d\ ^ Yj |V^„,|' + co«' + ci|f|2 + ^(|^i|2 + |^_,|2) 

^^^^ m=0,±l 



(490) 



where n - 2„, li/'ml and P - Tjmn ^m(f)mn4^n are the number density and spin density, re- 
spectively, and q is the quadratic Zeeman energy proportional to B^. For spin-1 BECs, q - 
(//BB)^/(4Ahf). The linear Zeeman term simply causes Larmor precession due to the spin con- 
servation, and therefore, it is omitted by going onto a rotating frame of reference in the spin 
space. 

We consider the ferromagnetic interaction, ci < 0. At zero temperature, the ground state 
phase is the same as that we have obtained for a three-dimensional system in Sec. [3] Here, 
we assume that the total longitudinal magnetization is conserved to be zero. Then, the order 
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parameter for q > = 2|ci |no is the polar phase whose order parameter is given by 



^1 



polar 



^Jr^)e"• 



(0\ 
1 . 



(491) 



whereas for q < q^, the broken-axisymmetry phase appears: 



BA 



V2(l+§) 



(492) 



where «o is the absolute square of the order parameter at zero temperature and q = q/q^- 

For q > qc, there exists only one type of vortex; a conventional phase vortex. Therefore, the 
BKT transition for <7 > <7c is the same as that in single-component gases. Unbinding of phase 
vortex pairs destroys the coherence of the U(l) phase that appears in the correlation function of 
the m = component: 



<tAS(r)-Ao(0)) r-\ 



(493) 



where rf is given by Eq. ( 14831 ) with renormalized K instead of Kq - h^n^)/M. The universal jump 
is given by Eq. ( 14891 ). 

On the other hand, for q < q^, there exist two types of vortices, i.e., a phase vortex and 
a spin vortex. The phase vortex is the same as that in the polar state, around which the U(l) 
phase (p changes by integer multiples of In. The spin vortex is defined as a vortex around which 
the direction of transverse magnetization (i.e., a) varies by integer multiple of 27!-. Here, the 
transverse magnetization appears atqKq^ and it is given by 



l<A- 



0. 



(494) 
(495) 



The effects of these vortices are independent because the energy for these fluctuations are decou- 
pled as 



where 



^0 = 



sp _ H^rioil - q) 



^0 = 



2M 



(496) 

(497) 
(498) 



Similar to the discussion for Eq. (14861 ), a rough estimation for the transition temperature for the 
appearance of a free phase vortex and a free spin vortex are given by 



cO 



^ cO 



2^B " 
2^B 



no 



nffrioil -q) _ 1 -§^p 
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respectively. Hence, as the temperature decreases, phase vortex pairs first bind and then spin vor- 
tices form pairs at lower temperature. More proper treatment following the Kosterlitz-Thouless's 
renormalization method gives the relation between the universal jump and the transition temper- 
ature: 

The phase transitions related to the binding of the phase and spin vortices can be observed 
as an appearance of the quasi-long-range order in the correlation function of m = and m - +\ 
components, respectively, which are given by 

{il/l{r)il,Qm oc (g-'f-AW-^W)) oc r'^^, (502) 

<^/'li(r)iA±i(0)> oc (e-'W'-)-0(O)))(e^''W'-)-«'(O))) oc r'^v-ri^^ (503) 

with 

'^P^^- '7sp = ^- (504) 
The binding of spin vortices also contributes to the ferromagnetic long-range order 

{F4r)F-m oc oc r"''*. (505) 

Figure [38] summarizes the above results. Depending on the temperature and the external 
magnetic field, there exist three ordered phases: (I) quasi-long-range order arises only in the 
m - component, where the gas has no local magnetization; (II) a quasi-long-range order arises 
only in the m - component, where the gas is locally magnetized, but with no long-range 
order of magnetization; (III) a quasi-long-range order appears in every m component, where the 
transverse magnetization also exhibits the quasi-long-range order. 

Here, we discuss the phase diagram at low magnetic fields. In the above system, spin vortices 
are stable because the quadratic Zeeman effect suppresses the fluctuation of longitudinal magne- 
tization. As the external field decreases, however, the contribution of such fluctuations becomes 
larger and induces coupling between the spin and the phase vortices. Then, both transition tem- 
peratures are expected to be suppressed. At zero magnetic field, the quasi-long-range order does 
not appear at finite temperature as in the case of the Hisenberg spin system. Understanding the 
detailed behavior in the low-magnetic-field region requkes further investigation. 



10. Summary and Future Prospects 

In the present paper, we have reviewed the basic knowledge concerning spinor Bose-Einstein 
condensates (BECs) that has been accumulated thus far. The fundamental characteristics of 
spinor BECs are rotational invariance, coupling between the spin and the gauge degrees of free- 
dom, and magnetism arising from the magnetic moment of the spin. 

The rotational invariance and gauge invariance alone uniquely determine the microscopic 
Hamiltonian of the spinor BEC, as discussed in Sec. |2l The mean-field theory described in 
Sec. [3]discribes the ground-state phases and the dynamics of spinor BECs. In a scalar BEC, 
the measurements of the collective modes are found to agree with the theoretical values with 
an accuracy better than 1% (see Ref. IJ50il for a comprehensive review). In the near future, 
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Figure 38: Schematic pliase diagram of a 2D ferromagnetic Bose gas. (I) a quasi-long-range order appears only in 
the m = component with no local magnetization; (II) a quasi-long-range order aiises only in the m = component, 
where the gas is locally magnetized, but without any long-range order of the magnetization; (III) a quasi-long-range 
order appears in every m component, where the transverse magnetization also exhibits a quasi-long-range order In the 
shaded region, the fluctuation of the longitudinal magnetization is expected to destroy the quasi-long-range order 

experimental tasks to be carried out include investigating the extent to which the predictions of 
the Bogoliubov theory described in Sec.|5]can be verified experimentally. 

The topological excitations described in Sec. Q are of cross-disciplinary interest; in particu- 
lar, the non-Abelian nature of vortices in the nematic and cyclic phases of spin-2 BECs will have 
fundamental implications on quantum turbulence and other subfields of physics. The topological 
excitations are usually considered within a given order parameter manifold. However, as exem- 
plified in Fig.[33]in Sec. 17.3.41 the order-parameter manifold can change dynamically from one to 
another, and during a transition, the topological charge is no longer conserved. It is an interesting 
theoretical question to investigate whether, if any, a more general topologically invariant quantity 
can be defined in such a general situation. The relative homotopy theory described in Sec. 17.2.21 
may be used as a guiding theoretical framework. 

A closely related subject is vortex nucleation dynamics and vortex lattice formation. Vortices 
are singularities of the order parameter In a scalar EEC, the only possible singularity is the 
density hole in the condensate. However, the spinor EEC has several different phases and the 
choice of creating a density hole is energetically most costly. Thus, in spinor EECs, the vortex 
core is filled by a state that does not belong to the order parameter manifold of the condensate. 
As discussed in Sec.|6l the spin-1 polar EEC has a ferromagnetic vortex core, and the spin-1 
ferromagnetic EEC can host a polar-core vortex. For spin-2 and higher-spin EECs, more than two 
possible phases can fill the vortex core. Eecause each phase has its own internal spin structure, 
an interesting question arises as to ho w the order parameter of the EEC smoothly connects the 



The many -body quantum correlations of spinor EECs described in Sec. [8] are fragile against 
an external magnetic field because the bosonic stimulation favors the mean field. The spinor- 
dipolar phenomena described in Sec. 14.31 are also vulnerable against a magnetic field because 
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the Larmor precession averages out the part of the dipole interaction that couples the different 
total spin components. It is noteworthy that both many-body spin correlations and spinor-dipolar 
effects become significant as the external magnetic field is reduced down to approximately 10 ;uG. 
We have good reasons to expect that further symmetry breaking occurs in such an environment 
of ultralow magnetic field. 

From the viewpoint of fundamental physics, the topological defect formation based on the 
Kibble-Zurek mechanism discussed in Sec. l9.1l merits further experimental and theoretical study. 
In particular, this effect combined with the spinor Berezinskii-Kosterlitz-Thouless transition in 
two dimensions described in Sec. 19.2. f] raises many interesting questions. In fact, the gauge, 
spin, and dipole-dipole interactions set very different energy scales that, we can conceive, lead 
to topological phase transitions at different energy and temperature scales. 

All the discussions presented in this review are restricted to absolute zero. However, experi- 
ments on spinor BECs have been carried out at temperatures much higher than the typical energy 
scale of the spin-exchange interaction. Nevertheless, once the system undergoes Bose-Einstein 
condensation, all the condensed particles respond to an external perturbation in exactly the same 
manner. When the number of condensate particles is of the order of 10^, the collective energy 
for the dipole interaction amounts to 1 /iK for '^^Rb and 100 /iK for ^^Cr Therefore, for temper- 
atures below 0. 1 fjK, the dynamics of the system can be understood, at least qualitatively, based 
on the zero-temperature theory. Nonetheless, it is likely that the precise measurements of collec- 
tive modes are sensitive to finite-temperature effects. Moreover, when the long-time dynamics 
is considered, the thermaUzation process is expected to play a crucial role in determining the 
magnetic structure of the system. In fact, very recent experiments on the '^^Rb / = 1 BEC show 
evidence of magnetic ordering |25] as a consequence of thermalization. To fully understand 
such phenomena and to explore new types of temperature-quenched symmetry-breaking phase 
transitions, it is obvious that we need to develop theoretical tools to investigate the dynamics of 
spinor-dipolar physics at finite temperatures. 
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